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Scilab numbering policy used in this document and the relation to the 
above book. 


Scilab numbering policy used in this document and the relation to 
the above book. Example (Solved example) Example (Solved ex- 
ample) Equation (Particular equation of the above book) Equation 
(Particular equation of the above book) Appendix to Example (Scilab 
Code that is an Appednix to a particular Example of the above 
book) 

For example, Exa 3.51 means solved example 3.51 of this book. Sec 2.3 
means a scilab code whose theory is explained in Section 2.3 of the 
book. 
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Chapter 2 

Signals and Systems 


Scilab code Exa 1.1 Time Shifting 

//clear// 

//Example 1.1: Time Shifting 

//SIGNALS & SYSTEMS, Second Edition 

/ /V . OPPENHEIM , S . WILLSKY , S . HAMID NAMWAB 

//PHI, 2008 Edition 

//Page 10 

clear ; 

clc ; 

close ; 

t = 0 : 1/100 : 1 ; 
for i = l:length(t) 
x ( i ) = 1 ; 

end 

for i = length (t) +1 : 2* length (t) 
x(i) = 1-t (i-length (t) ) ; 
end 

tl = 0 : 1/100 : 2 ; 
t2 = -1:1/100:1; 

// 1 3 = 0:1/100:4/3; 
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/ / 1 4 = 0:l/length(t3):l; 

//Mid =ceil(length(t3)/2); 

//for i = 1 : Mid 
// x3(i) = 1 ; 

/ / end 

//for i = Mid + 1: length ( t3 ) 

// x3 ( i ) = 1 — 1 4 ( i —Mid ) ; 

/ / end 

f igure 
a=gca () ; 

plot2d ( t 1 ,x(l:$-l)) 

a . thickness =2 ; 

xtitle ( ’ The signal x(t)’) 

f igure 

a=gca () ; 

plot2d(t2 j x ( 1 : $-1) ) 
a . thickness =2 ; 
a.y_location = ’’middle”; 
xtitle (’The signal x(t+l)’) 
f igure 
a=gca () ; 

plot2d(t2 ,x($ : -1 :2) ) 
a . thickness =2 ; 
a.y_location = ’’middle”; 
xtitle (’The signal x(-t+l)’) 

Scilab code Exa 1.1 

Time Shifting 

//clear// 

//Example 1.1: Time Shifting 

//SIGNALS & SYSTEMS, Second Edition 

/ /V . OPPENHEIM , S . WILLSKY , S . HAMID NAMWAB 

//PHI, 2008 Edition 

//Page 10 

clear ; 

clc ; 
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close ; 

t = 0 : 1/100 : 1 ; 
for i = l:length(t) 
x ( i ) = 1 ; 

end 

for i = length (t) +1 : 2* length (t) 
x(i) = 1-t (i-length (t) ) ; 
end 

tl = 0 : 1/100 : 2 ; 
t2 = -1:1/100:1; 

// 1 3 = 0:1/100:4/3; 

/ / 1 4 = 0:l/length(t3):l; 

//Mid =ceil (length ( 1 3 ) / 2 ) ; 

//for i = 1 : Mid 
// x3 ( i ) = 1 ; 

/ / end 

//for i = Mid + 1: length ( t3 ) 

// x3(i) = 1 — 1 4 ( i —Mid ) ; 

/ / end 

f igure 
a=gca () ; 

plot2d(tl ,x(l:$-l)) 

a . thickness =2 ; 

xtitle ( ’ The signal x(t)’) 

f igure 

a=gca () ; 

plot2d(t2 ,x(l : $-1) ) 
a . thickness =2 ; 
a.y_location = ” middle”; 
xtitle (’The signal x(t+l)’) 
f igure 
a=gca () ; 

plot2d(t2 ,x($ : -1 :2) ) 
a . thickness =2 ; 
a.y_location = ’’middle”; 
xtitle (’The signal x(-t+l)’) 
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Scilab code Exa 1.2 Time Scaling 

//clear// 

//Example 1.2: Time Scaling 

//SIGNALS & SYSTEMS, Second Edition 

/ /V . OPPENHEIM , S . WILLSKY , S . HAMID NAMWAB 

//PHI, 2008 Edition 

//Page 11 

clear ; 

clc ; 

close ; 

t3 = 0 : 1/100:4/3; 
t4 = 0 : 1/ length (t3) : 1 ; 

Mid =ceil (length (t3) /2) ; 
for i = 1 : Mid 
x3(i) = 1 ; 
end 

for i = Mid+1 : length (t3) 
x3(i) = l-t4 (i-Mid) ; 
end 

f igure 
a=gca () ; 
plot2d (t3 , x3 ) 
a . thickness =2 ; 

x t i 1 1 e ( ’ Time Scaling x(3t/2)’) 

Scilab code Exa 1.2 

Time Scaling 

//clear// 

//Example 1.2: Time Scaling 
//SIGNALS & SYSTEMS, Second Edition 
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/ / V . OPPENHEIM , S . WILLSKY , S . HAMID NAMWAB 
//PHI, 2008 Edition 
//Page 11 

clear ; 
clc ; 
close ; 

t3 = 0: 1/100:4/3; 
t4 = 0 : 1/ length (t3) : 1 ; 

Mid =ceil (length (t3) /2) ; 
for i = 1 : Mid 
x3(i) = 1 ; 
end 

for i = Mid+1 : length (t3) 
x3(i) = l-t4 (i-Mid) ; 
end 

f igure 
a = gca () ; 
plot2d (t3 , x3 ) 
a . thickness =2 ; 

x t i 1 1 e ( ’ Time Scaling x(3t/2)’) 


Scilab code Exa 1.3 Time Scaling and Time Shifting 

//clear// 

//Example 1.3: Time Scaling and Time Shifting 

//SIGNALS & SYSTEMS, Second Edition 

/ /V . OPPENHEIM , S . WILLSKY , S . HAMID NAMWAB 

//PHI, 2008 Edition 

//Page 11 

clear ; 

clc ; 

close ; 

t3 = 0: 1/100:4/3; 
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t4 = 0 : 1/ length (t3) : 1 ; 

Mid =ceil (length (t3) /2) ; 
for i = 1 : Mid 
x3(i) = 1 ; 
end 

for i = Mid+1 : length (t3) 
x3(i) = l-t4 (i-Mid) ; 
end 

t5 = -2/3:1/100:2/3; 
f igure 
a=gca () ; 
plot2d (t5 , x3 ) 
a . thickness =2 ; 
a.y_location =” middle”; 

xtitle(’Time Scaling and Time Shifting x ( ( 3 t /2 ) +1) ’ ) 

Scilab code Exa 1.3 

Time Scaling and Time Shifting 

//clear// 

//Example 1.3: Time Scaling and Time Shifting 

//SIGNALS & SYSTEMS, Second Edition 

/ /V . OPPENHEIM , S . WILLSKY , S . HAMID NAMWAB 

//PHI, 2008 Edition 

//Page 11 

clear ; 

clc ; 

close ; 

t3 = 0: 1/100:4/3; 
t4 = 0 : 1/ length (t3) : 1 ; 

Mid =ceil (length (t3) /2) ; 
for i = 1 : Mid 
x3(i) = 1 ; 
end 

for i = Mid+1 : length (t3) 
x3(i) = l-t4 (i-Mid) ; 
end 
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t5 = -2/3:1/100:2/3; 
f igure 
a=gca () ; 
plot2d (t5 , x3 ) 
a . thickness =2 ; 
a.y_location =” middle”; 

xtitle(’Time Scaling and Time Shifting x ( (3 t /2) +1) ’ ) 


Scilab code Exa 1.4 Combinationation two periodic signals 

//clear// 

//Example 1 . 4 : Combinationation two periodic signals 
// Aperiodic signal 
//Page 12 

clear ; 
clc ; 
close ; 

F = 1 ; //Frequency = 1 Hz 

tl = 0: -1/100: -2*°/„pi ; 
xl = cos (F*t 1 ) ; 
t2 = 0 : 1/100 : 2*%pi ; 
x2 = sin (F*t2) ; 
a=gca () ; 
plot (t2 , x2 ) ; 
plot (tl , xl ) ; 

a.y_location = ’’middle”; 
a.x_location = ’’middle”; 

xtitle(’The signal x(t) = cost for t < 0 and sint 
for t > 0: Aperiodic Signal’) 

Scilab code Exa 1.4 

Combinationation two periodic signals 
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//clear// 

//Example 1 . 4 : Combinationation two periodic signals 
// Aperiodic signal 
//Page 12 

clear ; 
clc ; 
close ; 

F = 1 ; //Frequency = 1 Hz 

tl = 0: -1/100: - 2 * %p i ; 
xl = cos (F*t 1 ) ; 
t2 = 0 : 1/100 : 2*%pi ; 
x2 = sin (F*t2) ; 
a=gca () ; 
plot (t2 , x2 ) ; 
plot (tl , xl ) ; 

a.y_location = ’’middle”; 
a.x.location = ’’middle”; 

xtitle(’The signal x(t) = cost for t < 0 and sint 
for t > 0: Aperiodic Signal’) 


Scilab code Exa 1.6 Fundamental period of composite discrete time sig- 
nal 

//clear// 

//Example 1 . 6 : Determine the fundamental period of 
composite 

// discrete time signal 

/ / x [ n ] = exp ( j (2* %pi /3) n)+exp ( j (3*%pi/4)n) 

clear ; 
clc ; 
close ; 

Oraegal = 2*°/„pi/3; //Angular frequency signal 1 
0mega2 = 3*°/„pi/4; //Angular frequency signal 2 
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N1 = (2* Zpi ) / Omegal ; //Peirod of signal 1 
N2 = (2* Zpi ) / 0mega2 ; //Period of signal 2 
//To find rational period of signal 1 

for ml = 1 : 100 
period = N1 *ml ; 
if (modulo (period , 1) ==0) 
periodl = period; 
int eger_value = ml 
break ; 
end 
end 

//To find rational period of signal 2 

for m2 = 1 : 100 
period = N2*m2; 
if (modulo (period , 1) ==0) 
period2 = period; 
int eger_value = m2 
break ; 
end 
end 

disp (periodl ) 
disp (period2 ) 

//To determine the fundamental period N 

N = periodl *period2 

Scilab code Exa 1.6 

Fundamental period of composite discrete time signal 

//clear// 

//Example 1 . 6 : Determine the fundamental period of 
composite 

// discrete time signal 

//x [ n ] = exp ( j (2* %pi /3) n)+exp ( j (3*%pi/4)n) 

clear ; 
clc ; 
close ; 

Omegal = 2*°/„pi/3; //Angular frequency signal 1 
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0raega2 = 3*°/„pi/4; //Angular frequency signal 2 
N1 = (2* °/ 0 pi ) / Omegal ; //Peirod of signal 1 
N2 = (2* °/ 0 pi ) / 0mega2 ; //Period of signal 2 
//To find rational period of signal 1 

for ml = 1 : 100 
period = N1 *ml ; 
if (modulo (period , 1) ==0) 
periodl = period; 
int eger_value = ml 
break ; 
end 
end 

//To find rational period of signal 2 

for m2 = 1 : 100 
period = N2*m2; 
if (modulo (period , 1) ==0) 
period2 = period; 
int eger_value = m2 
break ; 
end 
end 

disp (periodl ) 
disp (period2 ) 

//To determine the fundamental period N 

N = periodl *period2 


Scilab code Exa 1.12 Classification of system 

//clear// 

//Example 1 . 1 2 : C 1 as s i f i c a t i o n of system : Causality 
property 
//Page 47 

//To check whether the given discrete system is a 
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Causal System (or) Non— Causal System 
//Given discrete system y[n]= x[ — n] 

clear ; 
clc ; 

x = [2,4,6,8,10,0,0,0,1]; //Assign some value to 
input 

n = -length (x) /2 : length (x) /2 ; 
count = 0; 

raid = ceil ( length (x) /2) ; 
y = zeros ( 1 , length (x) ) ; 
y(mid+l:$) = x ($ : -1 : mid+1) ; 
for n = — 1 : — 1 : -mid 
y(n+l+mid) = x(-n); 
end 

for i = l:length(x) 
if (y ( i ) ==x ( i ) ) 
count = count + 1 ; 
end 
end 

if ( count == length (x) ) 

dispC’The given system is a causal system’) 

else 

dispC’ Since it depends on future input value’) 
dispC’The given system is a non— causal system’) 

end 

Scilab code Exa 1.12 

Classification of system 

//clear// 

//Example 1 . 1 2 : C 1 as s i f i c a t i o n of system : Causality 
property 
//Page 47 

//To check whether the given discrete system is a 
Causal System (or) Non— Causal System 
//Given discrete system y[n]= x[ — n] 

clear ; 
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clc ; 

x = [2,4,6,8,10,0,0,0,1]; //Assign some value to 
input 

n = -length (x) /2 : length (x) /2 ; 
count = 0; 

raid = ceil ( length (x) /2) ; 
y = zeros ( 1 , length (x) ) ; 
y(mid+l:$) = x ($ : -1 : mid+1) ; 
for n = — 1 : — 1 : -mid 
y(n+l+mid) = x(-n); 
end 

for i = l:length(x) 
if (y ( i ) ==x ( i ) ) 
count = count + 1 ; 
end 
end 

if ( count == length (x) ) 

dispC’The given system is a causal system’) 

else 

dispC’ Since it depends on future input value’) 
dispC’The given system is a non— causal system’) 

end 


Scilab code Exa 01.13 Determination of stability of a given system 

//clear// 

//Example 1 . 1 3 ( b ) : Deter rninat ion of stability of a 
given system 
//Page 50 

//given system y(t) = exp(x(t)) 

clear ; 
clc ; 

Maximum_Lirait = 10; 
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S = 0; 

for t = 0 : Maximum_Limit - 
x (t+1) = -2~t ; 

S = S+exp (x (t+1) ) ; 
end 

if (S >Maximum_Limit ) 
disp ( ’ Eventhough input 
unbounded ’ ) 

disp (’The given system 
disp ( ’ S = ’ ) ; 

S 

else 

disp(’The given system 
disp (S) ; 
end 


//Input some bounded value 

is bounded output is 
is unstable’); 

is stable ’) ; 


Scilab code Exa 01.13 

Determination of stability of a given system 

//clear// 

//Example 1 . 1 3 ( b ) : Deter rninat ion of stability of a 
given system 
//Page 50 

//given system y(t) = exp(x(t)) 

clear ; 
clc ; 

Maximum_Lirait = 10; 

S = 0; 

for t = 0 : Maximum_Lirait -1 

x(t+l)= -2~t; //Input some bounded value 

S = S+exp (x (t+1) ) ; 
end 

if (S >Maxiraum_Limit ) 

disp (’ Eventhough input is bounded output is 
unbounded ’ ) 

disp (’The given system is unstable’); 
disp ( ’ S = ’ ) ; 
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S 

else 

disp( ’The given system is stable ’) ; 
disp (S) ; 
end 


Scilab code Exa 1.13 Determination of stablility of a given system 


Scilab code Exa 1.H3 //clear// 

//Example 1 . 1 3 : Determination of stablility of a 
given system 
//Page 49 

//given system y(t) = t.x(t) 

clear ; 
clc ; 

x = [1,2,3,4,0,2,1,3,5,81; //Assign some input 
Maximum_Lirait = 10; 


S = 0; 

for t = 0 : Maximum_Lirait -1 
S = S+t*x (t+1) ; 
end 

if (S >Maximum_Limit ) 

disp ( ’ Eventhough input is 
unbounded ’ ) 

disp (’The given system is 
disp ( ’ S = ’ ) ; 

S 

else 

disp(’The given system is 
disp (’The value of S ='); 
S 


bounded output 
unstable ’ ) ; 


stable ’ ) ; 


i s 


end 
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Determination of stablility of a given system 


//clear// 

//Example 1 . 1 3 : Determination of stablility of a 
given system 
//Page 49 

//given system y(t) = t.x(t) 

clear ; 
clc ; 

x = [1,2,3,4,0,2,1,3,5,81; //Assign some input 
Maximum_Lirait = 10; 


S = 0; 

for t = 0 : Maximum_Limit -1 
S = S+t*x (t+1) ; 
end 

if (S >Maximum_Limit ) 

disp ( ’ Eventhough input is 
unbounded ’ ) 

disp(’The given system is 
disp ( ’ S = ’ ) ; 

S 

else 

dispf’The given system is 
disp (’The value of S = ’); 
S 


bounded output 
unstable ’ ) ; 


stable ’ ) ; 


i s 


end 


Scilab code Exa 1.14 Time Invariance Property 

//clear// 

//Example 1.14: classification of a system: Time 
Invariance Property 
//Page 51 
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//To check whether the given system is a Time 
variant (or) Time In — variant 
// The given discrete signal is y(t) = sin(x(t)) 

clear ; 
clc ; 

to = 2; //Assume the amount of time shift =2 
T = 10; //Length of given signal 

for t = 1 : T 

x(t) = (2*°/ 0 pi/T)*t; 
y(t) = sin (x (t ) ) ; 
end 

//First shift the input signal only 

Input_shift = sin (x (T-to ) ) ; 

Output_shift = y(T-to); 

if ( Input_shif t == Output_shif t ) 

dispC'The given discrete system is a Time In- 
variant system’); 

else 

dispC’The given discrete system is a Time Variant 
system ’ ) ; 

end 

Scilab code Exa 1.14 

Time Invariance Property 

//clear// 

//Example 1.14: classification of a system : Time 
Invariance Property 
//Page 51 

//To check whether the given system is a Time 
variant (or) Time In — variant 
// The given discrete signal is y(t) = sin(x(t)) 

clear ; 
clc ; 

to = 2; //Assume the amount of time shift =2 
T = 10; //Length of given signal 

for t = 1 : T 
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x(t) = (2*°/ 0 pi/T) *t ; 
y(t) = sin (x (t ) ) ; 
end 

//First shift the input signal only 

Input_shift = sin (x (T-to ) ) ; 

Output_shift = y(T-to); 

if ( Input_shif t == Output_shif t ) 

dispC'The given discrete system is a Time In- 
variant system’); 

else 

disp(’The given discrete system is a Time Variant 
system ’ ) ; 

end 


Scilab code Exa 0.15 Sum of two complex exponentials 

//clear// 

//Example 1 . 5 : To express sum of two complex 
exponentials 
//as a single sinusoid 

clear ; 
clc ; 
close ; 

t =0 : 1/100 : 2* %pi ; 

xl = exp (sqrt (-1) *2*t ) ; 
x 2 = exp (sqrt (-1) *3*t) ; 

x = x 1 +x2 ; 

for i = l:length(x) 

X ( i ) = sqrt ( (real (x(i) ). “2) +( iraag (x(i) ). ~2) ) ; 
end 

plot (t , X) ; 

xtitle(’Full wave rectified sinusoid time t ’ , ’ 
Magnitude ’ ) ; 
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Scilab code Exa 0.15 

Sum of two complex exponentials 

//clear// 

//Example 1 . 5 : To express sum of two complex 
exponentials 
//as a single sinusoid 

clear ; 
clc ; 
close ; 

t =0 : 1/100 : 2* °/„pi ; 

xl = exp (sqrt (-1) *2*t ) ; 

x 2 = exp (sqrt (-1) *3*t) ; 

x = x 1 +x2 ; 

for i = l:length(x) 

X ( i ) = sqrt (( real (x ( i) ). ~2) +( imag (x (i) ). “2) ) ; 
end 

plot (t , X) ; 

xtitle(’Full wave rectified sinusoid’, ’time t ’ , ’ 
Magnitude ’ ) ; 


Scilab code Exa 1.15 Classification of a System 

//clear// 

//Example 1.15: Classification of a System : Time 
Invariance Property 
//Page 51 

//To check whether the given system is a Time 
variant (or) Time In — variant 
// The given discrete signal is y[n] = n.x[n] 

clear ; 
clc ; 


37 


8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 


no = 2 ; //Assume the amount of time shift =2 
L = 10; //Length of given signal 

for n = 1 : L 
x (n) = n; 
y (n) = n*x (n) ; 
end 

//First shift the input signal only 

Input_shift = x(L-no); 

Output_shift = y(L-no); 

if ( Input_shif t == Output_shif t ) 

disp(’The given discrete system is a Time In- 
variant system’); 

else 

dispC’The given discrete system is a Time Variant 
system ’ ) ; 

end 

Scilab code Exa 1.15 

Classification of a System 

//clear// 

//Example 1.15: Classification of a System : Time 
Invariance Property 
//Page 51 

//To check whether the given system is a Time 
variant (or) Time In — variant 
// The given discrete signal is y[n] = n.x[n] 

clear ; 
clc ; 

no = 2 ; //Assume the amount of time shift =2 
L = 10; //Length of given signal 

for n = 1 : L 
x (n) = n; 
y (n) = n*x (n) ; 
end 

//First shift the input signal only 

Input_shift = x(L-no); 
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Output_shift = y(L-no); 

if ( Input_shif t == Output_shif t ) 

dispC’The given discrete system is a Time In- 
variant system’); 

else 

dispC’The given discrete system is a Time Variant 
system ’ ) ; 

end 


Scilab code Exa 1.16 Time Invariance Property 

//clear// 

//Example 1 . 1 6 : C 1 as s i f i c at i o n of system : Time 
Invariance Property 
//Page 52 

//To check whether the given system is a Time 
variant (or) Time In — variant 
// The given discrete signal is y(t) = x(2t) 

clear ; 
clc ; 

to = 2; //Assume the amount of time shift =2 
T = 10; //Length of given signal 

x = [1,2,3,4,5,6,7,8,9,10]; 
y = zeros ( 1 , length (x) ) ; 
for t = 1 : length (x) /2 
y (t ) = x (2*t ) ; 

end 

//First shift the input signal only 

Input_shift = x(T-to); 

Output_shift = y(T-to); 

if ( Input_shif t == Output_shif t ) 

dispC’The given discrete system is a Time In- 
variant system’); 
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else 

dispC’The given discrete system is a Time Variant 
system ’ ) ; 

end 

Scilab code Exa 1.16 

Time Invariance Property 

//clear// 

//Example 1 . 1 6 : C 1 as s i f i c a t i o n of system : Time 
Invariance Property 
//Page 52 

//To check whether the given system is a Time 
variant (or) Time In — variant 
// The given discrete signal is y(t) = x(2t) 

clear ; 
clc ; 

to = 2; //Assume the amount of time shift =2 
T = 10; //Length of given signal 

x = [1,2,3,4,5,6,7,8,9,10]; 
y = zeros ( 1 , length (x) ) ; 
for t = 1 : length (x) /2 
y (t) = x (2*t ) ; 

end 

//First shift the input signal only 

Input_shift = x(T-to); 

0utput_shift = y(T-to); 

if ( Input_shif t == 0utput_shif t ) 

dispC’The given discrete system is a Time In- 
variant system’); 

else 

dispC’The given discrete system is a Time Variant 
system ’ ) ; 

end 
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Scilab code Exa 1.17 Linearity Property 

//clear// 

/ / Example 1.17: Classification of system: Linearity 
Property 
//Page 54 

//To check whether the given discrete system is a 
Linear System (or) Non-Linear System 
//Given discrete system y(t)= t*x(t) 

clear ; 
clc ; 

xl = [1 , 1 ,1 ,1] ; 

x2 = [2 , 2 , 2 , 2] ; 

a = 1; 
b = 1; 

for t = l:length(xl) 

x3(t) = a*xl (t ) +b*x2 (t ) ; 
end 

for t = l:length(xl) 
yl (t) = t*xl (t) ; 
y2 (t) = t * x2 (t) ; 
y3 (t) = t*x3 (t) ; 
end 

for t = l:length(yl) 

z(t) = a*yl (t)+b*y2 (t) ; 
end 

count = 0; 
for n =1 : length (yl ) 
if ( y3 (t) == z(t) ) 
count = count + 1 ; 
end 
end 
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if(count == length(y3)) 

disp(’Since It satisifies the superposition 
principle ’ ) 

dispC’The given system is a Linear system’) 

y3 

z 

else 

dispC ’Since It does not satisify the 
superposition principle ’) 
dispC’The given system is a Non-Linear system’) 

end 

Scilab code Exa 1.17 

Linearity Property 


//clear// 

//Example 1 . 1 7 : C 1 as s i f i c at i o n of system : Linearity 
Property 
//Page 54 

//To check whether the given discrete system is a 
Linear System (or) Non-Linear System 
//Given discrete system y ( t )= t *x ( t ) 

clear ; 
clc ; 

xl = [1 , 1 ,1 ,1] ; 

x2 = [2 , 2 , 2 , 2] ; 

a = 1; 
b = 1; 
f or t = 
x3 (t ) 
end 

f or t = 
yl (t) 
y2(t) 
y3 (t) 
end 

for t = l:length(yl) 


length (xl ) 
a*xl ( t ) +b*x2 ( t ) ; 

length (xl ) 
t*xl (t) ; 
t * x2 (t) ; 
t*x3 (t) ; 
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z(t) = a*yl (t ) +b*y2 (t ) ; 
end 

count = 0; 
for n =1 : length (yl ) 
if ( y3 (t) == z(t) ) 
count = count+1 ; 
end 
end 

if(count == length(y3)) 

dispC’Since It satisifies the superposition 
principle ’ ) 

dispC’The given system is a Linear system’) 

y3 

z 

else 

dispC’Since It does not satisify the 
superposition principle ’) 
dispC’The given system is a Non-Linear system’) 

end 


Scilab code Exa 1.18 Linearity Property 

//clear// 

/ / Example 1.18:Classsification of a system: Linearity 
Property 
//Page 54 

//To check whether the given discrete system is a 
Linear System (or) Non-Linear System 
//Given discrete system y(t)= (x(t)~2) 

clear ; 
clc ; 

xl = [1 , 1 ,1 ,1] ; 

x2 = [2 ,2,2, 2] ; 


43 


10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

1 

2 


a = 1; 

b = 1; 

for t = l:length(xl) 

x3(t) = a*xl (t) +b*x2 (t) ; 
end 

for t = l:length(xl) 
yl (t) = (xl (t) ~2) ; 
y2 (t) = (x2(t)~2); 
y3 (t) = (x3 (t) ~2) ; 
end 

for t = l:length(yl) 

z(t) = a*yl (t) +b* y 2 (t) ; 
end 

count = 0; 
for n =1 : length (yl ) 
if ( y3 (t) == z ( t ) ) 
count = count + 1 ; 
end 
end 

if(count == length(y3)) 

dispC’Since It satisifies the superposition 
principle ’ ) 

disp(’The given system is a Linear system’) 

y3 

z 

else 

dispC’Since It does not satisify the 
superposition principle ’) 
disp(’The given system is a Non-Linear system’) 

end 

Scilab code Exa 1.18 

Linearity Property 

//clear// 

//Example 1 . 1 8 : C 1 as s s i f i c a t i o n of a system : Linearity 
Property 
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//Page 54 

//To check whether the given discrete system is a 
Linear System (or) Non-Linear System 
//Given discrete system y(t)= (x(t)~2) 

clear ; 
clc ; 

xl = [1 , 1 ,1 ,1] ; 

x2 = [2 , 2 , 2 , 2] ; 

a = 1; 
b = 1; 

for t = l:length(xl) 

x3(t) = a*xl (t)+b*x2 (t) ; 
end 

for t = l:length(xl) 
yl (t) = (xl (t) ~2) ; 
y2 (t) = (x2(t)~2) ; 
y3 (t) = (x3 (t) ~2) ; 
end 

for t = l:length(yl) 

z(t) = a*yl (t ) +b*y2 (t ) ; 
end 

count = 0; 
for n =1 : length (yl ) 
if ( y3 (t) == z ( t ) ) 
count = count+1 ; 
end 
end 

if(count == length(y3)) 

disp(’Since It satisifies the superposition 
principle ’ ) 

disp(’The given system is a Linear system’) 

y3 

z 

else 

disp( ’Since It does not satisify the 
superposition principle ’) 
disp(’The given system is a Non-Linear system’) 

end 
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Scilab code Exa 1.20 Linearity Property 

//clear// 

//Example 1 . 2 0 : C 1 a s s s i f i c a t i o n of a system : Linearity 
Property 
//Page 55 

//To check whether the given discrete system is a 
Linear System (or) Non-Linear System 
//Given discrete system y[n])= 2*x[n]+3 

clear ; 
clc ; 

xl = [1 , 1 ,1 , 1] ; 

x2 = [2 , 2 , 2 , 2] ; 

a = 1; 
b = 1; 

for n = l:length(xl) 

x3(n) = a*xl (n) +b*x2 (n) ; 
end 

for n = l:length(xl) 
yl (n) = 2*xl (n) +3 ; 
y2(n) = 2*x2 (n) +3 ; 
y3 (n) = 2*x3 (n) +3 ; 
end 

for n = l:length(yl) 

z (n) = a* y 1 (n)+b*y2(n) ; 
end 

count = 0; 
for n =1 : length (yl ) 
if (y3 (n) == z (n) ) 
count = count + 1 ; 
end 
end 
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if(count == length(y3)) 

dispC’Since It satisifies the superposition 
principle ’ ) 

dispC’The given system is a Linear system’) 

y3 

z 

else 

dispC ’Since It does not satisify the 
superposition principle ’) 
dispC’The given system is a Non-Linear system’) 

end 

Scilab code Exa 1.20 

Linearity Property 

//clear// 

//Example 1 . 2 0 : C 1 as s s i f i c a t i o n of a system : Linearity 
Property 
//Page 55 

//To check whether the given discrete system is a 
Linear System (or) Non-Linear System 
//Given discrete system y[n])= 2*x[n]+3 

clear ; 
clc ; 

xl = [1 , 1 ,1 ,1] ; 

x2 = [2 , 2 , 2 , 2] ; 

a = 1; 
b = 1; 

for n = l:length(xl) 

x3(n) = a*xl (n) +b*x2 (n) ; 
end 

for n = l:length(xl) 
yl(n) = 2*xl (n) +3 ; 
y2(n) = 2*x2 (n) +3 ; 
y3 (n) = 2*x3 (n) +3 ; 
end 

for n = l:length(yl) 
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z(n) = a*yl (n) +b*y2 (n) ; 
end 

count = 0; 
for n =1 : length (yl ) 
if (y3 (n) == z (n) ) 
count = count+1 ; 
end 
end 

if(count == length(y3)) 

disp(’Since It satisifies the superposition 
principle ’ ) 

disp(’The given system is a Linear system’) 

ys 

z 

else 

dispC’Since It does not satisify the 
superposition principle ’) 
disp(’The given system is a Non-Linear system’) 

end 
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Chapter 4 

Linear Time Invariant Systems 


Scilab code Exa 2.1 Linear Convolution Sum 

//clear// 

//Example 2.1: Linear Convolution Sum 
//page 80 

clear ; 
close ; 
clc ; 

h = [0 ,0 , 1 , 1 , 1 ,0 ,0] ; 

N 1 = -2:4; 

x = [0,0, 0.5, 2, 0,0,0]; 

N2 = -2:4; 
y = convol (x , h) ; 
for i = l:length(y) 
if (y (i) <=0 . 0001) 
y ( i) =0 ; 
end 
end 

N = -4:8; 
f igure 
a=gca () ; 
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plot2d3 ( ’ gnn 7 ,N1 ,h) 

xt it le ( 'Impulse Response 7 , 7 n 7 , 7 h [ n ] 7 ) ; 
a . thickness = 2 ; 
f igure 
a=gca () ; 

plot2d3 ( 7 gnn 7 , N2 , x) 

xt it le ( ’Input Response 7 , 7 n 7 , 7 x [ n ] 7 ) ; 
a . thickness = 2 ; 
f igure 
a=gca () ; 

plot2d3 ( 7 gnn 7 , N , y ) 

xt it le ( ’Output Response 7 , 7 n 7 , 7 y[n] 7 ) ; 
a . thickness = 2 ; 

Scilab code Exa 2.1 

Linear Convolution Sum 

//clear// 

//Example 2.1: Linear Convolution Sum 
//page 80 

clear ; 
close ; 
clc ; 

h = [0 ,0 , 1 , 1 , 1 ,0 ,0] ; 

N 1 = -2:4; 

x = [0,0, 0.5, 2, 0,0,0]; 

N2 = -2:4; 
y = convol (x , h) ; 
for i = l:length(y) 
if (y (i) <=0 . 0001) 
y ( i ) =0 ; 
end 
end 

N = -4:8; 
f igure 
a=gca () ; 

plot2d3( 7 gnn 7 ,N1 ,h) 
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xt it le ( ’Impulse Response ’ , ’ n ’ , ’ h [ n ] ’) ; 
a . thickness = 2 ; 
f igure 
a=gca () ; 

plot2d3 ( ’ gnn ’ , N2 , x) 

xtitle( "Input Response ’ , ’ n ’ , ’ x [ n ] ’) ; 
a . thickness = 2 ; 
f igure 
a=gca () ; 

plot2d3 ( ’ gnn ’ , N , y ) 

xt it le ( ’Output Response ’ , ’ n ’ , ’ y [ n ] ’) ; 
a . thickness = 2 ; 


Scilab code Exa 2.3 Scilab code Exa 2.3 Convolution of x[n] and Unit 
Impulse response h[n] 

//clear// 

//Example 2 . 3 : Convolution Sum : Convolution of x[n] 
and 

//Unit Impulse response h[n] 

clear ; 
close ; 
clc ; 

Max_Limit = 10; 
h = ones ( 1 , Max_Limit ) ; 

N1 = 0 : Max_Limit -1 ; 

Alpha = 0.5; //alpha < 1 
for n = l:Max_Lirait 

x (n) = (Alpha ~(n-l))*l; 
end 

N2 = 0 : Max_Limit -1 ; 
y = convol (x , h) ; 

N = 0 : 2* Max_Limit -2 ; 
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f igure 
a=gca () ; 

plot2d3 ( ’ gnn ’ ,N1 ,h) 

xtitle( ’Impulse Response Fig 2.5.(b)’,’n’,’h[n]’); 
a . thickness = 2 ; 
f igure 
a=gca () ; 

plot2d3( ’ gnn ’ , N2 , x) 

xt it le ( ’Input Response Fig 2.5. (a) ’ n ’ , ’ x [ n ] ’); 

a . thickness = 2 ; 
f igure 
a=gca () ; 

plot2d3 ( ’ gnn ’ , N ( 1 : Max _ Limit ) ,y(l:Max_Limit) ,5) 
xtitle (’ Output Response Fig 2 . 7 ’ , ’ n ’ , ’ y [ n ] ’ ) ; 
a . thickness = 2 ; 

Convolution of x[nj and Unit Impulse response h[nj 

//clear// 

//Example 2 . 3 : Convolution Sum : Convolution of x[n] 
and 

//Unit Impulse response h[n] 

clear ; 
close ; 
clc ; 

Max_Limit = 10; 
h = ones ( 1 , Max_Limit ) ; 

N1 = 0 : Max_Limit -1 ; 

Alpha = 0.5; //alpha < 1 
for n = l:Max_Lirait 

x(n)= ( Alpha ~ (n - 1 ))* 1 ; 
end 

N2 = 0 : Max_Limit -1 ; 
y = convol (x , h) ; 

N = 0 : 2*Max_Limit -2 ; 
f igure 
a=gca () ; 

plot2d3( ’gnn ’ , N1 ,h) 
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xtitle( ’Impulse Response Fig 2 . 5 . ( b ) ’ , ’ n ’ , ’ h [ n ] ’ ) ; 
a . thickness = 2 ; 
f igure 
a=gca () ; 

plot2d3 ( ’ gnn ’ , N2 , x) 

xtitle( "Input Response Fig 2 . 5 . ( a ) ’ , ’ n ’ , ’ x [ n ] ’ ) ; 
a . thickness = 2 ; 
f igure 
a=gca () ; 

plot2d3 ( ’gnn ’ , N (1 : Max_Lirait) , y ( 1 : Max_Limit) ,5) 
xtitle (’ Output Response Fig 2 . 7 ’ , ’ n ’ , ’ y [ n ] ’ ) ; 
a . thickness = 2 ; 


Scilab code Exa 2.4 Convolution Sum of finite duration sequences 

//clear// 

//Example 2 . 4 : Convolution Sum of finite duration 
sequences 

clear ; 
close ; 
clc ; 

x = ones (1,5); 

N1 =0 : length (x) -1 ; 

Alpha = 1.4; //alpha > 1 
for n = 1:7 

h(n)= ( Alpha ~ (n - 1 ))* 1 ; 
end 

N2 =0 : length (h) -1 ; 
y = convol (x , h) ; 

N = 0 : length (x) + length (h) -2 ; 
f igure 
a=gca () ; 

plot2d3( ’gnn ’ ,N2 ,h) 
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xt it le ( ’Impulse Response ’ , ’ n ’ , ’ h [ n ] ’) ; 
a . thickness = 2 ; 
f igure 
a=gca () ; 

plot2d3 ( ’ gnn ’ , N1 , x) 

xtitle( "Input Response ’ , ’ n ’ , ’ x [ n ] ’) ; 
a . thickness = 2 ; 
f igure 
a=gca () ; 

plot2d3 ( ’ gnn ’ , N , y ) 

xt it le ( ’Output Response ’ , ’ n ’ , ’ y [ n ] ’) ; 
a . thickness = 2 ; 

Scilab code Exa 2.4 

Convolution Sum of finite duration sequences 

//clear// 

//Example 2 . 4 : Convolution Sum of finite duration 
sequences 

clear ; 
close ; 
clc ; 

x = ones (1,5); 

N1 =0 : length (x) -1 ; 

Alpha = 1.4; //alpha > 1 
for n = 1:7 

h(n)= ( Alpha ~ (n - 1 ))* 1 ; 
end 

N2 =0 : length (h) -1 ; 
y = convol (x , h) ; 

N = 0 : length (x) + length (h) -2 ; 
f igure 
a=gca () ; 

plot2d3( ’gnn ’ ,N2 ,h) 

xtitle( "Impulse Response ’ , ’ n ’ , ’ h [ n ] ’) ; 
a . thickness = 2 ; 
f igure 
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a=gca () ; 

plot2d3 ( ’ gnn ’ ,N1 ,x) 

x t i 1 1 e ( 'Input Response ’ , ’ n ’ , ’ x [ n ] ’) ; 
a . thickness = 2 ; 
f igure 
a=gca () ; 

plot2d3 ( ’ gnn ’ , N , y ) 

xt it le ( ’Output Response ’ , ’ n ’ , ’ y [ n ] ’) ; 
a . thickness = 2 ; 


Scilab code Exa 2.5 Convolution Sum of input sequence 

//clear// 

//Example 2 . 5 : Convolution Sum of input sequence x[n 
] = (2 A n) • u [ — n ] 

/ / and h [ n ] = u [ n ] 

clear ; 
close ; 
clc ; 

Max_Limit = 10; 
h = ones ( 1 , Max_Limit ) ; 

N2 =0 : length (h) -1 ; 
for n = l:Max_Limit 

xl (n) = (2~ ( - (n-1) ) ) *1 ; 
end 

x = xl ($ : -1 : 1) ; 

N1 = -length (x) +1 : 0 ; 
y = convol (x , h) ; 

N = -length (x) +1 : length (h) -1 ; 
f igure 
a=gca () ; 

plot2d3 ( ’gnn ’ ,N2 ,h) 

xt it le ( "Impulse Response ’ , ’ n ’ , ’ h [ n ] ’) ; 
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2 ; 


a . thickness = 
f igure 
a=gca () ; 
a.y_location = ’’origin”; 
plot2d3 ( ’ gnn ’ ,N1 ,x) 

xt it le ( ’Input Response Fig 2.11(a) ’, ’ n ’ , ’ x [ n ] ’); 
a . thickness = 2 ; 
f igure 
a=gca () ; 

a.y_location =” origin”; 
plot2d3 ( ’ gnn ’ , N , y ) 

xt it le (’ Output Response Fig 2.11(b) ’ , ’ n ’ , ’ y [ n ] ’ ) ; 
a . thickness = 2 ; 

Scilab code Exa 2.5 

Convolution Sum of input sequence 

//clear// 

//Example 2 . 5 : Convolution Sum of input sequence x[n 
] = ( 2 ~ n ) ,u[-n] 

/ / and h [ n ] = u [ n ] 

clear ; 
close ; 
clc ; 

Max_Limit = 10; 
h = ones ( 1 , Max_Limit ) ; 

N2 =0 : length (h) -1 ; 
for n = l:Max_Limit 

xl (n) = (2~(-(n-l)))*l; 
end 

x = xl ($ : -1 : 1) ; 

N1 = -length (x) +1 : 0 ; 
y = convol (x , h) ; 

N = -length (x) +1 : length (h) -1 ; 
f igure 
a=gca () ; 

plot2d3( ’gnn ’ , N2 ,h) 
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xt it le ( ’Impulse Response ’ , ’ n ’ , ’ h [ n ] ’) ; 
a . thickness = 2 ; 
f igure 
a=gca () ; 

a.y_location = ’’origin”; 
plot2d3 ( ’ gnn ’ ,N1 ,x) 

xtitle (’ Input Response Fig 2.11(a) ’, ’n’, ’x[n] ’); 
a . thickness = 2 ; 
f igure 
a=gca () ; 

a.y_location =” origin”; 
plot2d3( ’gnn ’ , N , y ) 

xt it le (’ Output Response Fig 2.11(b) ’ , ’ n ’ , ’ y [ n ] ’ ) ; 
a . thickness = 2 ; 


Scilab code Exa 2.6 Convolution Integral of input 

//clear// 

//Example 2 . 6 : Convolution Integral of input x(t) = ( 
e " — at ) . u ( t ) 

/ / and h ( t ) =u ( t ) 

clear ; 
close ; 
clc ; 

Max_Limit = 10; 
h = ones ( 1 , Max_Limit ) ; 

N2 =0 : length (h) -1 ; 
a = 0.5; //constant a>0 
for t = l:Max_Lirait 
x (t ) = exp (-a* (t-1) ) ; 
end 

N1 =0 : length (x) -1 ; 
y = convol (x , h) -1 ; 
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N = 0 : length (x) + length (h) -2 ; 
f igure 
a=gca () ; 
plot2d (N2 , h) 

xt it le ( 'Impulse Response ’ , ’ t ’ , ’ h ( t ) ’ ) ; 
a . thickness = 2 ; 
f igure 
a=gca () ; 
plot2d(Nl , x) 

xt it le ( "Input Response ’,’t’,’x(t)’); 
a . thickness = 2 ; 
f igure 
a=gca () ; 

plot2d(N(l: Max_Limit) ,y(l: Max_Limit)) 
xt it le ( ’Output Response ’ , ’ t ’ , ’ y ( t ) ’) ; 
a . thickness = 2 ; 

Scilab code Exa 2.6 

Convolution Integral of input 

//clear// 

//Example 2 . 6 : Convolution Integral of input x(t) = ( 
e ~ — at ) • u ( t ) 

/ / and h ( t ) =u ( t ) 

clear ; 
close ; 
clc ; 

Max_Limit = 10; 
h = ones ( 1 , Max_Limit ) ; 

N2 =0 : length (h) -1 ; 
a = 0.5; //constant a>0 
for t = l:Max_Lirait 
x (t ) = exp (-a* (t-1) ) ; 
end 

N1 =0 : length (x) -1 ; 
y = convol (x , h) -1 ; 

N = 0 : length (x) + length (h) -2 ; 
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f igure 
a=gca () ; 
plot2d (N2 , h) 

xt it le ( ‘Impulse Response 7 , 7 t 7 , 7 h ( t ) 7 ) ; 
a . thickness = 2 ; 
f igure 
a=gca () ; 
plot2d(Nl ,x) 

xt it le ( ’Input Response 7 , 7 t 7 ,’x(t) 7 ); 
a . thickness = 2 ; 
f igure 
a=gca () ; 

plot2d(N(l: Max_Limit) , y ( 1 : Max_Limit)) 
xt it le ( ’Output Response 7 , 7 t 7 , 7 y ( t ) 7 ) ; 
a . thickness = 2 ; 


Scilab code Exa 2.7 Convolution Integral of Untie duration signals 


Scilab code Exa 217 //clear// 

//Example 2 . 7 : Convolution Integral of fintie 
duration signals 
/ / page99 

clear ; 
close ; 
clc ; 

T = 10; 

x = ones(l,T); //Input Response 
for t = 1 : 2* T 

h(t) = t — 1 ; //Impulse Response 

end 

N1 = 0 : length (x) -1 ; 

N2 = 0 : length (h) -1 ; 
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y = convol (x , h) ; 

N = 0 : length (x) + length (h) -2 ; 
f igure 
a=gca () ; 

a. x_location=” origin” ; 
plot2d (N2 , h) 

xtitle( ’Impulse Response ’ , ’ t ’ , ’ h ( t ) ’) ; 
a . thickness = 2 ; 
f igure 
a=gca () ; 
plot2d(Nl ,x) 

xtitle( ‘Input Response ’,’t’,’x(t)’); 
a . thickness = 2 ; 
f igure 
a=gca () ; 
plot 2d (N , y ) 

xt it le ( ’Output Response ’ , ’t ’ , ’y(t) ’) ; 
a . thickness = 2 ; 


Convolution Integral of fintie duration signals 

//clear// 

//Example 2 . 7 : Convolution Integral of fintie 
duration signals 
/ / page99 

clear ; 
close ; 
clc ; 

T = 10; 

x = ones(l,T); //Input Response 
for t = 1 : 2* T 

h(t) = t — 1 ; //Impulse Response 
end 

N1 = 0 : length (x) -1 ; 

N2 = 0 : length (h) -1 ; 
y = convol (x , h) ; 

N = 0 : length (x) + length (h) -2 ; 
f igure 
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a=gca () ; 

a. x_location=” origin” ; 
plot2d (N2 , h) 

xt it le ( ‘Impulse Response 7 , 7 t 7 , 7 h ( t ) 7 ) ; 
a . thickness = 2 ; 
f igure 
a=gca () ; 
plot2d(Nl ,x) 

xtitle ( ’Input Response 7 , 7 t 7 ,’x(t) 7 ); 
a . thickness = 2 ; 
f igure 
a=gca () ; 
plot2d (N , y) 

xt it le ( ’Output Response 7 , 7 t 7 , 7 y ( t ) 7 ) ; 
a . thickness = 2 ; 


Scilab code Exa 2.8 Convolution Integral of input 

//clear// 

//Example 2 . 8 : Convolution Integral of input x(t)=(e 
"2 t ) . u( — t ) and 

//MO = u(t-3) 

clear ; 
close ; 
clc ; 

Max_Limit = 10; 

h =[0,0,0, ones(l, M ax _Limit-3)]; //h(n— 3) 
a = 2 ; 
t = -9:0; 
x= exp ( a*t ) ; 

//x = xl ( $ : — 1:1) 

N2 = 0 : length (h) -1 ; 

N1 = -length (x) +1 : 0 ; 
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tl = -6:3; 

yl = (1/a) *exp (a* (tl -3) ) ; 
y2 = ( 1/a) * ones ( 1 , Max_Lirait ) ; 
y = [yl y2] 

N = -length (h) +1 : length (x) -1 ; 
f igure 
a=gca () ; 

a. x_location = ” origin” ; 

a. y_location = ” origin” ; 

plot2d(-Max_Limit+l : 0 ,h($ : -1 : 1) ) 

xt it le ( ’Impulse Response ’ , ’ t ’ , ’ h ( t — T ) ’ ) ; 

a . thickness = 2 ; 

f igure 

a=gca () ; 

a.y_location = ’’origin”; 
plot2d (t , x) 

xt it le ( 'Input Response ’t ’x(t) ’); 
a . thickness = 2 ; 
f igure 
a=gca () ; 

a.y_location = ’’origin”; 
a. x_location = ’’origin”; 
a. data_bounds = [-10 ,0; 13 , 1] ; 
plot2d(-Max_Limit+4:Max_Limit+3,y) 
xt it le ( ’Output Response ’ , ’ t ’ , ’y(t) ’) ; 
a . thickness = 2 ; 

Scilab code Exa 2.8 

Convolution Integral of input 

//clear// 

//Example 2 . 8 : Convolution Integral of input x(t)=(e 
"2 t ) . u( — t ) and 
// h ( t ) = u ( t —3) 

clear ; 
close ; 
clc ; 
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Max_Limit = 10; 

h =[0,0,0, ones(l, M ax _Limit-3)]; //h(n — 3) 
a = 2 ; 
t = -9:0; 
x= exp ( a*t ) ; 

//x = xl ( $ : — 1:1) 

N2 = 0 : length (h) -1 ; 

N1 = -length (x) +1 : 0 ; 
tl = -6:3; 

yl = (1/a) *exp (a* (tl -3) ) ; 
y2 = (1/a) *ones (1 , Max_Lirait ) ; 
y = [yl y2] 

N = -length (h) +1 : length (x) -1 ; 
f igure 
a=gca () ; 

a.x_location=” origin” ; 

a.y_location=” origin” ; 

plot2d(-Max_Limit+l : 0 , h ( $ : -1 : 1) ) 

xt it le ( ’Impulse Response ’ , ’ t ’ , ’ h ( t — T) ’ ) ; 

a . thickness = 2 ; 

f igure 

a=gca () ; 

a.y_location = ’’origin”; 
plot2d (t , x) 

xt it le ( 'Input Response ’, ’t ’, ’x(t) ’); 
a . thickness = 2 ; 
f igure 
a=gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a. data_bounds = [-10 ,0; 13 , 1] ; 
plot2d(-Max_Limit+4: Max_Limit+3 , y ) 
xt it le ( ’Output Response ’ , ’ t ’ , ’y(t) ’) ; 
a . thickness = 2 ; 
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Chapter 6 

Fourier Series Repreentation of 
Periodic Signals 


Scilab code Exa 3.2 CTFS of a periodic signal x(t) 

//clear// 

//Example 3 . 2 : CTFS of a periodic signal x(t) 
//Expression of continuous time signal 
//using continuous time fourier series 

clear ; 
close ; 
clc ; 

t = -3:0.01:3; 

// 1 1 = — %pi *4 : ( %pi *4) / 100 : %pi *4 ; 

// 1 2 =— %pi * 6 : ( %pi *6) / 1 00 : %pi *6 ; 

xot = ones ( 1 , length (t )) ; 
xlt = (1/2) * c o s (°/ 0 pi*2*t) ; 
xot_xlt = xot+xlt; 
x2t = cos ( °/ 0 pi *4*t ) ; 
xot_xlt_x2t = xot+xlt+x2t; 
x3t = (2/3)*cos(°/ 0 pi*6*t); 
xt = xot +xlt +x2t +x3t ; 
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// 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a . data_bounds = [-4 , 0 ; 2 4]; 
plot (t , xot ) 
ylabel ( ’ t ’ ) 
title ( ’ xot =1 ’ ) 

// 

f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a . data_bounds = [-4 , -3 ; 2 4]; 
plot (t , xlt ) 
ylabel ( ’ t ’ ) 

title(’xl(t) =1/2* cos (2 * pi * t ) ’ ) 

subplot (2,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a . data_bounds = [-4 , 0 ; 2 4]; 
plot (t , xot_xlt ) 
ylabel ( ’ t ’ ) 
t i 1 1 e ( ’xo(t) +x 1 ( t ) ’ ) 

// 

f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a . data_bounds = [-4 , -2 ; 4 2]; 
plot (t , x2t ) 
ylabel ( ’ t ’ ) 

title(’x2(t) =cos (4* pi * t ) ’ ) 

subplot (2,1,2) 
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a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a . data_bounds = [-4 , 0 ; 4 4]; 
plot (t ,xot_xlt_x2t) 
ylabel ( ’ t ’ ) 

title( ’xo(t) +x 1 ( t ) +x2 ( t ) ’ ) 

// 

f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 
a.x_location = ’’origin”; 
a . data_bounds = [-4 , -3 ; 4 3]; 
plot (t , x3t ) 
ylabel ( ’ t ’ ) 

title(’xl(t) =2/3*cos (6* pi*t ) ’ ) 
subplot (2,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a . dat a_bounds = [ -4 , -3 ; 4 3]; 
plot (t , xt ) 
ylabel ( ’ t ’ ) 

title ( ’x(t )=xo ( t )+xl ( t )+x2 ( t )+x3 ( t ) ’ ) 

Scilab code Exa 3.2 

CTFS of a periodic signal x(t) 

//clear// 

//Example 3 . 2 : CTFS of a periodic signal x(t) 
//Expression of continuous time signal 
//using continuous time fourier series 

clear ; 
close ; 
clc ; 

t = -3:0.01:3; 
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// 1 1 = — %pi *4 : ( %pi *4) / 100 : %pi *4 ; 
// 1 2 =— %pi * 6 : ( %pi * 6 ) / 1 0 0 : %pi * 6 ; 

xot = ones ( 1 , length (t )) ; 
xlt = (l/2)*cos(°/ 0 pi*2*t); 
xot_xlt = xot+xlt; 
x2t = cos ( %pi *4*t ) ; 
xot_xlt_x2t = xot+xlt+x2t; 
x3t = (2/3)*cos(°/ 0 pi*6*t); 
xt = xot +xlt +x2t +x3t ; 

// 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a . data_bounds = [-4 , 0 ; 2 4]; 
plot (t , xot ) 
ylabel ( ’ t ’ ) 
title ( ’ xot =1 ’ ) 

// 

f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a . data_bounds = [-4 , -3 ; 2 4]; 
plot (t ,xlt) 
ylabel ( ’ t ’ ) 

title(’xl(t) =1/2* cos (2 * pi * t ) ’ ) 

subplot (2,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a . data_bounds = [-4 , 0 ; 2 4]; 
plot (t , xot_xlt ) 
ylabel ( ’ t ’ ) 
t i 1 1 e ( ’xo(t) +x 1 ( t ) ’ ) 

11 

f igure 
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subplot (2,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a . data_bounds = [-4 , -2 ; 4 2]; 
plot (t , x2t ) 
ylabel ( ’ t ’ ) 

title(’x2(t) =cos (4* pi * t ) ’ ) 

subplot (2,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a . data_bounds = [-4 , 0 ; 4 4]; 
plot (t ,xot_xlt_x2t) 
ylabel ( ’ t ’ ) 

title ( ’xo ( t )+xl ( t )+x2 ( t ) ’ ) 

// 

f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 
a.x_location = ’’origin”; 
a . data_bounds = [-4 , -3 ; 4 3]; 
plot (t , x3t ) 
ylabel ( ’ t ’ ) 

title(’xl(t) =2/3* cos (6 * pi * t ) ’ ) 

subplot (2,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a . dat a_bounds = [ -4 , -3 ; 4 3]; 
plot (t , xt ) 
ylabel ( ’ t ’ ) 

title ( ’x(t )=xo ( t )+xl ( t )+x2 ( t )+x3 ( t ) ’ ) 
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Scilab code Exa 3.3 Continuous Time Fourier Series Coefficients 


//clear// 

// Example3 . 3 : Continuous Time Fourier Series 
Coefficients of 

//a periodic signal x ( t ) = sin( Wot ) 

clear ; 
close ; 
clc ; 

t = 0:0.01:1; 

T = 1; 

Wo = 2* °/„pi / T ; 
xt = sin (Wo*t ) ; 
for k =0:5 

C(k+1,:) = exp ( - sqrt ( - 1 ) * Wo *t . *k) ; 
a(k+l) = xt*C(k+l,:)’/length(t); 
if(abs(a(k + l)) < = 0 . 0 1 ) 
a (k + 1 ) =0 ; 
end 
end 

a =a ’ ; 

ak = [-a , a (2 : $ ) ] ; 

Scilab code Exa 3.3 

Continuous Time Fourier Series Coefficients 

//clear// 

// Example3 . 3 : Continuous Time Fourier Series 
Coefficients of 

//a periodic signal x ( t ) = sin( Wot ) 

clear ; 
close ; 
clc ; 

t = 0:0.01:1; 

T = 1; 
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Wo = 2*%pi/T; 
xt = sin (Wo*t ) ; 
for k =0:5 

C (k + 1 , : ) = exp(-sqrt(-l)*Wo*t.*k) ; 

a(k + l) = xt *C (k + 1 ,:) ’ / length (t ) ; 
if(abs(a(k + l)) < = 0 . 0 1 ) 
a (k + 1 ) =0 ; 
end 
end 

a =a ’ ; 

ak = [-a , a (2 : $ ) ] ; 


Scilab code Exa 3.4 CTFS coefficients of a periodic signal 

//clear// 

// Example3 . 4 : CTFS coefficients of a periodic signal 
//x(t) = 1+sin (Wot) +2cos (Wot) + cos (2Wot+%pi/4) 

clear ; 
close ; 
clc ; 

t = 0:0.01:1; 

T = 1; 

Wo = 2* °/„pi / T ; 

xt =ones (1 , length(t))+sin(Wo*t)+2*cos(Wo*t)+cos (2*Wo 
*t+%pi/4) ; 
for k =0:5 

C(k+1,:) = exp ( - sqrt ( - 1 ) * Wo *t ■ *k) ; 
a(k+l) = xt*C(k+l,:)’/length(t); 
if (abs(a(k+l) ) <=0. 1) 
a (k + 1 ) =0 ; 
end 
end 

a =a ’ ; 
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a_conj =conj (a) ; 

ak = [a_conj($:-l:l),a(2:$)]; 

Mag_ak = abs(ak); 
for i = l:length(a) 

Phase_ak(i) = atan ( iraag ( ak ( i ) ) / ( real ( ak ( i ) ) 

+ 0 . 0001 ) ) ; 

end 

Phase_ak = Phase_ak ’ 

Phase_ak = [Phase_ak ( 1 : $ ) -Phase_ak ( $ - 1 : - 1 : 1 ) ] ; 
f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 
a.x_location = ’’origin”; 
plot2d3( ’ gnn ’ , [-k : k] ,Mag_ak ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title ( ’ abs ( ak) ’ ) 
xlabel ( ’ 

k’) 

subplot (2,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
plot2d3( ’gnn’ , [-k : k] ,Phase_ak ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title ( ’ <(ak ) ’ ) 
xlabel ( ’ 

k’) 

Scilab code Exa 3.4 

CTFS coefficients of a periodic signal 

//clear// 
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// Example3 . 4 : CTFS coefficients of a periodic signal 
//x(t) = 1+sin (Wot) +2cos (Wot) + cos (2Wot+%pi/4) 

clear ; 
close ; 
clc ; 

t = 0:0.01:1; 

T = 1; 

Wo = 2* °/„pi / T ; 

xt =ones (1 , length(t))+sin(Wo*t)+2*cos(Wo*t)+cos (2*Wo 
*t + °/„pi /4) ; 
for k =0:5 

C (k + 1 , : ) = exp(-sqrt(-l)*Wo*t.*k) ; 

a(k+l) = xt *C (k+1 ,:) J / length (t ) ; 
if (abs(a(k+l) ) <=0 . 1) 
a (k+1 ) =0 ; 
end 
end 

a =a ’ ; 

a_conj =conj (a) ; 

ak = [a_conj($:-l:l),a(2:$)]; 

Mag_ak = abs(ak); 
for i = l:length(a) 

Phase_ak(i) = atan ( imag ( ak ( i ) ) / ( real ( ak ( i ) ) 

+ 0 . 0001 ) ) ; 

end 

Phase_ak = Phase_ak ’ 

Phase_ak = [Phase_ak ( 1 : $ ) -Phase_ak ( $ - 1 : - 1 : 1 ) ] ; 
f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location = ” origin”; 
a. x.location = ’’origin”; 
plot2d3( ’ gnn ’ , [-k : k] ,Mag_ak ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title ( ’ abs ( ak) ’ ) 
xlabel ( ’ 
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k’) 

subplot (2,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
plot2d3( ’ gnn ’ , [-k : k] ,Phase_ak ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 
polyl . thickness = 3; 
title ( ’ <(ak ) ’ ) 
xlabel ( ’ 

k’) 


Scilab code Exa 3.5 CTFS coefficients of a periodic signal 

//clear// 

// Example3 . 5 : CTFS coefficients of a periodic signal 
//x(t) = 1, | 1 1 <T1 , and 0, T1 < 1 1 1 <T / 2 

clear ; 
close ; 
clc ; 

T =4; 

T 1 = T/4 ; 

t = -T1 : Tl/100 : T1 ; 

Wo = 2* °/„pi / T ; 
xt =ones (1 , length (t) ) ; 

// 

for k =0:5 

C(k+1,:) = exp ( - sqrt ( - 1 ) * Wo *t . *k) ; 
a (k + 1 ) = xt*C(k + l,:)’/length(t); 
if (abs(a(k+l) ) <=0 . 1) 
a (k + 1 ) =0 ; 
end 
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end 

a =a ’ ; 

a_conj = real ( a ( : ) ) -sqrt ( -1) * imag (a ( : ) ) ; 
ak = [a_conj($:-l:l) , ,a(2:$)]; 

k = 0:5; 

k = [-k ( $ : -1:1) ,k(2: $)] ; 

Spectrum_ak = ( 1/2) *real ( ak) ; 

// 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a. data_bounds = [-2 , 0 ; 2 , 2 ] ; 
plot2d (t , xt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title ( ’x( t ) ’) 
xlabel ( ’ 

t ’) 

// 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
plot2d3( ’ gnn ’ ,k, Spectrum _ak ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title ( ’ abs ( ak) ’ ) 
xlabel ( ’ 

k’) 

Scilab code Exa 3.5 

CTFS coefficients of a periodic signal 

//clear// 
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// Example3 . 5 : CTFS coefficients of a periodic signal 
//x(t) = 1, | 1 1 <T1 , and 0, T1 < 1 1 1 <T / 2 

clear ; 
close ; 
clc ; 

T =4; 

T 1 = T/4 ; 

t = -T 1 : Tl/100 : T1 ; 

Wo = 2* °/„pi / T ; 
xt =ones (1 , length (t) ) ; 

// 

for k =0:5 

C (k+1 , : ) = exp(-sqrt(-l)*Wo*t.*k) ; 
a(k+l) = xt *C (k+1 length (t ) ; 
if (abs(a(k+l) ) <=0. 1) 
a (k+1 ) =0 ; 
end 
end 

a =a ’ ; 

a_conj = real ( a ( : ) ) -sqrt ( -1) * imag (a ( : ) ) ; 
ak = [a_conj($:-l:l) , ,a(2:$)]; 

k = 0:5; 

k = [-k ( $ : -1: 1) ,k(2: $)] ; 

Spectrum_ak = ( 1 /2) * real ( ak) ; 

// 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a.data_bounds=[-2,0;2,2] ; 
plot2d (t , xt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 
polyl . thickness = 3; 
t itle ( ’x ( t ) ’) 
xlabel ( ’ 

t ’) 

// 


77 


38 

39 

40 

41 

42 

43 

44 

45 

46 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 


f igure 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
plot2d3( ’ gnn ’ ,k, Spectrum _ak ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title ( ’ abs ( ak) ’ ) 
xlabel ( ’ 

k’) 


Scilab code Exa 3.6 Time Shift Property of CTFS 

//clear// 

// Example3 . 6 : Time Shift Property of CTFS 

clear ; 
close ; 
clc ; 

T =4; 

T 1 = T/2 ; 
t = 0:T1/100:T1; 

Wo = 2* °/„pi / T ; 

gt =( 1/2) * ones ( 1 , length (t )) ; 
a(l)=0; / / k = 0 , ak =0 

d (1) =0; 
for k =1:5 

a (k + 1 ) = (sin (°/„pi*k/2) /(k*°/„pi) ) ; 
if(abs(a(k+l)) <=0 . 0 1 ) 
a (k+ 1 ) =0 ; 
end 

d(k + l) = a (k + 1 ) * exp ( - sqrt ( -1 ) *k* °/ 0 pi /2) ; 

end 
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k = 0:5 

disp (’ Fourier Series Coefficients of Square Wave’) 

a 

d i s p ( ’ Fourier Series Coefficients of g(t) =x ( t — 1 ) — 0 . 5 
’) 
d 

// 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a. data_bounds = [-l , -2; 1 ,4] ; 
plot2d([-t($:-l:l) , t ( 1 : $) ] , [-gt , gt ] ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title( ’g(t) ’) 
xlabel ( ’ 

t ’) 

Scilab code Exa 3.6 

Time Shift Property of CTFS 

//clear// 

// Example3 . 6 : Time Shift Property of CTFS 

clear ; 
close ; 
clc ; 

T =4; 

T 1 = T/2 ; 
t = 0:T1/100:T1; 

Wo = 2* °/„pi / T ; 

gt = ( 1/2) * ones ( 1 , length (t )) ; 
a(l)=0; / / k =0 , ak =0 

d (1) =0; 
for k =1:5 

a (k + 1 ) = (sin(°/ 0 pi*k/2)/(k*°/ 0 pi)) ; 
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if(abs(a(k+l)) <=0 . 0 1 ) 
a (k+1 ) =0 ; 
end 

d (k + 1 ) = a (k + 1 ) * exp ( - sqrt (-l)*k*°/ 0 pi/2) ; 

end 

k = 0:5 

disp (’ Fourier Series Coefficients of Square Wave’) 
a 

disp( ’Fourier Series Coefficients of g(t) =x ( t — 1 ) — 0 . 5 
’) 
d 

// 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
a. x_location = ’’origin”; 
a. data_bounds = [- l , -2; 1 ,4] ; 
plot2d([-t($:-l:l) , t ( 1 : $) ] ,[-gt,gt] ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 
polyl . thickness = 3; 
title (’ g ( t ) ’) 
xlabel ( ’ 

t ’) 


Scilab code Exa 3.7 Derivative Property of CTFS 

//clear// 

//Example3 . 7 : Derivative Property of CTFS 

clear ; 
clc ; 
close ; 

T =4; 
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T 1 = T/2 ; 

t = 0:T1/100:T1; 

xt = [t ($ : -1 : 1) t] /T1 ; 

gt = ( 1/2) * ones ( 1 , length (t )) ; 

e ( 1 ) = 1/2; //k =0, eO = 1/2 

for k =1:5 

a (k + 1 ) = (sin(°/ 0 pi*k/2) / (k*°/ 0 pi) ) ; 

if(abs(a(k+l)) <=0 . 0 1 ) 
a (k+ 1 ) =0 ; 
end 

d (k+1 ) = a (k+1 ) *exp (-sqrt (-1) *k*%pi/2) ; 
e (k+1 ) = 2*d(k+l) /(sqrt (-1) *k*%pi) ; 
end 

k = 0:5 

disp (’ Fourier Series Coefficients of Square Wave’) 
a 

disp( ’Fourier Series Coefficients of g(t) =x ( t — 1 ) — 0 . 5 
’) 
d 

disp( ’Fourier Series Coefficients of Triangular Wave 
’) 
e 

//Plotting the time shifted square waveform 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a. data_bounds = [- l , -2; 1 ,2] ; 
plot2d([-t($ : -1: 1) , t ( 1 : $) ] ,[-gt,gt] ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 
polyl . thickness = 3; 
title ( ’g(t) ’) 
xlabel ( ’ 

t ’) 

//Plotting the Triangular waveform 

f igure 
a = gca () ; 


81 


41 

42 

43 

44 

45 

46 

47 

48 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 


a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a. data_bounds = [-l ,0; 1 ,2] ; 
plot2d([-t($:-l:l) ,t (1: $)] ,xt ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title ( ’x( t ) ’) 
xlabel ( ’ t ’ ) 

Scilab code Exa 3.7 

Derivative Property of CTFS 

//clear// 

//Example3 . 7 : Derivative Property of CTFS 

clear ; 
clc ; 
close ; 

T =4; 

T 1 = T/2 ; 

t = 0:T1/100:T1; 

xt = [t ( $ : - 1 : 1 ) t ] / T 1 ; 

gt =( 1/2) * ones ( 1 , length (t )) ; 

e ( 1 ) = 1/2; //k =0, eO = 1/2 

for k =1:5 

a (k + 1 ) = (sin(°/ 0 pi*k/2)/(k*°/ 0 pi)) ; 
if(abs(a(k+l)) <=0 . 0 1 ) 
a (k+1 ) =0 ; 
end 

d (k+1 ) = a(k+l) *exp (-sqrt (-1) *k*%pi/2) ; 
e (k + 1 ) = 2*d(k + l) / ( sqrt ( -1 ) *k* °/„pi ) ; 
end 

k = 0:5 

disp (’ Fourier Series Coefficients of Square Wave’) 
a 

disp( ’Fourier Series Coefficients of g(t) =x ( t — 1 ) — 0 . 5 
’) 
d 
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disp( ’Fourier Series Coefficients of Triangular Wave 
’) 
e 

//Plotting the time shifted square waveform 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a. data_bounds = [- l , -2; 1 ,2] ; 
plot2d([-t($ : -1: 1) , t ( 1 : $) ] ,[-gt,gt] ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 
polyl . thickness = 3; 
title( ’g(t) ’) 
xlabel ( ’ 

t ’) 

//Plotting the Triangular waveform 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
a. x_location = ’’origin”; 
a. data_bounds = [- l ,0; 1 ,2] ; 
plot2d([-t($:-l:l) ,t (1: $)] ,xt ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title ( ’x ( t ) ’) 
xlabel ( ’ t ’ ) 


Scilab code Exa 3.8 Fourier Series Representation of Periodic Impulse Train 


//clear// 

// Example3 . 8 : Fourier Series Representation of 
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Periodic Impulse Train 

clear ; 
clc ; 
close ; 

T =4; 

T 1 = T/4; 
t = [-T , 0 , T] ; 

xt = [1,1,1]; //Generation of Periodic train of 
Impulses 

tl = -T1 : Tl/100 : T1 ; 

gt = ones (1 , length (tl ));// Generation of periodic 
square wave 

t2 = [ — T 1 ,0 , Tl] ; 

qt = [1 , 0 , -1] ;// Derivative of periodic square wave 

Wo = 2* °/„pi / T ; 
ak = 1/T ; 
b ( 1 ) = 0; 
c(l) = 2*T1 /T ; 
for k =1:5 

b(k + l) = ak * ( exp ( sqrt ( - 1 ) *k* Wo *T1 ) - exp ( - sqrt ( - 1 ) *k 
*Wo *T1 ) ) ; 

if (abs(b(k+l) ) <=0. 1) 
b (k+1 ) =0; 
end 

c (k+1 ) = b(k+l) /(sqrt (-1) *k*Wo) ; 
if (abs(c(k+l) ) <=0. 1) 
c (k+1 ) =0 ; 
end 
end 

k = 0:5 

disp (’ Fourier Series Coefficients of periodic Square 
Wave ’ ) 
disp (b) 

d i s p ( ’ Fourier Series Coefficients of derivative of 
periodic square wave’) 

di sp ( c ) 

//Plotting the periodic train of impulses 

f igure 
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subplot (3,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

a. data_bounds = [-6 ,0;6 ,2] ; 

plot2d3 ( ’ gnn ’ , t , xt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

title ( ’x( t ) ’) 

//Plotting the periodic square waveform 

subplot (3,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a. data_bounds = [-6 ,0;6 ,2] ; 
plot2d (tl , gt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 
polyl . thickness = 3; 
plot2d(T + tl , gt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 
polyl . thickness = 3; 
plot2d(-T+tl , gt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title ( ’ g ( t ) ’) 

//Plotting the periodic square waveform 

subplot (3,1,3) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

a. data_bounds = [- 6 , — 2 ; 6 ,2] ; 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

plot2d3( ’gnn ’ ,t2 ,qt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

plot2d3( ’gnn ’ , T + t2 , qt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 
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polyl . thickness = 3; 

plot2d3( ’ gnn ’ ,-T + t2,qt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

t it le ( ’ q ( t ) ’) 

Scilab code Exa 3.8 

Fourier Series Representation of Periodic Impulse Train 

//clear// 

// Example3 . 8 : Fourier Series Representation of 
Periodic Impulse Train 
clear ; 
clc ; 
close ; 

T =4; 

T 1 = T/4; 
t = [-T , 0 , T] ; 

xt = [1,1,1]; //Generation of Periodic train of 
Impulses 

tl = -T1 : Tl/100 : T1 ; 

gt = ones (1 , length (tl ));// Generation of periodic 
square wave 

t2 = [ — T 1 ,0 , Tl] ; 

qt = [1 , 0 , -1] ;// Derivative of periodic square wave 

Wo = 2* °/„pi / T ; 
ak = 1/T ; 
b ( 1 ) = 0; 
c ( 1 ) = 2*T1 /T ; 
for k =1:5 

b(k + l) = ak * ( exp ( sqrt ( - 1 ) *k* Wo *T1 ) - exp ( - sqrt ( - 1 ) *k 
*Wo *T1 ) ) ; 

if (abs(b(k+l) ) <=0. 1) 
b (k+1 ) =0; 
end 

c (k+1 ) = b(k+l) /(sqrt (-1) *k*Wo) ; 
if (abs(c(k+l) ) <=0. 1) 
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c (k+1) =0 ; 
end 
end 

k = 0:5 

disp (’ Fourier Series Coefficients of periodic Square 
Wave ’ ) 

disp (b) 

d i s p ( ’Fourier Series Coefficients of derivative of 
periodic square wave’) 
di sp ( c ) 

//Plotting the periodic train of impulses 

f igure 

subplot (3,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 

a.x_location = ’’origin”; 

a. data_bounds = [- 6 ,0;6 ,2] ; 

plot2d3 ( ’gnn ’ ,t , xt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

title ( ’x( t ) ’) 

//Plotting the periodic square waveform 

subplot (3,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
a. data_bounds = [-6 ,0;6 ,2] ; 
plot2d (tl , gt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 
polyl . thickness = 3; 
plot2d(T + tl ,gt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 
polyl . thickness = 3; 
plot2d(-T+tl ,gt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title ( ’ g ( t ) ’) 

//Plotting the periodic square waveform 
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subplot (3,1,3) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

a. data_bounds = [-6 , — 2 ; 6 ,2] ; 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

plot2d3( ’ gnn ’ ,t2 ,qt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

plot2d3( ’gnn ’ ,T+t2,qt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

plot2d3( ’gnn ’ ,-T + t2,qt ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

title(’q(t) ’) 


Scilab code Exa 3.10 Scilab code Exa 3.10 DTFS of x[n] 

//clear// 

// Example3 . 1 0 : DTFS of x[n] =sin (Won) 

clear ; 
close ; 
clc ; 

n = 0:0.01:5; 

N = 5; 

Wo = 2* °/„pi / N ; 
xn = sin (Wo*n) ; 
for k =0 : N -2 

C ( k + 1 , : ) = exp(-sqrt(-l)*Wo*n.*k) ; 
a (k + 1 ) = xn*C(k + l ,:) ’/length (n) ; 
if(abs(a(k+l)) <=0 . 0 1 ) 
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a (k+ 1 ) =0 ; 
end 
end 
a =a J 

a_con j = con j ( a) ; 

ak = [a_conj($:-l:l),a(2:$)] 

k = -(N-2) : (N-2) ; 

// 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

a.data_bounds=[-8,-l;8,l] ; 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

plot2d3( ’ gnn ’ , k , - imag ( ak ) ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

plot2d3( ’gnn ’ ,N+k,-imag(ak) ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

plot 2d3 (’ gnn ’ , - (N + k) ,-imag(ak($:-l:l)) ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title (’ ak ’ ) 

DTFS of x[n| 

//clear// 

// Example3 . 1 0 : DTFS of x[n] =sin (Won) 

clear ; 
close ; 
clc ; 

n = 0:0.01:5; 

N = 5; 

Wo = 2* °/„pi / N ; 
xn = sin (Wo*n) ; 
for k =0 : N -2 
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C(k+1,:) = exp ( - sqrt ( - 1 ) * Wo *n . *k) ; 
a(k+l) = xn*C (k+1 length (n) ; 
if(abs(a(k+l)) <=0 . 0 1 ) 
a (k+1 ) =0 ; 


end 
end 
a = a ’ 

a_con j = conj ( a) ; 

ak = [a_conj($:-l:l),a(2:$)] 

k = -(N-2) : (N-2) ; 


f igure 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

a.data_bounds=[-8,-l;8,l] ; 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

plot2d3( ’ gnn ’ , k , - imag ( ak ) ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

plot2d3( ’gnn ’ ,N+k,-imag(ak) ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

plot2d3(’gnn’ , - (N + k) ,-imag(ak($:-l:l)) ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title (’ ak ’ ) 


Scilab code Exa 3.11 Scilab code Exa 3.11 DTFS of x[n] 

//clear// 

// Example3 . 1 1 : DTFS of 
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//x[n] = 1+sin (2* %pi/N) n+3cos (2 * %pi/N) n+cos [(4* %pi/N 
) n+%pi / 2] 

clear ; 
close ; 
clc ; 

N = 10; 
n = 0 : 0 . 01 : N ; 

Wo = 2*°/„pi/N; 

xn =ones (1 , length(n))+sin(Wo*n)+3*cos(Wo*n)+cos (2*Wo 
*n + 7„pi /2) ; 
for k =0 : N -2 

C ( k+ 1 , : ) = exp (-sqrt (-1) *Wo*n . *k) ; 
a(k + l) = xn*C (k + 1 ,:) J / length (n) ; 
if (abs(a(k+l) ) <=0 . 1) 
a (k+1 ) =0 ; 
end 
end 

a =a ’ ; 

a_conj =conj (a) ; 

ak = [a_conj($:-l:l),a(2:$)]; 

Mag_ak = abs(ak); 
for i = l:length(a) 

Phase_ak(i) = atan ( imag ( ak ( i ) ) / ( real ( ak ( i ) ) 

+ 0 . 0001 ) ) ; 

end 

Phase_ak = Phase_ak ’ 

Phase_ak = [Phase_ak ( 1 : $-1) -Phase_ak ($ : -1 : 1) ] ; 
k = -(N-2) : (N-2) ; 

// 

f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
plot2d3( ’ gnn ’ ,k,real(ak) ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title(’Real part of ( ak ) ’ ) 
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xlabel ( ’ 


k’) 

subplot (2,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

plot2d3( ’ gnn ’ ,k, imag(ak) ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

title ( ’imaginary part of(ak)’) 

xlabel ( ’ 

k’) 

// 

f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x_location = ’’origin”; 

plot2d3( ’gnn ’ ,k,Mag_ak ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

title ( ’ abs ( ak) ’ ) 

xlabel ( ’ 

k’) 

subplot (2,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

plot2d3( ’gnn ’ ,k,Phase_ak ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

title ( ’ <(ak ) ’ ) 

xlabel ( ’ 

k’) 
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DTFS of x[n| 


//clear// 

// Example3 . 1 1 : DTFS of 

//x[n] = 1+sin (2* %pi/N) n+3cos (2* %pi/N) n+cos [ (4 * %pi/N 
)n+%pi/2] 

clear ; 
close ; 
clc ; 

N = 10; 
n = 0 : 0 . 01 : N ; 

Wo = 2* °/„pi / N ; 

xn =ones (1 , length(n))+sin(Wo*n)+3*cos(Wo*n)+cos (2* Wo 
*n + °/„pi /2) ; 
for k =0 : N -2 

C(k+1,:) = exp ( - sqrt ( - 1 ) * Wo *n . *k) ; 
a(k + l) = xn*C (k + 1 ,:) ’ / length (n) ; 
if (abs(a(k+l) ) <=0. 1) 
a (k+1 ) =0 ; 
end 
end 

a =a ’ ; 

a_conj =conj (a) ; 

ak = [a_conj($:-l:l),a(2:$)]; 

Mag_ak = abs(ak); 
for i = l:length(a) 

Phase_ak(i) = atan ( imag ( ak ( i ) ) / ( real ( ak ( i ) ) 

+ 0 . 0001 ) ) ; 

end 

Phase_ak = Phase_ak ’ 

Phase_ak = [Phase_ak ( 1 : $-1) -Phase_ak ($ : -1 : 1) ] ; 
k = -(N-2) : (N-2) ; 

// 

f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 
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a.x_location = ’’origin”; 
plot2d3( ’ gnn ’ ,k,real(ak) ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title(’Real part of(ak)’) 
xlabel ( ’ 

k’) 

subplot (2,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

plot2d3( ’gnn ’ ,k, imag(ak) ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

title ( ’imaginary part of(ak)’) 

xlabel ( ’ 

k’) 

// 

f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x_location = ’’origin”; 

plot2d3 ( ’gnn ’ , k , Mag_ak ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

title ( ’ abs ( ak) ’ ) 

xlabel ( ’ 

k’) 

subplot (2,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
plot2d3( ’gnn ’ ,k,Phase_ak ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 


94 


65 

66 

67 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 


polyl . thickness = 3; 
title ( ’ <(ak ) ’ ) 
xlabel ( ’ 

k’) 


Scilab code Exa 3.12 DTFS coefficients of periodic square wave 

//clear// 

// Example3 . 1 2 : DTFS coefficients of periodic square 
wave 

clear ; 
close ; 
clc ; 

N = 10; 

N 1 = 2; 

Wo = 2* °/„pi / N ; 
xn = ones (1 , length (N) ) ; 
n = -(2*N1+1) : (2*N1+1) ; 
a ( 1 ) = (2*N1+1)/N; 
for k = 1 : 2* N1 

a(k + l) = sin((2* 0 / 0 pi*k*(Nl+0.5))/N)/sin( 0 / 0 pi*k/N); 
a (k+1 ) = a (k+1 ) /N ; 
if (abs(a(k+l) ) <=0. 1) 
a (k+ 1 ) =0; 
end 
end 

a =a ’ ; 

a_conj =conj (a) ; 

ak = [a_conj($:-l:l),a(2:$)]; 

k = -2*N1 : 2* N1 ; 

// 

f igure 
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a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
plot2d3( ’ gnn ’ ,k,real(ak) ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
title(’Real part of ( ak ) ’ ) 
xlabel ( ’ 

k’) 

Scilab code Exa 3.12 

DTFS coefficients of periodic square wave 

//clear// 

// Example3 . 1 2 : DTFS coefficients of periodic square 
wave 

clear ; 
close ; 
clc ; 

N = 10; 

N 1 = 2; 

Wo = 2* °/„pi / N ; 
xn = ones ( 1 , length (N) ) ; 
n = -(2*N1 + 1) : (2*N1+1) ; 
a ( 1 ) = (2*N1+1)/N; 
for k = 1 : 2* N1 

a (k + 1 ) = sin((2*°/ 0 pi*k*(Nl+0.5) )/N)/sin(°/ 0 pi*k/N) ; 
a (k+1 ) = a (k+1 ) /N ; 
if (abs(a(k+l) ) <=0 . 1) 
a (k+ 1 ) =0; 
end 
end 

a =a ’ ; 

a_conj =conj (a) ; 

ak = [a_conj($:-l:l),a(2:$)]; 

k = -2*N1 : 2* N1 ; 
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// 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
plot2d3( ’ gnn ’ ,k,real(ak) ,5) 
polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 
polyl . thickness = 3; 
title(’Real part of ( ak ) ’ ) 
xlabel ( ’ 

k’) 


Scilab code Exa 3.13 Periodic sequence 

//clear// 

/ / Example3 .13: DTFS 

//Expression of periodic sequence using 
//the summation two different sequence 

clear ; 
close ; 
clc ; 

N = 5; 

n = 0 : N - 1 ; 

xl = [1 , 1 , 0 ,0 ,1] ; 

xl = [xl ( $ : -1 : 1 ) xl (2 : $ ) ] ; // Square Wave xl[n] 

x2 = [1,1,1,1,1] ; 

x2 = [x2($:-l:l) x2 (2 : $ ) ] ; //DC sequence of x2[n] 
x = xl+x2;//sum of xl [n] & x2[n] 

//Zeroth DTFS coefficient of dc sequence 

c(l) = 1; 

//Zeroth DTFS coefficient of square waveform 

b ( 1 ) = 3/5; 
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//Zeroth DTFS coefficient of sum of xl [n] & x2[n] 

a ( 1 ) = b(l)+c(l) ; 

it 

Wo = 2*%pi/N; 
for k =1 : N-l 

a (k + 1 ) = sin((3* , /.pi*k)/N)/sin( , /.pi*k/N) ; 
a (k+1 ) = a (k+1 ) /N ; 
if (abs(a(k+l) ) <=0. 1) 
a (k+ 1 ) =0; 
end 
end 

a =a ’ ; 

a_conj =conj (a) ; 
ak = [a_conj($:-l:l),a(2:$)]; 
k = -(N-l) : (N-l) ; 
n = -(N-l) : (N-l) ; 

// 

f igure 

subplot (3,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 

a.x_location = ’’origin”; 

plot2d3( ’ gnn ’ , n , x , 5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

title ( ’x [n] ’ ) 

xlabel ( ’ 

n ’ ) 

subplot (3,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

plot2d3( ’gnn ’ ,n,xl ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

title ( ’xl [n] ’) 

xlabel ( ’ 
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n ’ ) 

subplot (3,1,3) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

plot2d3( ’ gnn ’ ,n,x2 ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

title(’x2[n] ’) 

xlabel ( ’ 

n’) 

Scilab code Exa 3.13 

Periodic sequence 

//clear// 

/ / Example3 .13: DTFS 

//Expression of periodic sequence using 
//the summation two different sequence 

clear ; 
close ; 
clc ; 

N = 5; 

n = 0 : N - 1 ; 

xl = [1 , 1 ,0 ,0 , 1] ; 

xl = [xl ( $ : - 1 : 1 ) xl (2 : $ ) ] ; // Square Wave xl[n] 

x2 = [1 , 1 , 1 , 1 , 1] ; 

x2 = [x2($:-l:l) x2 (2 : $ ) ] ; //DC sequence of x2[n] 
x = xl+x2;//sum of xl [n] & x2[n] 

//Zeroth DTFS coefficient of dc sequence 

c(l) = 1; 

//Zeroth DTFS coefficient of square waveform 

b ( 1 ) = 3/5; 

//Zeroth DTFS coefficient of sum of xl [n] & x2[n] 

a ( 1 ) = b(l)+c(l) ; 
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// 

Wo = 2*%pi/N; 
for k =1 : N-l 

a(k + l) = sin((3*°/ 0 pi*k)/N)/sin(°/„pi*k/N); 

a(k+l) = a(k+l)/N; 
if (abs(a(k+l) ) <=0. 1) 
a (k+ 1 ) =0; 
end 
end 

a =a ’ ; 

a_conj =conj (a) ; 
ak = [a_conj($:-l:l),a(2:$)]; 
k = -(N-l) : (N-l) ; 
n = -(N-l) : (N-l) ; 

// 

f igure 

subplot (3,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x_location = ’’origin”; 

plot2d3( ’ gnn ’ , n , x , 5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

title ( ’x [n] ’ ) 

xlabel ( ’ 

n ’ ) 

subplot (3,1,2) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

plot2d3( ’gnn ’ ,n,xl ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

title ( ’xl [n] ’) 

xlabel ( ’ 

n ’) 
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subplot (3,1,3) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

plot2d3 ( ’ gnn ’ , n , x2 ,5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

title(’x2[n] ’) 

xlabel ( ’ 

n ’) 


Scilab code Exa 3.14 Parseval’s relation of DTFS 

//clear// 

/ / Example3 . 1 4 : DTFS 

//Finding x [ n ] using parseval’s relation of DTFS 
clear ; 
close ; 
clc ; 

N = 6; 
n = 0 : N - 1 ; 
a ( 1 ) = 1/3; 
a (2) =0; 
a (4) =0; 
a (5) =0; 

al = (1/6) *((-1) ~n) ; 
x =0 ; 

for k = 0 : N -2 
if (k==2) 
x = x+al ; 
else 

x = x+a (k+1 ) ; 
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end 

end 

x = [x($ : -1 : 1) ,x(2: $)] ; 
n = - (N - 1 ) : (N-l) ; 

// 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

plot2d3( ’ gnn ’ , n , x , 5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

title ( ’x [n] ’ ) 

xlabel ( ’ 

n’) 

Scilab code Exa 3.14 

Parseval’s relation of DTFS 

//clear// 

/ / Example3 .14: DTFS 

//Finding x [ n ] using parseval’s relation of DTFS 
clear ; 
close ; 
clc ; 

N = 6; 
n = 0 : N - 1 ; 
a ( 1 ) = 1/3; 
a (2) =0; 
a (4) =0; 
a (5) =0; 

al = (1/6) *((-1) ~n) ; 
x =0 ; 

for k = 0 : N -2 
if (k==2) 
x = x+al ; 
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else 

x = x+a (k+1 ) ; 
end 
end 

x = [x($ : -1 : 1) ,x(2: $)] ; 
n = - (N - 1 ) : (N-l) ; 

// 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

plot2d3( ’ gnn ’ , n , x , 5) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

title ( ’x [n] ’ ) 

xlabel ( ’ 

n ’) 


Scilab code Exa 3.15 DTFS:Periodic Convolution Property 

//clear// 

//Example3 . 1 5 : DTFS : Per iodic Convolution Property 

clear ; 
clc ; 
close ; 

x = [1 , 1 ,0 ,0 ,0 ,0 , 1] ; 

X = f ft (x) ; 

w = x . *x ; 

w = if f t (W) ; 
w = abs (w) ; 
for i =l:length(x) 
if ( abs ( w ( i ) ) < = 0 . 1 ) 
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w ( i ) = 0 ; 
end 

end 

w = [w($:— 1:1) w ( 2 : $ ) ] ; 

N = length (x) ; 
f igure 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

plot2d3( ’ gnn ’ , [ — ( N — 1 ) : 0 , 1 : N — 1 ] , w , 5) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

title(’w[n] ’) 

xlabel ( ’ 

n’) 

Scilab code Exa 3.15 

DTFS: Periodic Convolution Property 

//clear// 

//Example3 . 1 5 : DTFS : Per iodic Convolution Property 

clear ; 
clc ; 
close ; 

x = [1 , 1 ,0 ,0 ,0 ,0 , 1] ; 

X = f ft (x) ; 

w = x . *x ; 

w = if f t (W) ; 
w = abs (w) ; 
for i =l:length(x) 
if ( abs ( w ( i ) ) < = 0 . 1 ) 
w ( i ) = 0 ; 
end 

end 

w = [w($:— 1:1) w(2:$)]; 

N = length (x) ; 
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18 figure 

19 a = gca () ; 

20 a.y_location = ’’origin”; 

21 a. x.location = ’’origin”; 

22 plot2d3 ( ’ gnn ’ , [ - ( N - 1 ) : 0 , 1 : N-l] , w ,5) 

23 polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

24 polyl . thickness = 3; 

25 title(’w[n]’) 

26 xlabel ( ’ 

n ’) 
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Chapter 8 

The Cotntinuous Time Fourier 
Transform 


Scilab code Exa 4.1 clear 

//clear// 

//Example 4 . 1 : Continuous Time Fourier Transform of a 
//Continuous Time Signal x(t)= exp(— A* t ) u ( t ) , t >0 

clear ; 
clc ; 
close ; 

// Analog Signal 
A =1; //Amplitude 

Dt = 0.005; 
t = 0 : Dt : 10 ; 
xt = exp ( - A*t ) ; 

// 

// Continuous —time Fourier Transform 

Wraax = 2*%pi*l; //Analog Frequency = 1Hz 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 
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XW = xt* exp (- sqrt ( -1) *t ’ *W) * Dt ; 

XW_Mag = abs(XW); 

W = [-mtlb_fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 

XW.Mag = [mtlb.f liplr (XW.Mag) , XW.Mag (2 : 1001) ] ; 
[XW_Phase , db] = phaseraag (XW) ; 

XW.Phase = [-mtlb_f liplr (XW.Phase ), XW_Phase (2 : 1001) 

] ; 

//Plotting Continuous Time Signal 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 

plot (t , xt ) ; 

xlabel ( ’ t in sec.’); 


ylabel ( ’x ( t ) ’ ) 

title (’ Continuous Time Signal’) 
f igure 

//Plotting Magnitude Response of CTS 

subplot (2,1.1) ; 
a = gca () ; 

a.y_location = ’’origin”; 
plot (W , XW.Mag ) ; 

xlabel (’ Frequency in Radians / Seconds > W’ ) ; 

ylabel (’ abs (X(jW) ) ’ ) 

title (’ Magnitude Response (CTFT) ’ ) 

//Plotting Phase Reponse of CTS 

subplot (2,1,2); 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
plot(W,XW_Phase*°/ 0 pi/180) ; 

xlabel (’ Frequency in 

Radians / Seconds > W’ ) ; 

ylabel ( ’ 


<X 


(jW)’) 

title (’Phase Response (CTFT) in Radians’) 
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Scilab code Exa 4.1 

clear 

//clear// 

//Example 4 . 1 : Continuous Time Fourier Transform of a 
//Continuous Time Signal x(t)= exp(— A* t ) u ( t ) , t >0 

clear ; 
clc ; 
close ; 

// Analog Signal 
A =1; //Amplitude 

Dt = 0.005; 
t = 0 : Dt : 10 ; 
xt = exp ( - A*t ) ; 

// 

// Continuous —time Fourier Transform 

Wmax = 2* °/„pi * 1 ; //Analog Frequency = 1Hz 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 

XW = xt* exp (- sqrt ( -1) *t ’ *W) * Dt ; 

XW.Mag = abs (XW) ; 

W = [-mtlb.fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 

XW.Mag = [mtlb.f liplr (XW.Mag) , XW_Mag (2 : 1001) ] ; 
[XW_Phase , db] = phaseraag (XW) ; 

XW_Phase = [-mtlb_f liplr (XW_Phase ), XW_Phase (2 : 1001) 

] ; 

//Plotting Continuous Time Signal 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
plot (t , xt ) ; 
x label (’ t in sec.’); 
ylabel ( ’x ( t ) ’ ) 

title (’ Continuous Time Signal’) 
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f igure 

//Plotting Magnitude Response of CTS 

subplot (2,1,1) ; 
a = gca () ; 

a.y_location = ’’origin”; 
plot (W , XW_Mag ) ; 

xlabel (’ Frequency in Radians / Seconds > W’ ) ; 

ylabel ( ’ abs (X(jW) ) ’ ) 

title (’ Magnitude Response (CTFT) ’ ) 

//Plotting Phase Reponse of CTS 

subplot (2,1,2); 
a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
plot (W,XW_Phase*°/ 0 p i/180) ; 

xlabel ( ’ Frequency in 

Radians / Seconds > W’ ) ; 

ylabel ( ’ 

(jW) ’ ) 

title (’Phase Response (CTFT) in Radians’) 


<X 


Scilab code Exa 4.2 clear 

//clear// 

//Example 4 . 2 : Continuous Time Fourier Transform of a 
//Continuous Time Signal x(t)= exp (— A* abs ( t ) ) 

clear ; 
clc ; 
close ; 

// Analog Signal 
A =1; //Amplitude 

Dt = 0.005; 
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t = -4 . 5 : Dt : 4 . 5 ; 
xt = exp ( - A* abs ( t ) ) ; 

// 

// Continuous —time Fourier Transform 

Wmax = 2* °/„pi * 1 ; //Analog Frequency = 1Hz 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 

XW = xt* exp (- sqrt ( -1) *t ’ *W) * Dt ; 

XW = real (XW) ; 

W = [-mtlb_fliplr(W), W(2:1001)]; // Omega from 
Wmax t o Wmax 

XW = [mtlb.f liplr (XW) , XW(2:1001)]; 
subplot (1,1,1) 
subplot (2,1,1) ; 
a = gca () ; 

a.y_location = ’’origin”; 
plot (t , xt ) ; 
x label (’ t in sec.’); 
ylabel ( ’x ( t ) ’ ) 

title (’ Continuous Time Signal’) 
subplot (2,1,2); 
a = gca () ; 

a.y_location = ’’origin”; 
plot (W , XW) ; 

xlabel (’ Frequency in Radians / Seconds W’ ) ; 
ylabel ( ’X(jW) ’) 

title (’ Continuous— time Fourier Transform’) 

Scilab code Exa 4.2 

clear 

//clear// 

//Example 4 . 2 : Continuous Time Fourier Transform of a 
//Continuous Time Signal x(t)= exp (— A* abs ( t ) ) 

clear ; 
clc ; 
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close ; 

// Analog Signal 
A =1; //Amplitude 

Dt = 0.005; 

t = -4 . 5 : Dt : 4 . 5 ; 

xt = exp ( - A* abs ( t ) ) ; 

// 

// Continuous —time Fourier Transform 

Wraax = 2*%pi*l; //Analog Frequency = 1Hz 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 

XW = xt* exp ( - sqrt ( - 1 ) *t 5 * W) * Dt ; 

XW = real (XW) ; 

W = [-mtlb_fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 

XW = [mtlb_f liplr (XW) , XW(2:1001)]; 
subplot (1,1,1) 
subplot (2,1,1) ; 
a = gca () ; 

a.y_location = ’’origin”; 
plot (t , xt ) ; 
xlabel ( ’ t in sec.’); 
ylabel ( ’x ( t ) ’ ) 

title (’ Continuous Time Signal’) 
subplot (2,1,2); 
a = gca () ; 

a.y_location = ’’origin”; 
plot (W , XW) ; 

xlabel (’ Frequency in Radians / Seconds W’ ) ; 
ylabel ( ’X(jW) ’) 

title (’ Continuous— time Fourier Transform’) 
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Scilab code Exa 4.4 clear 


//clear// 

//Example 4.4 

// Continuous Time Fourier Transform 

//and Frequency Response of a Square Waveform 

// x(t)= A, from — T1 to T1 

clear ; 
clc ; 
close ; 

// CTS Signal 
A =1; //Amplitude 
Dt = 0.005; 

T1 = 4; //Time in seconds 

t = -Tl/2: Dt :Tl/2; 
for i = l:length(t) 
xt ( i ) = A ; 
end 

// 

// Continuous —time Fourier Transform 

Wraax = 2*%pi*l; //Analog Frequency = 1Hz 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 
xt = xt 1 ; 

XW = xt* exp ( - sqrt ( - 1 ) *t 1 * W) * Dt ; 

XW.Mag = real (XW) ; 

W = [-mtlb_fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 

XW.Mag = [mtlb.f liplr (XW_Mag) , XW_Mag (2 : 1001) ] ; 

// 

subplot (2,1,1) ; 
a = gca () ; 

a.data_bounds=[-4,0;4,2] ; 

a.y_location =” origin”; 

plot (t , xt ) ; 

xlabel ( ’ t in msec . ’ ) ; 

title (’ Contiuous Time Signal x(t)’) 
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subplot (2,1,2); 
a = gca () ; 

a.y_location =” origin”; 
plot (W , XW_Mag ) ; 

xlabel (’ Frequency in Radians / Seconds ’) ; 

title (’ Continuous— time Fourier Transform X(jW) ’) 

Scilab code Exa 4.4 

clear 

//clear// 

//Example 4.4 

// Continuous Time Fourier Transform 

//and Frequency Response of a Square Waveform 

// x(t)= A, from — T1 to T1 

clear ; 
clc ; 
close ; 

// CTS Signal 
A =1; //Amplitude 

Dt = 0.005; 

T1 = 4; //Time in seconds 

t = -Tl/2: Dt :Tl/2; 
for i = l:length(t) 
xt ( i ) = A ; 
end 

// 

// Continuous —time Fourier Transform 

Wraax = 2*%pi*l; //Analog Frequency = 1Hz 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 
xt = xt ’ ; 

XW = xt* exp ( - sqrt ( - 1 ) *t 1 * W) * Dt ; 

XW.Mag = real (XW) ; 

W = [-mtlb_fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 
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XW.Mag = [mtlb.f liplr (XW_Mag) , XW_Mag (2 : 1001) ] ; 

// 

subplot (2,1,1) ; 
a = gca () ; 

a.data_bounds=[-4,0;4,2] ; 

a.y_location = ” origin”; 

plot (t , xt ) ; 

xlabel ( ’ t in msec . ’ ) ; 

title (’ Contiuous Time Signal x(t)’) 

subplot (2,1,2); 

a = gca () ; 

a.y_location =” origin”; 
plot (W , XW_Mag ) ; 

xlabel (’ Frequency in Radians / Seconds ’) ; 

title (’ Continuous— time Fourier Transform X(jW) ’) 


Scilab code Exa 4.5 clear 

//clear// 

//Example 4.5 

// Inverse Continuous Time Fourier Transform 
// X(jW)= 1, from — T1 to T1 

clear ; 
clc ; 
close ; 

// CTFT 

A =1; //Amplitude 

Dw = 0.005; 

W1 = 4; //Time in seconds 

w = -Wl/2:Dw:Wl/2; 
for i = l:length(w) 

XW ( i ) = A; 
end 
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XW = XW 5 ; 

// 

//Inverse Continuous— time Fourier Transform 

t = -%pi : %pi/ length (w) : %pi ; 

xt = ( 1/ (2* %pi ) ) *XW * exp ( sqrt ( - 1 ) *w ’ * t ) *Dw ; 
xt = real ( xt ) ; 
f igure 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot (t , xt ) ; 

x label ( ’ t time 

in Seconds 7 ) ; 

title (' Inverse Continuous Time Fourier Transform x(t 
) ’) 


Scilab code Exa 4.5 

clear 

//clear// 

//Example 4.5 

// Inverse Continuous Time Fourier Transform 
// X(jW)= 1, from — T1 to T1 

clear ; 
clc ; 
close ; 

// CTFT 

A =1; //Amplitude 

Dw = 0.005; 

W1 = 4; //Time in seconds 

w = -Wl/2:Dw:Wl/2; 
for i = l:length(w) 

XW(i) = A; 
end 

XW = XW ’ ; 

// 

//Inverse Continuous —time Fourier Transform 
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t = -%pi : %pi/ length (w) : °/„pi ; 

xt = ( 1/ (2* %pi ) ) *XW * exp ( sqrt ( - 1 ) *w ’ * t ) *Dw ; 
xt = real ( xt ) ; 
f igure 
a = gca () ; 

a.y_location = ” origin”; 
a. x.location =” origin”; 
plot (t , xt ) ; 

x label ( ’ t time 

in Seconds ’ ) ; 

title (' Inverse Continuous Time Fourier Transform x(t 
) ’) 


Scilab code Exa 4.6 clear 

//clear// 

//Example 4.6 

// Continuous Time Foirrier Transform of Symmetric 
// periodic Square waveform 

clear ; 
clc ; 
close ; 

// CTFT 
T 1 = 2; 

T = 4* T1 ; 

Wo = 2* °/„pi / T ; 

W = -%pi : Wo : °/ 0 pi ; 

delta = ones (1 , length (W) ) ; 

XW(1) = ( 2 * °/„p i * Wo * T 1 / %p i ) ; 
raid_value = ceil ( length (W) /2) ; 
for k = 2:mid_value 

XW(k) = (2* °/ 0 pi * sin ( (k- 1 ) * Wo *T1) / ( °/ 0 pi * (k-1 ) ) ) ; 
end 
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f igure 
a = gca () ; 

a.y_location = ” origin”; 

a. x.location =” origin”; 

plot2d3( ’ gnn ’ ,W(mid_value :$) ,XW,2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

plot2d3( ’gnn’ ,W(l:mid_value-l) , XW ( $ : -1 : 2) ,2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

xlabel ( ’W in radians / Seconds ’) ; 

title (’ Continuous Time Fourier Transform of Periodic 
Square Wave’) 

Scilab code Exa 4.6 

clear 

//clear// 

//Example 4.6 

// Continuous Time Foirrier Transform of Symmetric 
// periodic Square waveform 

clear ; 
clc ; 
close ; 

// CTFT 
T 1 = 2; 

T = 4* T1 ; 

Wo = 2* °/„pi / T ; 

W = -%pi : Wo : °/ 0 pi ; 

delta = ones (1 , length (W) ) ; 

XW(1) = ( 2 * %p i * Wo * T 1 / %p i ) ; 
raid_value = ceil ( length (W) /2) ; 
for k = 2:mid_value 

XW(k) = (2* °/ 0 pi * sin ( (k-1) *Wo *T1 ) / ( °/ 0 pi * (k-1) ) ) ; 

end 

f igure 
a = gca () ; 
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a.y_location =” origin”; 

a. x.location =” origin”; 

plot2d3( ’ gnn ’ ,W(mid_value :$) ,XW,2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

plot2d3( , gnn’,W(l:mid_value-l) , XW ( $ : -1 : 2) ,2); 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

xlabel ( ’W in radians / Seconds ’) ; 

title (’ Continuous Time Fourier Transform of Periodic 
Square Wave’) 


Scilab code Exa 4.7 clear 

//clear// 

//Example 4.7 

// Continuous Time Fourier Transforms of 
// Sinusoidal waveforms (a) sin (Wot) (b) cos (Wot) 

clear ; 
clc ; 
close ; 

// CTFT 
T 1 = 2; 

T = 4* T1 ; 

Wo = 2* °/„pi / T ; 

W = [-Wo ,0 , Wo] ; 

ak = (2* Zpi * Wo *T1 /°/«pi ) / sqrt ( - 1 ) ; 

XW = [-ak ,0 , ak] ; 

akl = (2* %pi * Wo *T1 /%pi ) ; 

XW1 = [akl ,0, akl] ; 

// 

f igure 
a = gca () ; 
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a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d3( ’ gnn ’ , W , imag(XW) ,2) ; 
polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 
xlabel ( ’ 

W’) ; 

title ( ’CTFT of sin (Wot) ’) 

// 

f igure 
a = gca () ; 

a.y_location =” origin”; 

a.x_location =” origin”; 

plot2d3( ’gnn ’ , W , XW1 ,2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

xlabel ( ’ 

W’) ; 

title (’CTFT of cos (Wot)’) 

Scilab code Exa 4.7 

clear 

//clear// 

//Example 4.7 

// Continuous Time Fourier Transforms of 
// Sinusoidal waveforms (a) sin (Wot) (b) cos (Wot) 

clear ; 
clc ; 
close ; 

// CTFT 
T 1 = 2; 

T = 4* T1 ; 

Wo = 2* °/„pi / T ; 

W = [-Wo ,0 , Wo] ; 
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ak = (2* Zpi * Wo *T1 /%pi ) / sqrt ( - 1 ) ; 

XW = [-ak ,0 , ak] ; 

akl = (2* %pi * Wo *T1 /%pi ) ; 

XW1 = [akl ,0 , akl] ; 

// 

f igure 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d3( ’ gnn ’ , W , imag(XW) ,2) ; 
polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 
polyl . thickness = 3; 
xlabel ( ’ 

W’) ; 

title (’CTFT of sin (Wot)’) 

// 

f igure 
a = gca () ; 

a.y_location =” origin”; 

a. x.location =” origin”; 

plot2d3( ’gnn ’ , W , XW1 ,2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

xlabel ( ’ 

W’) ; 

title (’CTFT of cos (Wot)’) 


Scilab code Exa 4.8 clear 

//clear// 

//Example 4.8 
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// Continuous Time Fourier Transforms of 

/ / Periodic Impulse Train 

clear ; 

clc ; 

close ; 

// CTFT 
T = -4:4; ; 

T1 = 1; //Sampling Interval 

xt = ones ( 1 , length (T) ) ; 
ak = 1/T1 ; 

XW = 2*°/„pi *ak*ones (1 , length (T) ) ; 

Wo = 2* %pi /T1 ; 

W = Wo *T ; 
f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location =” origin”; 

a.x_location =” origin”; 

plot2d3( ’ gnn ’ ,T,xt ,2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

xlabel ( ’ 

t ’) ; 

t i 1 1 e ( ’ Periodic Impulse Train’) 
subplot (2,1,2) 
a = gca () ; 

a.y_location =” origin”; 

a. x.location =” origin”; 

plot2d3( ’gnn ’ , W , XW , 2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

xlabel ( ’ 

t ’) ; 

title (’CTFT of Periodic Impulse Train’) 


Scilab code Exa 4.8 clear 
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//clear// 

//Example 4.8 

// Continuous Time Fourier Transforms of 

/ / Periodic Impulse Train 

clear ; 

clc ; 

close ; 

// CTFT 
T = -4:4; ; 

T1 = 1; //Sampling Interval 

xt = ones ( 1 , length (T) ) ; 
ak = 1/T1 ; 

XW = 2* °/„pi * ak* ones ( 1 , length (T) ) ; 

Wo = 2* %pi /T1 ; 

W = Wo *T ; 
f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location =” origin”; 

a. x_location =” origin”; 

plot2d3( ’gnn ’ ,T,xt ,2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

xlabel ( ’ 

t ’) ; 

t i 1 1 e ( ’ Periodic Impulse Train’) 
subplot (2,1,2) 
a = gca () ; 

a.y_location =” origin”; 

a. x.location =” origin”; 

plot2d3( ’gnn ’ , W , XW , 2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

xlabel ( ’ 

t ’) ; 

title (’CTFT of Periodic Impulse Train’) 
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Scilab code Exa 4.9 clear 

//clear// 

//Example 4 . 9 : Continuous Time Fourier Transform 
Properties : 

//Linearity and Time Shift Property 
clear ; 
clc ; 
close ; 

// CTFT 

tl = -1/2 : 0 . 1 : 1/2 ; 

t2 = -3/2 : 0 . 1 : 3/2 ; 

xl = ones (1 , length (tl )) ; 

x2 = ones (1 , length (t2) ) ; 

t3 = t 1 +2 . 5 ; 

t4 = t2+2.5; 

xl = (1/2) *xl ; 

x = [x2 ( 1 : f loor ( length ( x2 ) /3) ), xl +x2 ( ce i 1 ( length ( x2 ) 

/3) :$-floor(length(x2)/3)) ,x2(($-ceil(length(x2) 
/3))+2:$)] ; 
subplot (3,1,1) 
a = gca () ; 

a. x.location = ’’origin”; 
a.y_location = ’’origin”; 
plot (tl , xl ) 
xt it le ( ’ xl ( t ) ’) 
subplot (3,1,2) 
a = gca () ; 

a.x_location = ’’origin”; 
a.y_location = ’’origin”; 
plot (t2 , x2 ) 
xt it le ( ’ x2 ( t ) ’) 
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subplot (3,1,3) 
a = gca () ; 

a. x.location = ’’origin”; 
a.y_location = ’’origin”; 
plot (t4 , x) 
xt it le ( ’ x ( t ) ’) 

Scilab code Exa 4.9 

clear 

//clear// 

//Example 4 . 9 : Continuous Time Fourier Transform 
Properties : 

//Linearity and Time Shift Property 

clear ; 
clc ; 
close ; 

// CTFT 

tl = -1/2 : 0 . 1 : 1/2 ; 

t2 = -3/2 : 0 . 1 : 3/2 ; 

xl = ones (1 , length (tl )) ; 

x2 = ones (1 , length (t2) ) ; 

t3 = tl+2.5; 

t4 = t2+2.5; 

xl = (1/2) *xl ; 

x = [x2 ( 1 : f loor ( length ( x2 ) /3) ), xl +x2 ( ce i 1 ( length ( x2 ) 

/3) :$-floor(length(x2)/3)) ,x2(($-ceil(length(x2) 
/3))+2:$)] ; 
subplot (3,1,1) 
a = gca () ; 

a.x_location = ’’origin”; 
a.y_location = ’’origin”; 
plot (tl , xl ) 
xt it le ( ’ xl ( t ) ’) 
subplot (3,1,2) 
a = gca () ; 

a. x.location = ’’origin”; 
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a.y_location = ’’origin”; 
plot (t2 , x2 ) 
xt it le ( ’ x2 ( t ) ’) 
subplot (3,1,3) 
a = gca () ; 

a. x.location = ’’origin”; 
a.y_location = ’’origin”; 
plot (t4 , x) 
xt it le ( ’ x ( t ) ’) 


Scilab code Exa 4.12 clear 

//clear// 

//Example 4 . 1 2 : Continuous Time Fourier Transform: 
//Derivative property 

clear ; 
clc ; 
close ; 

// CTFT 

t = -1 : 0 . 1 : 1 ; 

xl = ones ( 1 , length (t )) ; 

x2 = [-l,zeros(l,length(t)-2),-l]; 

x = t ; 

//differentiation of x can be expressed as 
//summation of xl and x2 

subplot (3,1,1) 
a = gca () ; 

a.x_location = ’’origin”; 
a.y_location = ’’origin”; 
plot (t ,xl) 
xt it le ( ’ xl ( t ) ’) 
subplot (3,1,2) 
a = gca () ; 
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a.x_location = ’’origin”; 
a.y_location = ’’origin”; 
plot2d3 ( ’ gnn ’ ,t , x2) 
xt it le ( ’ x2 ( t ) ’) 
subplot (3,1,3) 
a = gca () ; 

a. x.location = ’’origin”; 
a.y_location = ’’origin”; 
plot (t , x) 
xt it le ( ’ x ( t ) ’) 


Scilab code Exa 4.12 clear 

//clear// 

//Example 4 . 1 2 : Continuous Time Fourier Transform: 
//Derivative property 

clear ; 
clc ; 
close ; 

// CTFT 

t = -1 : 0 . 1 : 1 ; 

xl = ones ( 1 , length (t )) ; 

x2 = [-l,zeros(l,length(t)-2),-l]; 

x = t ; 

//differentiation of x can be expressed as 
//summation of xl and x2 

subplot (3,1,1) 
a = gca () ; 

a. x.location = ’’origin”; 
a.y_location = ’’origin”; 
plot (t , xl ) 
xt it le ( ’ xl ( t ) ’) 
subplot (3,1,2) 
a = gca () ; 

a. x.location = ’’origin”; 
a.y_location = ’’origin”; 
plot2d3 ( ’gnn ’ ,t , x2) 
xt it le ( ’ x2 ( t ) ’) 
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subplot (3,1,3) 
a = gca () ; 

a. x.location = ’’origin”; 
a.y_location = ’’origin”; 
plot (t , x) 
xt it le ( ’ x ( t ) ’) 


Scilab code Exa 4.18 clear 

//clear// 

//Example 4 . 1 8 : Frequency Response of Ideal Low pass 
Filter 

// X(jW)= 1, from — T1 to T1 

clear ; 
clc ; 
close ; 

Wc = 10 ; //I rad /sec 

W = -Wc:0.1:Wc; //Passband of filter 
HWO = 1; //Magnitude of Filter 

HW = HW0*ones (1 , length (W) ) ; 

//Inverse Continuous —time Fourier Transform 

t = -%pi : %pi / length (W) : %pi ; 

Dw = 0.1; 

ht = ( 1/ (2* 0 /,pi ) ) *HW * exp ( sqrt ( -1 ) *W ’ * t ) *Dw ; 
ht = real ( ht ) ; 
f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot (W , HW) ; 

xt it le (’ Frequency Response of Filter H(jW)’) 
subplot (2,1,2) 
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a = gca () ; 

a.y_location = ” origin”; 
a. x.location =” origin”; 
plot (t , ht ) ; 

xt it le (' Impulse Response of Filter h(t)’) 


Scilab code Exa 4.18 clear 

//clear// 

//Example 4 . 1 8 : Frequency Response of Ideal Low pass 
Filter 

// X(jW)= 1, from — T1 to T1 

clear ; 
clc ; 
close ; 

Wc = 10 ; //I rad /sec 

W = -Wc:0.1:Wc; //Passband of filter 
HWO = 1; //Magnitude of Filter 

HW = HW0*ones (1 , length (W) ) ; 

//Inverse Continuous— time Fourier Transform 

t = -%pi : %pi / length (W) : %pi ; 

Dw = 0.1; 

ht = ( 1/ (2* °/ 0 pi ) ) *HW * exp ( sqrt ( - 1 ) *W ’ * t ) *Dw ; 
ht = real ( ht ) ; 
f igure 

subplot (2,1)1) 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot (W , HW) ; 

xt it le (’ Frequency Response of Filter H(jW)’) 
subplot (2,1,2) 
a = gca () ; 

a.y_location =” origin”; 
a.x_location =” origin”; 
plot (t , ht ) ; 

xt it le (’ Impulse Response of Filter h(t)’) 
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Scilab code Exa 4.22 clear 

//clear// 

//Figure 4.22 

//Plotting Continuous Time Fourier Transform of 
//Impulse Response h(t)= exp (— A* t ) u ( t ) , t >0 

clear ; 
clc ; 
close ; 

// Analog Signal 
A =1; //Amplitude 

Dt = 0.005; 
t = 0 : Dt : 10 ; 
ht = exp ( - A*t ) ; 

// Continuous —time Fourier Transform 

Wraax = 2* °/„pi * 1 ; //Analog Frequency = 1Hz 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 

HW = ht* exp ( - sqrt ( - 1 ) *t 1 * W) * Dt ; 

HW.Mag = abs (HW) ; 

W = [-mtlb_fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 

HW_Mag = [mt lb_f 1 iplr ( HW_Mag ) , HW_Mag (2 : 1001) ] ; 

//Plotting Continuous Time Signal 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 

plot (t , ht ) ; 

x label (’ t in sec.’); 

title (' Impulse Response h ( t ) ’ ) 

f igure 
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//Plotting Magnitude Response of CTS 

a = gca () ; 

a.y_location = ’’origin”; 
plot (W , HW_Mag ) ; 

xlabel (’ Frequency in Radians / Seconds > W’ ) ; 

title (’ Frequency Response H(jW)’) 

Scilab code Exa 4.22 

clear 

//clear// 

//Figure 4.22 

//Plotting Continuous Time Fourier Transform of 
//Impulse Response h(t)= exp (— A* t ) u ( t ) , t >0 

clear ; 
clc ; 
close ; 

// Analog Signal 
A =1; //Amplitude 

Dt = 0.005; 
t = 0 : Dt : 10 ; 
ht = exp ( - A*t ) ; 

// Continuous —time Fourier Transform 

Wmax = 2* °/„pi * 1 ; //Analog Frequency = 1Hz 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 

HW = ht* exp ( - sqrt ( - 1 ) *t ’ * W) * Dt ; 

HW.Mag = abs (HW) ; 

W = [-mtlb_fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 

HW.Mag = [mtlb.f liplr (HW_Mag) , HW_Mag (2 : 1001) ] ; 

//Plotting Continuous Time Signal 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
plot (t , ht ) ; 
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x label (’ t in sec.’); 

title (’ Impulse Response h(t)’) 

f igure 

//Plotting Magnitude Response of CTS 

a = gca () ; 

a.y_location = ’’origin”; 
plot (W , HW_Mag ) ; 

xlabel (’ Frequency in Radians / Seconds > W’ ) ; 

title (’ Frequency Response H(jW)’) 


Scilab code Exa 4.23 clear 

//Figure 4 . 2 3 : Mult ip 1 icat ion Property of CTFT 
clear ; 
clc ; 
close ; 

W1 = -1 : 0 . 1 : 1 ; 

W2 = -2 : 0 . 1 : 2 ; 

W = -3 : 0 . 1 : 3 ; 

//Fourier Transform of sine funcion is square wave 

XW1 = ( 1/ %pi )* ones ( 1 , length (W1 )) ; //CTFT of xl ( t ) 

XW2 = ( 1 / (2* °/ 0 pi ))* ones ( 1 , length (W2 )); //CTFT of x2 ( t ) 
XW = (1/2) *convol (XW1 , XW2) ; //CTFT of x ( t )=xl ( t ) *x2 ( t 

^ ) 

//X(jw) = linear convolution of Xl(jw)and X2(jw) 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
a.x_location = ’’origin”; 
plot (W , XW) ; 

xlabel (’ Frequency in Radians / Seconds > W’ ) ; 

title (’Multiplication Property X(jW) ’ ) 
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Scilab code Exa 4.23 

clear 

//Figure 4 . 2 3 : Mult ip 1 icat ion Property of CTFT 
clear ; 
clc ; 
close ; 

W1 = -1 : 0 . 1 : 1 ; 

W2 = -2 : 0 . 1 : 2 ; 

W = -3 : 0 . 1 : 3 ; 

//Fourier Transform of sine funcion is square wave 

XW1 = ( 1/ %pi )* ones ( 1 , length (W1 )) ; //CTFT of xl ( t ) 

XW2 = ( 1 / (2* °/ 0 pi ))* ones ( 1 , length (W2 )); //CTFT of x2 ( t ) 
XW = (1/2) *convol (XW1 , XW2) ; //CTFT of x ( t )=xl ( t ) *x2 ( t 

^ ) 

//X(jw) = linear convolution of Xl(jw)and X2(jw) 

f igure 
a = gca () ; 

a.y_location = ’’origin”; 
a.x_location = ’’origin”; 
plot (W , XW) ; 

xlabel (’ Frequency in Radians / Seconds > W’ ) ; 

title (’ Multiplication Property X(jW) ’ ) 
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Chapter 9 

The Discreet Time Fourier 
Transform 
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Chapter 10 

The Discreet Time Fourier 
Transform 


Scilab code Exa 5.1 Discrete Time Fourier Transform of discrete sequence 


//clear// 

//Example 5.1: Discrete Time Fourier Transform of 
discrete sequence 
//x[n]= (a~n).u[n], a>0 and a<0 

clear ; 
clc ; 
close ; 

/ / DTS Signal 

al = 0.5; 
a2 = -0.5; 
raax_limit = 10; 
for n = 0 : raax_l irait - 1 
xl(n+l) = ( al ~n) ; 
x2 (n + 1 ) = ( a2 ~n) ; 
end 

n = 0 : raax_l imit - 1 ; 
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// D iscrete — time Fourier Transform 

Wmax = 2*%pi; 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 
xl = xl J ; 
x2 = x2 5 ; 

XW1 = xl* exp ( - sqrt ( - 1 ) *n J * W) ; 

XW2 = x2* exp (- sqrt (- 1 ) *n J *W) ; 

XW1 _Mag = abs (XW1) ; 

XW2_Mag = abs (XW2) ; 

W = [-mtlb_fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 

XW1 _Mag = [mtlb.f liplr (XWl.Mag) , XWl_Mag (2 : 1001) ] ; 
XW2_Mag = [mtlb_f liplr (XW2_Mag) , XW2_Mag (2 : 1001) ] ; 
[XWl_Phase , db] = phasemag (XW1 ) ; 

[XW2_Phase , db] = phasemag ( XW2 ) ; 

XWl_Phase = [-mt lb_f 1 iplr ( XW1 .Phase ), XW1 _Phase 
( 2 : 1001 )] ; 

XW2_Phase = [-mt lb_f 1 iplr ( XW2_Phase ) , XW2_Phase 
( 2 : 1001 )] ; 

//plot for a>0 

f igure 

subplot (3,1,1) ; 
plot2d3( ’ gnn ’ , n , xl ) ; 

xtitle (’ Discrete Time Sequence x[n] for a>0 ’ ) 
subplot (3,1,2); 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d (W , XWl_Mag) ; 

title (’ Magnitude Response abs(X(jW))’) 
subplot (3,1,3); 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d (W , XWl_Phase ) ; 
title (’Phase Response <(X(jW))’) 
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//plot for a<0 

f igure 

subplot (3,1,1) ; 
plot2d3( ’ gnn ’ ,n,x2) ; 

xtitle (’ Discrete Time Sequence x[n] for a>0 ’ ) 
subplot (3,1,2); 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d (W , XW2_Mag ) ; 

title (’ Magnitude Response abs(X(jW))’) 
subplot (3,1,3); 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d (W , XW2_Phase ) ; 
title (’Phase Response <(X(jW))’) 

Scilab code Exa 5.1 

Discrete Time Fourier Transform of discrete sequence 

//clear// 

//Example 5.1: Discrete Time Fourier Transform of 
discrete sequence 
//x[n] = (a'n) ,u[n] , a>0 and a<0 

clear ; 
clc ; 
close ; 

// DTS Signal 

al = 0.5; 
a2 = -0.5; 
raax_limit = 10; 
for n = 0 : raax_l irait - 1 
xl (n+1) = (al ~n) ; 
x2 (n + 1 ) = ( a2 ~n) ; 
end 

n = 0 : raax_l imit - 1 ; 
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// D iscrete — time Fourier Transform 

Wmax = 2*%pi; 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 
xl = xl J ; 
x2 = x2 5 ; 

XW1 = xl* exp ( - sqrt ( - 1 ) *n J * W) ; 

XW2 = x2* exp (- sqrt (- 1 ) *n J *W) ; 

XW1 _Mag = abs (XW1) ; 

XW2_Mag = abs (XW2) ; 

W = [-mtlb_fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 

XW1 _Mag = [mtlb.f liplr (XWl.Mag) , XWl_Mag (2 : 1001) ] ; 
XW2_Mag = [mtlb_f liplr (XW2_Mag) , XW2_Mag (2 : 1001) ] ; 
[XWl_Phase , db] = phasemag (XW1 ) ; 

[XW2_Phase , db] = phasemag ( XW2 ) ; 

XWl_Phase = [-mt lb_f 1 iplr ( XW1 .Phase ), XW1 _Phase 
( 2 : 1001 )] ; 

XW2_Phase = [-mt lb_f 1 iplr ( XW2_Phase ) , XW2_Phase 
( 2 : 1001 )] ; 

//plot for a>0 

f igure 

subplot (3,1,1) ; 
plot2d3( ’ gnn ’ , n , xl ) ; 

xtitle (’ Discrete Time Sequence x[n] for a>0 ’ ) 
subplot (3,1,2); 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d (W , XWl_Mag) ; 

title (’ Magnitude Response abs(X(jW))’) 
subplot (3,1,3); 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d (W , XWl_Phase ) ; 
title (’Phase Response <(X(jW))’) 
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//plot for a<0 

f igure 

subplot (3,1,1) ; 
plot2d3( ’ gnn ’ ,n,x2) ; 

xtitle (’ Discrete Time Sequence x[n] for a>0 ’ ) 
subplot (3,1,2); 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d (W , XW2_Mag ) ; 

title (’ Magnitude Response abs(X(jW))’) 
subplot (3,1,3); 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d (W , XW2_Phase ) ; 
title (’Phase Response <(X(jW))’) 


Scilab code Exa 5.2 Discrete Time Fourier Transform 

//clear// 

//Example 5.2: Discrete Time Fourier Transform of 
//x[n] = (a"abs(n)) a>0 and a<0 

clear ; 
clc ; 
close ; 

// DTS Signal 

a = 0.5; 
raax_limit = 10; 
n = -max_l imit +1 : max_l imit - 1 ; 
x = a~ abs (n) ; 

// D iscrete — time Fourier Transform 

Wmax = 2*%pi; 
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K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 

XW = x* exp (- sqrt ( -1) *n ’ *W) ; 

XW.Mag = real (XW) ; 

W = [-mtlb_fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 

XW.Mag = [mtlb.f liplr (XW_Mag) , XW_Mag (2 : 1001) ] ; 

//plot for abs(a)<l 

f igure 

subplot (2,1,1) ; 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d3( ’ gnn ’ ,n,x) ; 

xtitle(’ Discrete Time Sequence x[n] for a>0 ’ ) 
subplot (2,1,2); 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d (W , XW_Mag ) ; 

title (’Discrete Time Fourier Transform X(exp(jW))’) 

Scilab code Exa 5.2 

Discrete Time Fourier Transform 

//clear// 

//Example 5.2: Discrete Time Fourier Transform of 
//x[n] = (a'abs(n)) a>0 and a<0 

clear ; 
clc ; 
close ; 

/ / DTS Signal 

a = 0.5; 
raax_limit = 10; 
n = -max_l imit +1 : max_l imit - 1 ; 
x = a~ abs (n) ; 


140 


12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

1 

2 

3 

4 

5 

6 


// D iscrete — time Fourier Transform 

Wmax = 2*%pi; 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 

XW = x* exp (- sqrt ( -1) *n ’ *W) ; 

XW.Mag = real (XW) ; 

W = [-mtlb_fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 

XW.Mag = [mtlb_f liplr (XW_Mag) , XW_Mag (2 : 1001) ] ; 

//plot for abs(a)<l 

f igure 

subplot (2,1,1) ; 
a = gca () ; 

a.y_location =” origin”; 
a. x_location =” origin”; 
plot2d3( ’ gnn ’ , n , x) ; 

xtitle (’ Discrete Time Sequence x[n] for a>0 ’ ) 
subplot (2,1,2); 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d (W , XW.Mag) ; 

title (’ Discrete Time Fourier Transform X(exp(jW))’) 


Scilab code Exa 5.3 Discrete Time Fourier Transform 

//clear// 

//Example 5.3: Discrete Time Fourier Transform of 
//x[n] = 1 , abs(n)<=Nl 

clear ; 
clc ; 
close ; 
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// DTS Signal 

N 1 = 2; 
n = -N1 : N 1 ; 

x = ones ( 1 , length (n) ) ; 

// D iscrete — time Fourier Transform 

Wmax = 2*%pi; 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 

XW = x* exp ( - sqrt ( - 1 ) *n ’ * W) ; 

XW.Mag = real (XW) ; 

W = [-mtlb_fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 

XW_Mag = [mt lb_f 1 iplr ( XW_Mag ) , XW_Mag (2 : 1001) ] ; 

//plot for abs(a)<l 

f igure 

subplot (2,1,1) ; 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d3( ’ gnn ’ ,n,x) ; 

xtitle (’ Discrete Time Sequence x[n] ’) 
subplot (2,1,2); 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d (W , XW.Mag) ; 

title (’Discrete Time Fourier Transform X(exp(jW))’) 

Scilab code Exa 5.3 

Discrete Time Fourier Transform 

//clear// 

//Example 5.3: Discrete Time Fourier Transform of 
//x[n]= 1 , abs(n)<=Nl 

clear ; 
clc ; 
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close ; 

/ / DTS Signal 

N 1 = 2; 
n = -N1 : N 1 ; 

x = ones ( 1 , length (n) ) ; 

// D iscrete — time Fourier Transform 

Wraax = 2*%pi; 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 

XW = x* exp ( - sqrt ( - 1 ) *n ’ * W) ; 

XW.Mag = real (XW) ; 

W = [-mtlb_fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 

XW_Mag = [mt lb_f 1 iplr ( XW_Mag ) , XW_Mag (2 : 1001) ] ; 

//plot for abs(a)<l 

f igure 

subplot (2,1,1) ; 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d3( ’ gnn ’ ,n,x) ; 

xtitle(’ Discrete Time Sequence x[n] ’) 
subplot (2,1,2); 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
plot2d (W , XW_Mag ) ; 

title (’Discrete Time Fourier Transform X(exp(jW))’) 


Scilab code Exa 5.5 Scilab code Exa 5.5 Time Fourier Transform:x[n] = 
cos(nWo) 
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//clear// 

// Example5 . 5 : D is cr et e Time Fourier Transform : x [ n] = 
cos (nWo) 

clear ; 
clc ; 
close ; 

N = 5; 

Wo = 2* °/„pi / N ; 

W = [-Wo ,0 , Wo] ; 

XW =[°/.pi ,0,'/.pi] ; 

// 

f igure 
a = gca () ; 

a.y_location =” origin”; 

a. x.location =” origin”; 

plot2d3 ( ’ gnn ’ , W , XW , 2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

xlabel ( ’ 

W’) ; 

title ( ’DTFT of cos(nWo)’) 
disp (Wo) 

Time Fourier Transform:x[n|= cos(nWo) 

//clear// 

// Example5 . 5 : D is cr et e Time Fourier Transform : x [ n] = 
cos (nWo) 

clear ; 
clc ; 
close ; 

N = 5; 

Wo = 2* °/„pi / N ; 

W = [-Wo ,0 , Wo] ; 

XW =[*/.pi ,0,*/.pi] ; 

it 

f igure 
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a = gca () ; 

a.y_location = ” origin”; 

a. x.location =” origin”; 

plot2d3 ( ’ gnn ’ ,W,XW ,2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

xlabel ( ’ 

W’) ; 

title ( ’DTFT of cos(nWo)’) 
disp (Wo) 


Scilab code Exa 5.6 Discrete Time Fourier Transform 

//clear// 

//Example5 . 6 : Discrete Time Fourier Transform of 

// Periodic Impulse Train 

clear ; 

clc ; 

close ; 

N = 5; 

N 1 = -3*N : 3*N ; 
xn = [zeros(l,N-l),l]; 
x = [1 xn xn xn xn xn xn] ; 

ak = 1/N; 

XW = 2* °/„pi * ak* ones (1 , 2*N) ; 

Wo = 2* °/„pi / N ; 
n = -N : N - 1 ; 

W = Wo *n ; 
f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location =” origin”; 
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a.x_location =” origin”; 

plot2d3( ’ gnn ’ ,N1 ,x,2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

xlabel ( ’ 

n ’ ) ; 

title (’ Periodic Impulse Train’) 
subplot (2,1,2) 
a = gca () ; 

a.y_location =” origin”; 

a. x.location =” origin”; 

plot2d3 ( ’gnn ’ , W , XW , 2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

xlabel ( ’ 

W’) ; 

title (’DTFT of Periodic Impulse Train’) 
disp (Wo) 

Scilab code Exa 5.6 

Discrete Time Fourier Transform 

//clear// 

// Example5 . 6 : D is cr et e Time Fourier Transform of 

/ / Periodic Impulse Train 

clear ; 

clc ; 

close ; 

N = 5; 

N 1 = -3*N : 3*N ; 
xn = [zeros(l,N-l),l]; 
x = [1 xn xn xn xn xn xn] ; 

ak = 1/N; 

XW = 2* °/„pi * ak* ones (1 , 2*N) ; 

Wo = 2* °/„pi /N ; 
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n = -N : N - 1 ; 

W = Wo *n ; 
f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location = ” origin”; 

a. x.location =” origin”; 

plot2d3( ’ gnn ’ ,N1 ,x,2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

xlabel ( ’ 

n ’ ) ; 

t i 1 1 e ( ’ Periodic Impulse Train’) 
subplot (2,1,2) 
a = gca () ; 

a.y_location =” origin”; 

a.x_location =” origin”; 

plot2d3( ’gnn ’ , W , XW , 2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

polyl . thickness = 3; 

xlabel ( ’ 

W’) ; 

title (’DTFT of Periodic Impulse Train’) 
disp (Wo) 


Scilab code Exa 5.7 Frequency Shifting Property of DTFT 

//clear// 

//Example 5 . 7 : Frequency Shifting Property of DTFT: 
Frequency Response of Ideal Low pass Filter and 
HPF 
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clear ; 
clc ; 
close ; 

Wc = 1 ; //I rad /sec 

W = -Wc:0.1:Wc; //Passband of filter 
HO = 1; //Magnitude of Filter 

HlpW = H0*ones (1 , length (W) ) ; 

Whp 1 = W+%pi; 

Whp2 = -W-°/„pi; 
f igure 

subplot (2,1)1) 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
a . data_bounds = [ — %p i , 0 ; %pi ,2] ; 
plot2d (W , HlpW) ; 

xt it le (’ Frequency Response of LPF H(exp(jW))’) 
subplot (2,1,2) 
a = gca () ; 

a.y_location =” origin”; 
a.x_location =” origin”; 
a.data_bounds = [- 2*°/ 0 pi , 0 ; 2* %pi ,2] ; 
plot2d ( Whpl , HlpW) ; 
plot 2d (Whp 2 , HlpW) ; 

xt it le (’ Frequency Response of HPF H(exp(jW))’) 

Scilab code Exa 5.7 

Frequency Shifting Property of DTFT 

//clear// 

//Example 5 . 7 : Frequency Shifting Property of DTFT: 
Frequency Response of Ideal Low pass Filter and 
HPF 
clear ; 
clc ; 
close ; 

Wc = 1 ; //I rad /sec 
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W = -Wc:0.1:Wc; //Passband of filter 
HO = 1; //Magnitude of Filter 

HlpW = H0*ones (1 , length (W) ) ; 

Whp 1 = W + %pi ; 

Whp2 = -W-°/ 0 pi ; 
f igure 

subplot (2,1)1) 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
a . data_bounds = [ — %p i , 0 ; %pi , 2] ; 
plot2d (W , HlpW) ; 

xt it le (’ Frequency Response of LPF H(exp(jW))’) 
subplot (2,1,2) 
a = gca () ; 

a.y_location =” origin”; 
a.x.location =” origin”; 
a.data_bounds = [- 2*%pi , 0 ; 2* %pi ,2] ; 
plot2d ( Whpl , HlpW) ; 
plot2d ( Whp2 , HlpW) ; 

xtitle (’ Frequency Response of HPF H( exp (jW) ) ’ ) 


Scilab code Exa 5.9 Expansion Property of DTFT 

//clear// 

//Example 5.9: Time Expansion Property of DTFT 
clear ; 
close ; 
clc ; 

n = -1:11; 

x = [0 , 1 , 2 , 1 , 2 , 1 , 2 , 1 , 2 , 1 , 2 , 0 , 0] ; 

y = [1 ,1 ,1 ,1 ,1] ; 

y_2_n =zeros ( 1 , 2* length (y) +1 ) ; 
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y_2_n ( 1 : 2 : 2* length (y) ) = y; 
y_2_n = [0 y_2_n 0] ; 

y_2_n_ 1 = [0 , y_2_n ( 1 : $ - 1 ) ] ; 

x_r = y_2_n+2* y_2_n_ 1 ; 
y = [0 , y , zeros (1,7)] ; 
f igure 

subplot (4,1,1) 
plot2d3 ( ’ gnn ’ ,n,y) 
title ( ’y [n] 5 ) 
subplot (4,1,2) 
plot2d3 ( ’gnn ’ ,n , y_2_n) 
title ( ’y (2) [n] ’) 
subplot (4,1,3) 
plot2d3( ’ gnn ’ ,n,y_2_n_l) 
title ( ’ y ( 2 ) [n-1] ’) 
subplot (4,1,4) 
plot2d3( ’gnn’ ,n,x) 

title(’x[n] = y(2) [ n] + 2*y ( 2 ) [n — 1] ’) 

Scilab code Exa 5.9 

Expansion Property of DTFT 

//clear// 

//Example 5.9: Time Expansion Property of DTFT 
clear ; 
close ; 
clc ; 

n = -1:11; 

x = [0 , 1 , 2 , 1 , 2 , 1 , 2 , 1 , 2 , 1 , 2 , 0 , 0] ; 

y = [1 ,1 ,1 ,1 ,1] ; 

y_2_n =zeros ( 1 , 2* length (y) +1 ) ; 

y_2_n ( 1 : 2 : 2* length (y) ) = y; 

y_2_n = [0 y_2_n 0] ; 

y_2_n_ 1 = [0 , y_2_n ( 1 : $ - 1 ) ] ; 

x_r = y_2_n+2* y_2_n_ 1 ; 

y = [0 , y , zeros (1,7)] ; 

f igure 
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subplot (4,1,1) 
plot2d3 ( ’ gnn ’ ,n , y) 
title ( ’y [n] ’ ) 
subplot (4,1,2) 
plot2d3 ( ’gnn ’ ,n , y_2_n) 
title ( ’y (2) [n] ’) 
subplot (4,1,3) 
plot2d3( ’ gnn ’ ,n,y_2_n_l) 
title ( ’y (2) [n-1] ’) 
subplot (4,1,4) 
plot2d3( ’gnn ’ ,n,x) 

title( ’x[n] = y(2) [ n] + 2*y ( 2 ) [ n — 1] ’) 


Scilab code Exa 5.12 IDTFT:Impulse Response of Ideal Low pass Filter 

//clear// 

//Example 5 . 1 2 : IDTFT : Impulse Response of Ideal Low 
pass Filter 
clear ; 
clc ; 
close ; 

Wc = 1 ; //I rad /sec 

W = -Wc:0.1:Wc; //Passband of filter 
HO = 1; //Magnitude of Filter 
HlpW = H0*ones (1 , length (W) ) ; 

//Inverse D iscrete — time Fourier Transform 
t = -2* °/„pi : 2* °/ 0 pi / length ( W) : 2* °/„pi ; 
ht = (1/ (2*%pi ) ) *HlpW * exp ( sqrt ( -1) *W ’ *t ) ; 
ht = real ( ht ) ; 
f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location =” origin”; 
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a.x_location =” origin”; 
a . data_bounds = [ — %p i , 0 ; %pi , 2] ; 
plot2d (W , HlpW ,2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 

xt it le (’ Frequency Response of LPF H(exp(jW))’) 
subplot (2,1,2) 
a = gca () ; 

a.y_location =” origin”; 

a. x.location =” origin”; 

a. data_bounds=[-2*%pi , - 1 ; 2* %pi ,2] ; 

plot2d3 ( ’ gnn ’ , t , ht ) ; 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

xt it le (' Impulse Response of LPF h(t)’) 

Scilab code Exa 5.12 

IDTFT:Impulse Response of Ideal Low pass Filter 

//clear// 

//Example 5 . 1 2 : IDTFT : Impulse Response of Ideal Low 
pass Filter 
clear ; 
clc ; 
close ; 

Wc = 1 ; //I rad /sec 

W = -Wc:0.1:Wc; //Passband of filter 
HO = 1; //Magnitude of Filter 
HlpW = H0*ones (1 , length (W) ) ; 

//Inverse D iscrete — time Fourier Transform 
t = -2* °/„pi : 2* °/„pi / length ( W) : 2* °/ 0 pi ; 
ht = (1/ (2*%pi ) ) *HlpW * exp ( sqrt ( -1) *W ’ *t ) ; 
ht = real ( ht ) ; 
f igure 

subplot (2,1,1) 
a = gca () ; 

a.y_location =” origin”; 
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a.x_location =” origin”; 
a . data_bounds = [ — %p i , 0 ; %pi , 2] ; 
plot2d (W , HlpW ,2) ; 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
polyl . thickness = 3; 

xt it le (’ Frequency Response of LPF H(exp(jW))’) 
subplot (2,1,2) 
a = gca () ; 

a.y_location =” origin”; 

a. x.location =” origin”; 

a. data_bounds=[-2*%pi , - 1 ; 2* %pi ,2] ; 

plot2d3 ( ’ gnn ’ , t , ht ) ; 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

polyl . thickness = 3; 

xt it le (' Impulse Response of LPF h(t)’) 


Scilab code Exa 5.15 Multiplication Property of DTFT 

//clear// 

//Example5 . 1 5 : Multiplication Property of DTFT 
clear ; 
clc ; 
close ; 
n = 1:100; 

x2 = [3/4 , sin (0 . 75* °/ 0 pi *n) . / ( °/ 0 pi *n) ] ; 
xl = [1/2 , sin (0 . 5* °/„pi *n) . / ( %pi *n) ] ; 

x = xl . *x2 ; 

Wraax = °/„pi ; 

K = 1; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 
n = 0:100; 

XW1 = xl* exp ( - sqrt ( - 1 ) *n J * W) ; 
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XW2 = x2* exp ( - sqrt ( - 1 ) *n ’ * W) ; 

XW = x* exp (- sqrt (- 1 ) *n ’ *W) ; 

XW1 _Mag = real (XW1) ; 

XW2_Mag = real (XW2) ; 

XW.Mag = real (XW) ; 

W = [-mtlb_f liplr (W) , W ( 2 : $ ) ] ; // Omega from — Wmax 
to Wmax 

XW1 _Mag = [mtlb_f liplr (XWl.Mag) , XWl_Mag (2 : $) ] ; 
XW2_Mag = [mtlb_f liplr (XW2_Mag) , XW2_Mag (2 : $) ] ; 
XW.Mag = [mtlb.f liplr (XW_Mag) , XW_Mag (2 : $ ) ] ; 
f igure 

subplot (3,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

plot(W,XWl_Mag) ; 

title ( ’DTFT XI ( exp (jW) ) ’) ; 

subplot (3,1,2) 

a = gca () ; 

a.y_location = ’’origin”; 

a.x_location = ’’origin”; 

plot (W , XW2_Mag ) ; 

title ( ’DTFT X2 ( exp (jW) ) ’) ; 

subplot (3,1,3) 

a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
plot (W , XW_Mag ) ; 

title (’ Multiplication Property of DTFT’); 

Scilab code Exa 5.15 

Multiplication Property of DTFT 

//clear// 

//Example5 . 1 5 : Multiplication Property of DTFT 
clear ; 
clc ; 
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close ; 
n = 1:100; 

x2 = [3/4 , sin (0 . 75* %pi *n) . / ( °/ 0 pi *n) ] ; 
xl = [1/2 , sin (0 . 5*‘/.pi*n) . / (°/.pi*n) ] ; 

x = xl . *x2 ; 

Wmax = °/„pi ; 

K = 1; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 
n = 0:100; 

XW1 = xl* exp ( - sqrt ( - 1 ) *n J * W) ; 

XW2 = x2* exp (- sqrt (- 1 ) *n J *W) ; 

XW = x* exp (- sqrt (- 1 ) *n ’ *W) ; 

XW1 _Mag = real (XW1) ; 

XW2_Mag = real (XW2) ; 

XW.Mag = real (XW) ; 

W = [-mtlb_f liplr (W) , W ( 2 : $ ) ] ; // Omega from —Wmax 

to Wmax 

XW1 _Mag = [mtlb_f liplr (XWl.Mag) , XWl_Mag (2 : $) ] ; 
XW2_Mag = [mtlb_f liplr (XW2_Mag) , XW2_Mag (2 : $) ] ; 
XW.Mag = [mtlb.f liplr (XW_Mag) , XW_Mag (2 : $ ) ] ; 
f igure 

subplot (3,1,1) 
a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

plot(W,XWl_Mag) ; 

title ( ’DTFT Xl ( exp (jW) ) ’) ; 

subplot (3,1,2) 

a = gca () ; 

a.y_location = ’’origin”; 

a. x.location = ’’origin”; 

plot(W,XW2_Mag) ; 

title ( ’DTFT X2 ( exp (jW) ) ’) ; 

subplot (3,1,3) 

a = gca () ; 

a.y_location = ’’origin”; 
a. x.location = ’’origin”; 
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42 plot (W , XW_Mag ) ; 

43 title (’ Multiplication Property of DTFT ’ ) ; 
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Chapter 12 

Time and Frequency 
Characterization of Signals and 
Systems 


Scilab code Exa 6.1 Phase Response and Group Delay 

//clear// 

// Example6 . 1 : Phase Response and Group Delay 

clear ; 
clc ; 
close ; 
fl = 50; 
f 2 = 150; 
f 3 = 300; 
wl = 315; 
tuo 1 = 0.066; 
w2 = 943; 
tuo2 = 0.033; 
w3 = 1888; 
tuo3 = 0.058; 
f = 0:0.1:400; 
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W = 2* °/ 0 pi *f ; 

for i =l:length(f) 

numl(i) = (l+(sqrt (-1) *f (i)/f 1) ~2-2*sqrt (-1) *tuol*( 

f (i)/fl)) ; 

denl(i) = ( 1+ ( sqrt ( -1) *f ( i ) /f 1 ) “2 + 2* sqrt ( -1) * tuo 1 * ( 
f (i)/fl)) ; 

HlW(i) = numl ( i ) /deni ( i ) ; 

num2 ( i ) = ( 1+ ( sqrt ( -1) *f ( i ) /f 2 ) “2 -2* sqrt ( -1) *tuo2 * ( 
f (i)/f2)) ; 

den2 ( i ) = ( 1+ ( sqrt ( -1) *f ( i ) /f 2 ) “2 + 2* sqrt ( -1) * tuo2 * ( 
f (i)/f2)) ; 

H2W(i) = num2 ( i ) /den2 ( i ) ; 

num3 ( i ) = ( 1+ ( sqrt ( -1) *f ( i ) /f 3) “2 -2* sqrt ( -1) * tuo3 * ( 

f (i)/f3)) ; 

den3 ( i ) = ( 1+ ( sqrt ( - 1 ) *f ( i ) / f 3 ) “ 2 + 2* sqrt ( -1 ) * tuo3 * ( 

f (i)/f3)) ; 

H3W(i) = num3 ( i ) /den3 ( i ) ; 

H_W (i) = H 1 W ( i ) * H2W ( i ) ; 

HW ( i ) = H_W ( i ) *H3W ( i ) ; 

phasel(i) = -2* at an ( (2* tuo 1 * ( f ( i ) /f 1 ) ) / ( 1 . 00 1 - ( f ( i 
)/f 1) “2) ) ; 

phase2 ( i ) = -2* at an ( (2* tuo2 * ( f ( i ) /f 2 ) ) / ( 1 . 00 1 - ( f ( i 
) / f 2 ) “2) ) ; 

phase3(i) = -2* atan ( (2* tuo3 * (f ( i ) /f 3) ) / ( 1 . 001 - ( f ( i 
) / f 3 ) “2) ) ; 

phase_total ( i ) = phase 1 ( i ) +phase2 ( i ) +phase3 ( i ) ; 
if (f (i) <=50 ) 

W_phasel(i) = -2* atan ( (2* tuo 1 * ( f ( i ) /f 1 ) ) 

/ ( 1 . 001 - (f (i)/f 1) “2) ) ; 

W_phase2(i) = -2* atan ( (2* tuo2 * ( f ( i ) /f 2 ) ) 

/ ( 1 . 001 - (f (i)/f2) “2) ) ; 

W_phase3(i) = -2* atan ( (2* tuo3 * ( f ( i ) /f 3 ) ) 

/ ( 1 . 001 - (f (i)/f3)“2)) ; 

gr oup_de lay ( i ) = -phase_total ( i ) *0 . 1/ °/ 0 pi ; // 

d e 1 1 a _f= 0.1 

elseif (f (i) >=50 & f(i)<=150) 

W_phasel(i)= -2*7„pi-2*atan((2*tuol*(f(i)/fl)) 

/ ( 1 . 001 - (f (i)/f 1) “2) ) ; 
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W_phase2(i)= -2*atan((2*tuo2*(f(i)/f2))/(1.001-( 
f (i)/f2)~2)) ; 

41 W_phase3(i)= -2* atan ( (2* tuo3 * ( f ( i ) /f 3 ) ) / ( 1 . 001 - ( 

f (i)/f3) ~ 2 ) ) ; 

42 gr oup_de lay ( i ) = -phase_total ( i ) *0 . 1/ (2* °/„pi ) ; 

43 elseif (f (i) > = 150 & f(i)< = 300) 

44 W_phasel(i)= -2* atan ( (2* tuo 1 * ( f ( i ) /f 1 ) ) / ( 1 . 001 - ( 

f (i)/fl)~2)) ; 

45 W_phase2(i)= -4*°/„pi-2*atan((2*tuo2*(f(i)/f2)) 

/ ( 1 . 001 - (f (i)/f2) ~ 2 ) ) ; 

46 W_phase3(i)= -2* atan ( (2* tuo3 * ( f ( i ) /f 3 ) ) / ( 1 . 001 - ( 

f (i)/f3) ~ 2 ) ) ; 

47 gr oup_de lay ( i ) = -phase_total ( i ) *0 . 1/ (4* °/„pi ) ; 

48 elseif (f (i) >300 & f(i)< = 400) 

49 W_phasel(i)= -2* atan ( (2* tuo 1 * ( f ( i ) /f 1 ) ) / ( 1 . 001 - ( 

f (i)/fl)~2)) ; 

50 W_phase2(i)= -2* atan ( (2* tuo2 * ( f ( i ) /f 2 ) ) / ( 1 . 001 - ( 

f (i)/f2) ~ 2 ) ) ; 

51 W_phase3(i)= -6* °/„pi -2* at an ( (2* tuo3 * ( f ( i ) / f 3 ) ) 

/ ( 1 . 001 - (f (i)/f3) ~ 2 ) ) ; 

52 group_delay ( i ) = -phase_t ot al ( i ) *0 . 1 / (4* %pi ) ; 

53 end 

54 if (f (i) ==300 . 1) 

55 W_phase_total ( i ) = 2* °/„pi + W_phase 1 ( i ) + W_phase2 ( i ) + 

W_phase3 ( i ) ; 

56 else 

57 W_phase_total ( i ) = W_phase 1 ( i ) +W_phase2 ( i ) + 

W_phase3 ( i ) ; 

58 end 

59 end 

60 figure 

61 plot2d (f , phase_total , 2) 

62 xtitle ( ’ Principal phase ’ , ’ Frequency (Hz) ’ , ’ Phase (rad) 

63 figure 

64 plot2d (f , W_phase_total , 2) 

65 xt it le (’ unwrapped phase Frequency (Hz) Phase ( rad ) 

’) ; 
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f igure 

plot2d (f , abs (group_delay ) ,2) 

xtitle (’ group delay 7 Frequency (Hz) Group Delay ( 
sec ) ’ ) ; 

Scilab code Exa 6.1 

Phase Response and Group Delay 

//clear// 

// Example6 . 1 : Phase Response and Group Delay 

clear ; 
clc ; 
close ; 
fl = 50; 
f 2 = 150; 
f 3 = 300; 
wl = 315; 
tuo 1 = 0.066; 
w2 = 943; 
tuo2 = 0.033; 
w3 = 1888; 
tuo3 = 0.058; 
f = 0:0.1:400; 

W = 2* °/ 0 pi *f ; 

for i =l:length(f) 

numl(i) = (l+(sqrt (-1) *f (i)/f 1) “2-2*sqrt (-1) *tuol*( 
f (i)/fl)) ; 

denl(i) = ( 1+ ( sqrt ( -1) *f ( i ) /f 1 ) “2 + 2* sqrt ( -1) * tuo 1 * ( 
f (i)/fl)) ; 

HlW(i) = numl ( i ) /deni ( i ) ; 

num2 ( i ) = ( 1+ ( sqrt ( -1) *f ( i ) /f 2 ) “2 -2* sqrt ( -1) * tuo2 * ( 

f (i) /f 2) ) ; 

den2 ( i ) = ( 1+ ( sqrt ( -1) *f ( i ) /f 2 ) “2 + 2* sqrt ( -1) * tuo2 * ( 
f (i)/f2)) ; 

H2W(i) = num2 ( i ) /den2 ( i ) ; 

num3 ( i ) = ( 1+ ( sqrt ( -1) *f ( i ) /f 3) “2 -2* sqrt ( -1) * tuo3 * ( 

f (i) /f 3) ) ; 
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den3 ( i ) = ( 1+ ( sqrt ( -1) *f ( i ) /f 3) “2 + 2* sqrt ( -1) * tuo3 * ( 

f (i)/f3)) ; 

H3W(i) = num3 ( i ) /den3 ( i ) ; 

H_W (i) = H1W (i) *H2W (i) ; 

HW(i) = H_W(i)*H3W(i) ; 

phasel(i) = -2* atan ( (2* tuol * ( f ( i ) /f 1 ) ) / ( 1 . 00 1 - ( f ( i 
)/f 1) “2) ) ; 

phase2 ( i ) = -2* atan ( (2* tuo2 * ( f ( i ) /f 2 ) ) / ( 1 . 00 1 - ( f ( i 
) / f 2 ) “2) ) ; 

phase3(i) = -2* atan ( (2* tuo3 * (f ( i ) /f 3) ) / ( 1 . 001 - (f ( i 
) / f 3 ) “2) ) ; 

phase_total ( i ) = phase 1 ( i ) +phase2 ( i ) +phase3 ( i ) ; 
if (f (i) <=50 ) 

W_phasel(i) = -2* atan ( (2* tuol * (f ( i ) /f 1 ) ) 

/ ( 1 . 001 - (f (i)/fl)“2)) ; 

W_phase2(i) = -2* atan ( (2* tuo2 * ( f ( i ) /f 2 ) ) 

/ ( 1 . 001 - (f (i)/f2) “2) ) ; 

W_phase3(i) = -2* atan ( (2* tuo3 * ( f ( i ) /f 3 ) ) 

/ ( 1 . 001 - (f (i)/f3) “2) ) ; 

group_delay ( i ) = -phase_total (i) *0 . l/°/„pi ; // 

delta_f= 0.1 

elseif (f (i) >=50 & f(i)<=150) 

W_phasel(i)= -2*7„pi-2*atan((2*tuol*(f(i)/fl)) 

/ ( 1 . 001 - (f (i)/fl)“2)) ; 

W_phase2(i)= -2*atan((2*tuo2*(f(i)/f2))/(1.001-( 
f (i)/f2)“2)) ; 

W_phase3(i)= -2*atan((2*tuo3*(f(i)/f3))/(1.001-( 
f (i)/f3)“2)) ; 

group_delay ( i ) = -phase_t ot al ( i ) *0 . 1 / (2* °/„pi ) ; 
elseif (f (i) >=150 & f(i)<=300) 

W_phasel (i)= -2*atan((2*tuol*(f (i)/fl)) / ( 1 . 0 0 1 — ( 
f (i)/fl)“2)) ; 

W_phase2 (i) = -4*%pi -2* atan ( (2*tuo2* (f (i) /f 2) ) 

/ ( 1 . 001 - (f (i)/f2)“2)) ; 

W_phase3(i)= -2*atan((2*tuo3*(f(i)/f3))/(1.001-( 
f (i) /f 3) “2) ) ; 

gr oup_de lay ( i ) = -phase_total ( i ) *0 . 1/ (4* °/„pi ) ; 
elseif (f (i) >300 & f(i)<=400) 
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W_phasel ( i ) = -2*atan((2*tuol*(f (i)/fl)) / ( 1 . 0 0 1 — ( 
f (i)/fl)~2)) ; 

50 W_phase2(i)= -2* atan ( (2* tuo2 * ( f ( i ) /f 2 ) ) / ( 1 . 001 - ( 

f (i)/f2) ~2) ) ; 

51 W_phase3(i)= -6*%pi-2*atan((2*tuo3*(f(i)/f3)) 

/ ( 1 . 001 - (f (i)/f3)~2)) ; 

52 group_delay ( i ) = -phase_t ot al ( i ) *0 . 1 / (4* °/„pi ) ; 

53 end 

54 if (f (i) ==300 . 1) 

55 W_phase_total ( i ) = 2* °/ 0 pi + W_phase 1 ( i ) + W_phase2 ( i ) + 

W_phase3 ( i ) ; 

56 else 

57 W_phase_total ( i ) = W_phase 1 ( i ) +W_phase2 ( i ) + 

W_phase3 ( i ) ; 

58 end 

59 end 

60 figure 

61 plot2d (f , phase_total , 2) 

62 xtitle ( ’Principal phase ’ , ’Frequency(Hz) ’ , ’ Phase ( rad ) 

’) ; 

63 figure 

64 plot2d (f , W_phase_total , 2) 

65 xt it le (’ unwrapped phase Frequency (Hz) Phase ( rad ) 

’) ; 

66 figure 

67 plot2d (f , abs ( group_delay ) , 2) 

68 xtitle (’ group delay Frequency (Hz) Group Delay ( 

sec ) ’ ) ; 


Scilab code Exa 6.3 Analog Lowpass HR filter design 

1 //clear// 

2 // Example6 . 3 : Analog Lowpass HR filter design 
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//Cutoff frequency Fc = 500 Hz 

//Passband ripple 1—0.05 and stopband ripple = 0.05 
clear ; 
close ; 
clc ; 

hs_butt = analpf (5 , ’butt ’ , [0.05,0.05] ,500) ; 
hs_ellip = analpf (5 , ’ellip ’ , [0.05,0.05] ,500) ; 
f r = 0 : . 1 : 2000; 

hf_butt = freq(hs_butt(2) ,hs_butt(3) ,°/ 0 i*fr) ; 
hra_butt = abs ( hf _butt ) ; 

hf_ellip = freq(hs_ellip(2) ,hs_ellip(3) ,°/„i*fr) ; 
hm_ellip = abs (hf _ellip) ; 

//Plotting Magnitude Response of Analog HR Filters 

a = gca () ; 
plot2d(fr ,hm_butt) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 

poly 1 . f oreground = 2; 

polyl . thickness = 2; 

polyl . line_style = 3; 

plot2d(fr , hm_ellip) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 
poly 1 . f oreground = 5; 
polyl . thickness = 2; 
xlabel ( ’Frequency (Hz) ’ ) 

ylabel (’ Magnitude of frequency response’) 
legend([’Butterworth Filter’; ’Elliptic Filter’]) 

Scilab code Exa 6.3 

Analog Lowpass HR filter design 

//clear// 

// Example6 . 3 : Analog Lowpass HR filter design 
//Cutoff frequency Fc = 500 Hz 

//Passband ripple 1—0.05 and stopband ripple = 0.05 
clear ; 
close ; 
clc ; 
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hs_butt = analpf (5 , ’butt ’ , [0.05,0.05] ,500) ; 
hs_ellip = analpf (5 , ’ellip ’ , [0.05,0.05] ,500) ; 
f r = 0 : . 1 : 2000; 

hf_butt = freq(hs_butt(2) ,hs_butt(3) , °/ 0 i*f r) ; 
hm_butt = abs ( hf _butt ) ; 

hf_ellip = freq(hs_ellip(2) ,hs_ellip(3) , °/„ i * f r ) ; 
hra_ellip = abs (hf _ellip) ; 

//Plotting Magnitude Response of Analog HR Filters 

a = gca () ; 
plot2d(fr ,hm_butt) 

polyl = a . chi ldr en ( 1 ) . chi ldren ( 1 ) ; 

poly 1 . f oreground = 2; 

polyl . thickness = 2; 

polyl . line_style = 3; 

plot2d(fr , hm_ellip) 

polyl = a . chi ldr en ( 1 ) . chi ldr en ( 1 ) ; 
poly 1 . f oreground = 5; 
polyl . thickness = 2; 
xlabel ( ’Frequency (Hz) ’ ) 

ylabel (’ Magnitude of frequency response’) 
legend ([’ Butterworth F i 1 1 e r ’ ; ’ E 1 1 i p t i c Filter’]) 


Scilab code Exa 6.4 Bode Plot 

//clear// 

//Example 6 . 4 : Bode Plot 
s = y„s ; 

//Open Loop Transfer Function 

H = syslin( ’c ’ , [20000/ ( s ~ 2 + 1 00* s + 10000 ) ] ) ; //jw 
replaced by s 
elf ; 

bode (H , 0 . 01 , 10000) 
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Scilab code Exa 6.4 

Bode Plot 

//clear// 

//Example 6 . 4 : Bode Plot 
s = °/ 0 s ; 

//Open Loop Transfer Function 

H = syslin( ’c ’ , [20000/ ( s ~ 2 + 1 00* s + 10000 ) ] ) ; //jw 
replaced by s 
elf ; 

bode (H , 0 . 01 , 10000) 


Scilab code Exa 6.5 Bode Plot 

//clear// 

//Example 6 . 5 : Bode Plot 
s = °lo s ; 

//Open Loop Transfer Function 

H = syslinC’c’ , [(100*(l + s))/((10+s)*(100 + s))]) ;//jw 
replaced by s 
elf ; 

bode (H , 0 . 01 , 10000) 


Scilab code Exa 6.5 Bode Plot 

//clear// 

//Example 6 . 5 : Bode Plot 
s = °/o s ; 

//Open Loop Transfer Function 

H = syslin( ’c ’ , [ ( 100* ( 1 + s ) ) / ( ( 10+ s ) * ( 100 + s ) ) ] ) ; //jw 
replaced by s 
elf ; 
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7 bode (H , 0 . 01 , 10000) 
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Chapter 13 
Sampling 
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Chapter 14 
Sampling 


Scilab code Exa 7.1 Sinusoidal signal 

//clear// 

/ / Example 7 .1: Sinusoidal signal 
clear ; 
close ; 
clc ; 

Wm = 2* °/„pi ; 

Ws = 2*Wm; 

t = -2:0.01:2; 

phi = — %p i / 2 ; 

x = cos ( ( Ws /2) *t+phi ) ; 

y = sin ( ( Ws /2) *t ) ; 

subplot (2,1,1) 

a = gca () ; 

a. x.location = ’’origin”; 
a.y_location = ’’origin”; 
plot (t , x) 

title(’cos (Ws/2*t+phi ) ’ ) 

subplot (2,1,2) 
a = gca () ; 
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a.x_location = ’’origin”; 
a.y_location = ’’origin”; 
plot (t , y) 

title(’sin (Ws/2* t ) ’ ) 

Scilab code Exa 7.1 

Sinusoidal signal 

//clear// 

/ / Example 7 .1: Sinusoidal signal 
clear ; 
close ; 
clc ; 

Wm = 2* °/„pi ; 

Ws = 2*Wm; 

t = -2:0.01:2; 

phi = — %p i / 2 ; 

x = cos ( ( Ws /2) *t+phi ) ; 

y = sin ( ( Ws /2) *t ) ; 

subplot (2,1,1) 

a = gca () ; 

a. x.location = ’’origin”; 
a.y_location = ’’origin”; 
plot (t , x) 

title ( ’cos (Ws/2* t+phi ) ’ ) 
subplot (2,1,2) 
a = gca () ; 

a. x.location = ’’origin”; 
a.y_location = ’’origin”; 
plot (t , y) 

title(’sin (Ws/2* t ) ’ ) 


Scilab code Exa 7.2 Digital Differentiator 
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//clear// 

// Example7 . 2 : Digital Differentiator 

syms t n; 

T = 0.1; //Sampling time in seconds 

xct = sin ( °/ 0 pi *t /T) / ( °/ 0 pi *t ) ; 

yet = dif f (xct , t ) ; 
disp (yet , ’ yc ( t )= ' ) ; 
t = n*T ; 

xdn = sin ( %pi *t /T) / ( °/ 0 pi *t ) ; 

ydn = dif f (xdn , n) ; 
disp (ydn , ’ yd [ n] = ’ ) ; 
hdn = T*ydn ; 
disp (hdn , ’hd [n]= ’ ) ; 

/ / Result 

// yc ( t ) = (10*cos(31.415927*t)/t) — (0.3183099* sin 
(31 .415927* t ) /(t ~ 2 ) ) 

//yd [n] = (10* cos (3. 141592 7* n) /n) -3.183* sin (3.1415927* 
n) / ( n ~ 2) 

//hd [n] = ( cos (3. 141 5927* n) /n) -0.3183* sin (3. 14 15 927* n) 
/ ( n 2 ) 

Scilab code Exa 7.2 

Digital Differentiator 

//clear// 

// Example7 . 2 : Digital Differentiator 

syms t n; 

T = 0.1; //Sampling time in seconds 

xct = sin ( °/„pi *t /T) / ( °/ 0 pi *t ) ; 
yet = dif f (xct , t ) ; 
disp (yet , ’ yc ( t )= ’ ) ; 
t = n*T ; 

xdn = sin ( %pi *t /T) / ( °/ 0 pi *t ) ; 

ydn = diff(xdn,n); 
disp (ydn , ’ yd [ n]= ’ ) ; 
hdn = T*ydn ; 
disp (hdn , ’hd [n]= ’ ) ; 
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//Result 

// yc ( t ) = (10*cos (31.415927*t) / 1 ) — (0.3183099* sin 
(31 .415927* t ) /(t A 2 ) ) 

//yd [n] = (10*cos (3. 141 592 7* n) /n) -3.183* sin (3.1415927* 
n) /(n'2) 

/ / lid [n] = ( cos (3. 1415927*n) /n) -0.3183* sin (3. 14 15 927* n) 

/ ( n 2 ) 


Scilab code Exa 7.3 Half Sample Delay system 

//clear// 

// Example7 . 3 : Half Sample Delay system 

syms t n T ; 

//T = 0.1; //Sampling time in seconds 

xct = sin ( %pi *t /T) / ( °/ 0 pi *t ) ; 
t = t -T/2 ; 

yct_del = sin(°/„pi*t/T)/(°/ 0 pi*t); 

disp (yct_del , ’Output of Half Sample delay system 
continuous = ' ) ; 
t = n*T-T/2; 

xdn = sin ( °/„pi *t /T) / ( °/ 0 pi *t ) ; 
ydn_del = xdn; 

disp (ydn_del , ’Output of Half Sample delay system 
discrete ='); 

hdn = T*ydn_del ; 

disp(hdn, ’Impulse Response of discrete time half 
sample delay system=’); 

//Result 

//Output of Half Sample delay system continuous = 
//sin ( 3 . 1 4 * ( t — T / 2 ) /T) / (3. 14* ( t-T/2) ) 

//Output of Half Sample delay system discrete = 

// si n ( 3 . 1 4 * ( n*T— T / 2 ) /T) / ( 3 . 1 4 * ( n*T— T /2 ) ) 

// Impulse Response of discrete time half sample 
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delay system= 

/ / T* sin (3 . 14* ( n* T— T / 2 ) /T ) / ( 3 . 1 4 * ( n *T— ' T / 2 ) ) 

Scilab code Exa 7.3 

Half Sample Delay system 

//clear// 

// Example7 . 3 : Half Sample Delay system 

syms t n T ; 

//T = 0.1; //Sampling time in seconds 

xct = sin ( °/ 0 pi *t /T) / ( °/ 0 pi *t ) ; 
t = t - 1/2 ; 

yct_del = s in ( %pi *t /T) / ( °/ 0 pi *t ) ; 

disp (yct_del , ’Output of Half Sample delay system 
continuous = ' ) ; 
t = n*T-T/2; 

xdn = sin ( %pi *t /T) / ( °/ 0 pi *t ) ; 

ydn_del = xdn; 

disp (ydn_del , ’Output of Half Sample delay system 
discrete =’); 

hdn = T*ydn_del ; 

disp(hdn, ’Impulse Response of discrete time half 
sample delay system=’); 

/ / Result 

//Output of Half Sample delay system continuous = 
//sin ( 3 . 1 4 * ( t — T / 2 ) /T) / (3. 14* ( t-T/2) ) 

//Output of Half Sample delay system discrete = 

// si n ( 3 . 1 4 * ( n*T— ' T / 2 ) /T) / ( 3 . 1 4 * ( n*T— ' T /2 ) ) 

// Impulse Response of discrete time half sample 
delay system= 

// T*sin ( 3 . 1 4 * ( n* T— T / 2 ) /T ) / ( 3 . 1 4 * ( n *T— T / 2 ) ) 


Scilab code Exa 7.4 Period of the sampled signal and Sampling frequency 
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//clear// 

// Example7 . 4 : Finding the period of the sampled 
signal 

//and Sampling frequency 

clear ; 
close ; 
clc ; 

Wm = 2*%pi/9; 

N = floor(2*°/ 0 pi/(2*Wm)) 

disp(N, ’Period of the discrete signal ’) 

Ws = 2* °/„pi / N ; 

disp(Ws, ’The Sampling frequency corresponding to the 

period N ’ ) 

Scilab code Exa 7.4 

Period of the sampled signal and Sampling frequency 

//clear// 

// Example7 . 4 : Finding the period of the sampled 
signal 

//and Sampling frequency 

clear ; 
close ; 
clc ; 

Wm = 2*°/„pi/9; 

N = floor(2*°/ 0 pi/(2*Wm)) 

disp(N, ’Period of the discrete signal’) 

Ws = 2* °/„pi / N ; 

disp(Ws, ’The Sampling frequency corresponding to the 
period N ’ ) 


Scilab code Exa 7.5 Multirate Signal Processing 
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//clear// 

// Example7 . 5 : Multirate Signal Processing : Sampling 
Rate Conversion 
//( 1 ) Downsampling by 4 
//( 2 ) Upsampling by 2 

// (3) Upsampling by 2 and followed by downsampling by 

9 

clear ; 
close ; 
clc ; 

Wm = 2*°/„pi/9; // Maximum frequency of signal 
Ws = 2*Wm; //Sampling frequency 
N = floor (2*°/ 0 pi/Ws); //period of discrete signal 
//Original discrete time signal generation and 
Magnitude response 
n = 0 : 0 . 01 : N ; 
x = sin(Wm*n); 

Wraax = 2*%pi/9; 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 

XW = x* exp (- sqrt ( -1) *n ’ *W) ; 

XW.Mag = real (XW) ; 

XW_Mag = XW_Mag/max (XW_Mag) ; 

W = [-mtlb.fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 

XW.Mag = [mtlb.fliplr (XW.Mag) , XW_Mag (2 : 1001) ] ; 

//( 1 ) downsampling by 4 and corresponding magnitude 
response 

nl = 0:0.01: N/4; 
y = x ( 1 : 4 : length (x) ) ; 
kl = 0: (K/2000) : K; 

W1 = kl*4*Wmax/K; 

XW4 = y* exp (- sqrt (- 1 ) *nl J *W1 ) ; 

XW4_Mag = real (XW4) ; 

XW4_Mag = XW4_Mag/max ( XW4_Mag ) ; 

W1 = [-mtlb.fliplr(Wl), W 1 ( 2 : $ ) ] ; // Omega from - 

Wmax t o Wmax 
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XW4_Mag = [mtlb_f liplr (XW4_Mag) , XW4_Mag (2 : $) ] ; 

//( 2 ) Upsampling by 2 and corresponding magnitude 
response 

n2 = 0 : 0 . 01 : 2*N ; 
z = zeros (1 , length (n2) ) ; 
z ([ 1 : 2 : length ( z )] ) = x; 
k2 = 0 : (K/500) : K ; 

W2 = k2*Wmax/ (2*K) ; 

XW2 = z* exp(-sqrt(-l)*n2 J *W2) ; 

XW2_Mag = real (XW2) ; 

XW2_Mag = XW2_Mag/raax ( XW2_Mag ) ; 

W2 = [-mtlb.f liplr (W2) , W2 ( 2 : $ ) ] ; // Omega from 

Wmax t o Wmax 

XW2_Mag = [mt lb_f 1 iplr ( XW2_Mag ) , XW2_Mag (2 : $) ] ; 
//( 3 ) Upsampling by 2 and Downsampling by 9 
corresponding magnitude response 
n3 = 0:0.01:2*N/9; 
g = z ( [1 : 9 : length (z) ]) ; 
k3 = 0 : K/ (9*500) :K; 

W3 = k3 *9* Wmax / ( 2* K) ; 

XW3 = g* exp ( - sqrt ( - 1 ) *n3 J * W3 ) ; 

XW3_Mag = real (XW3) ; 

XW3_Mag = XW3_Mag/max ( XW3_Mag ) ; 

W3 = [-mtlb_f liplr (W3) , W3(2: $)] ; // Omega from - 
Wmax t o Wmax 

XW3_Mag = [mt lb_f 1 iplr ( XW3_Mag ) , XW3_Mag (2 : $) ] ; 

// 

f igure 

subplot (2,2,1) 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
a.data_bounds = [ — %p i ,0;°/ 0 pi ,1.5] ; 
plot2d(W,XW_Mag ,5) ; 

title (’ Spectrum of Discrete Signal X(exp(jW))’) 
subplot (2,2,2) 
a = gca () ; 

a.y_location =” origin”; 
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a.x_location =” origin”; 
a. data.bounds = [ — %p i ,0;°/ 0 pi ,1.5] ; 
plot2d (W1 , XW4_Mag ,5) ; 

title (’ Spectrum of downsampled signal by 4 X(exp(jW 

/ 4 ))’) 

subplot (2,2,3) 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
a . data_bounds = [ — %p i , 0 ; °/ 0 pi ,1.5]; 
plot 2d (W2 , XW2_Mag , 5) ; 

title (’ Spectrum of Upsampled signal by 2 X(exp(2jW) 
)’) 

subplot (2,2,4) 
a = gca () ; 

a.y_location =” origin”; 
a.x.location =” origin”; 
a . data_bounds = [ — %p i , 0 ; °/ 0 pi ,1.5]; 
plot 2d (W3 , XW3_Mag , 5) ; 

title (’ Spectrum of Upsampled by 2 and Downsampled by 
9 X( exp (2jW /9) ) ’ ) 

Scilab code Exa 7.5 

Multirate Signal Processing 

//clear// 

// Example7 . 5 : Multirate Signal Processing : Sampling 
Rate Conversion 
//( 1 ) Downsampling by 4 
//( 2 ) Upsampling by 2 

//( 3 ) Upsampling by 2 and followed by downsampling by 
9 

clear ; 
close ; 
clc ; 

Wm = 2*°/„pi/9; / / Maximum frequency of signal 
Ws = 2*Wm; //Sampling frequency 
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N = floor(2*°/ 0 pi/Ws);//period of discrete signal 
//Original discrete time signal generation and 
Magnitude response 

n = 0 : 0 . 0 1 : N ; 
x = sin(Wm*n); 

Wmax = 2*%pi/9; 

K = 4; 

k = 0: (K/ 1000) : K ; 

W = k*Wmax/K; 

XW = x* exp ( - sqrt ( - 1 ) *n ’ * W) ; 

XW.Mag = real (XW) ; 

XW_Mag = XW_Mag/max (XW_Mag) ; 

W = [-mtlb_fliplr(W), W(2:1001)]; // Omega from - 
Wmax t o Wmax 

XW.Mag = [mtlb.f liplr (XW_Mag) , XW_Mag (2 : 1001) ] ; 

//( 1 ) downsampling by 4 and corresponding magnitude 
response 

nl = 0:0.01: N/4; 
y = x ( 1 : 4 : length (x) ) ; 
kl = 0 : (K/2000) : K; 

W1 = kl*4*Wmax/K; 

XW4 = y* exp (- sqrt (- 1 ) *nl J *W1 ) ; 

XW4_Mag = real (XW4) ; 

XW4_Mag = XW4_Mag/raax ( XW4_Mag ) ; 

W1 = [-mtlb_fliplr(Wl), W 1 ( 2 : $ ) ] ; // Omega from - 

Wmax t o Wmax 

XW4_Mag = [mt lb_f 1 iplr ( XW4_Mag ) , XW4_Mag (2 : $) ] ; 

//( 2 ) Upsampling by 2 and corresponding magnitude 
response 

n2 = 0 : 0 . 01 : 2*N ; 
z = zeros (1 , length (n2) ) ; 
z ([ 1 : 2 : length ( z )] ) = x; 
k2 = 0: (K/500) :K; 

W2 = k2*Wmax/ (2*K) ; 

XW2 = z* exp (- sqrt (- 1 ) *n2 J *W2 ) ; 

XW2_Mag = real (XW2) ; 

XW2_Mag = XW2_Mag/max ( XW2_Mag ) ; 

W2 = [-mtlb_f liplr (W2) , W2(2:$)]; // Omega from 
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Wmax t o Wmax 

XW2_Mag = [mt lb_f 1 iplr ( XW2_Mag ) , XW2_Mag (2 : $) ] ; 

// (3) Upsampling by 2 and Downsampling by 9 
corresponding magnitude response 
n3 = 0:0.01:2*N/9; 
g = z ( [1 : 9 : length (z) ]) ; 
k3 = 0 : K/ (9*500) :K; 

W3 = k3 *9* Wmax / ( 2* K) ; 

XW3 = g* exp ( - sqrt ( - 1 ) *n3 J * W3 ) ; 

XW3_Mag = real (XW3) ; 

XW3_Mag = XW3_Mag/max (XW3_Mag) ; 

W3 = [-mtlb_f liplr (W3) , W3(2: $)] ; // Omega from 
Wmax t o Wmax 

XW3_Mag = [mt lb_f 1 iplr ( XW3_Mag ) , XW3_Mag (2 : $) ] ; 

// 

f igure 

subplot (2,2,1) 
a = gca () ; 

a.y_location =” origin”; 
a. x_location =” origin”; 
a . data_bounds = [-°/ 0 pi , 0 ; °/ 0 pi ,1.5]; 
plot2d(W,XW_Mag ,5) ; 

title (’ Spectrum of Discrete Signal X(exp(jW))’) 
subplot (2,2,2) 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
a . data_bounds = [-°/ 0 pi , 0 ; °/ 0 pi ,1.5]; 
plot 2d (W1 , XW4_Mag , 5) ; 

title (’ Spectrum of downsampled signal by 4 X(exp(jW 

/4))’) 

subplot (2,2,3) 
a = gca () ; 

a.y_location =” origin”; 
a. x.location =” origin”; 
a.data_bounds = [ — %p i , 0 ; %p i ,1.5] ; 
plot 2d (W2 , XW2_Mag , 5) ; 

title (’ Spectrum of Upsampled signal by 2 X(exp(2jW) 
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) ’) 

78 subplot (2,2,4) 

79 a = gca () ; 

80 a.y_location =” origin”; 

81 a. x.location = ” origin”; 

82 a. data_bounds =[ — %pi ,0;%pi ,1.5] ; 

83 plot2d(W3,XW3_Mag ,5) ; 

84 title (’ Spectrum of Upsampled by 2 and Downsampled by 

9 X( exp (2jW / 9) ) ’) 
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Chapter 9 

The Laplace Transform 


Scilab code Exa 9.1 Lapalce Transform x(t) 

//clear// 

// Example9 . 1 : Lapalce Transform x(t) = exp (— at ) . n ( t ) 

syms t s ; 
a = 3 ; 

y = laplace ( ’%e~ ( — a* t ) ’ , t , s) ; 
disp (y) 

//Result 
// l/( s+a) 

Scilab code Exa 9.1 

Lapalce Transform x(t) 

//clear// 

//Example9 . 1 : Lapalce Transform x(t) = exp (— at ) . u ( t ) 

syms t s ; 
a = 3 ; 

y = laplace ( ’%e~ ( — a* t ) ’ , t , s) ; 
disp (y) 
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//Result 
// l/( s+a) 


Scilab code Exa 9.2 Lapalce Transform x(t) 

//clear// 

// Example9 . 2 : Lapalce Transform x(t) = — exp (— at ) . u(— t 

) 

syms t s ; 
a =3 ; 

y = laplace ( ’%e ~ ( a*— t ) ’ , t , s ) ; 
disp (y) 

/ / Result 
// l/( s+a) 

Scilab code Exa 9.2 

Lapalce Transform x(t) 

//clear// 

// Example9 . 2 : Lapalce Transform x(t) = — exp (— at ) . u(— t 

) 

syms t s ; 
a =3 ; 

y = laplace ( ’%e ~ ( a*— t )’, t , s ) ; 
disp (y) 

/ / Result 
// l/( s+a) 


Scilab code Exa 9.3 Lapalce Transform x(t) 
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//clear// 

// Example9 . 3 : Lapalce Transform x(t) = 3exp( — 2t ) u ( t ) 
— 2exp ( — t ) u ( t ) 

syms t s ; 

y = laplace ( ’ 3*%e" ( — 2 * t ) — 2*%e" ( — t ) ’ , t , s ) ; 
disp (y) 

//Result 

//(3/(s+2))-(2/(s + l)) 

Scilab code Exa 9.3 Lapalce Transform x(t) 

//clear// 

// Example9 . 3 : Lapalce Transform x(t) = 3exp( — 2t ) u ( t ) 
— 2exp ( — t ) u ( t ) 
syms t s ; 

y = laplace (’ 3*%e" ( — 2 * t ) — 2*%e" ( — t )’, t , s ) ; 
disp (y) 

/ / Result 

// (3 / ( s + 2) ) — (2/(s + l)) 


Scilab code Exa 9.4 clear 

//clear// 

// Example9 . 4 : Lapalce Transform x(t) = exp (— 2t ) u ( t ) + 
exp (— t ) (cos3t)u(t) 

syms t s ; 

y = laplace ( ’%e~ ( — 2* t )+%e~ ( — t ) * cos (3 * t ) ’ , t , s ) ; 
disp (y) 

/ / Result 

// [ ( s + 1) /( s A 2 + 2*s +10) ] + [1 / ( s+2) ] refer equation 
9.29 

//Equivalent to (2* s "2 + 5* s +12) / ( ( s ~ 2 + 2* s + 10) * ( s +2) ) 
refer equation 9.30 
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Scilab code Exa 9.4 

clear 

//clear// 

// Example9 . 4 : Lapalce Transform x(t) = exp (— 2t ) u ( t ) + 
exp (— t ) (cos3t)u(t) 

syms t s ; 

y = laplace ( ’%e~ ( — 2* t )+%e~ ( — t ) * cos (3 * t ) ’ , t , s ) ; 
disp (y) 

//Result 

// [ ( s + 1) / ( s ~2 + 2*s +10) ] + [1 / ( s+2) ] refer equation 
9.29 

//Equivalent to (2* s "2 + 5* s +12) / ( ( s ~ 2 + 2* s + 10) * ( s +2) ) 
refer equation 9.30 


Scilab code Exa 9.5 clear 

//clear// 

// Example9 . 5 : Lapalce Transform of x ( t )=s ( t ) — (4/3) exp 
( — t)u(t)+(l/3) exp (2 t ) u ( t ) 

syms t s ; 

y = laplace ( ’ — (4/3) *%e~ ( — t ) + (1/3) *%e ~ (2* t ) ’ , t , s) ; 

y = l+y; 

disp (y) 

/ / Result 

//[-4/(3*(s + l))] + [l/(3*(s— 2))] + l 

Scilab code Exa 9.5 

clear 

//clear// 
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// Example9 . 5 : Lapalce Transform of x ( t )=s ( t ) — (4/3) exp 
( — t)u(t)+(l/3) exp (2t)u(t) 

syms t s ; 

y = laplace ( ’ — (4/3) *%e~ ( — t ) + (1/3) *%e ~ (2* t ) ’ , t , s) ; 

y = l+y ; 

disp (y) 

//Result 

/ / [ 4 / (3*(s + l)) ] + [! / (3*(s— 2))] + l 


Scilab code Exa 9.6 clear 

//clear// 

/ / Example9 . 6 

//Lapalce Transform x(t) = exp(— at)u(t) , 0<t<T 

syms t s ; 
a = 3 ; 

T = 10; 

//t = T; 

y = laplace ( ’%e~ ( — a* t )— %e~ ( — a* t ) *%e~ ( — ( s+a ) *T) ’ , t , s ) 

) 

disp (y) 

/ / Result 

// [ 1 / ( s+a) ] - [(exp(( — s— a)*T) ) / (s+a) ] 

Scilab code Exa 9.6 

clear 

//clear// 

/ / Example9 . 6 

//Lapalce Transform x(t) = exp(— at)u(t) , 0<t<T 

syms t s ; 
a = 3 ; 

T = 10; 
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//t = T; 

y = laplace ( ’%e~ ( — a* t )— %e~ ( — a* t ) *%e~ ( — ( s+a ) *T) ’ , t , s ) 

> 

disp (y) 

/ / Result 

// [ 1 / ( s+a) ] - [(exp(( — s— a)*T) ) / (s+a) ] 


Scilab code Exa 9.7 clear 

//clear// 

/ / Example9 . 7 

//Lapalce Transform x(t) = exp ( — b . abs ( t ) ) . u ( t ) , 0<t< 

T 

//x(t) = exp(— bt ) . u ( t )+exp ( bt ) . u( — t ) 

syms t s ; 
b = 3; 

y = laplace ( ’ %e~ ( — b* t )— %e ~( b* t )’, t , s ) ; 
disp (y) 

//Result 

// [ 1 / ( s +b ) ] — [ 1 / ( s — b ) ] 

Scilab code Exa 9.7 

clear 

//clear// 

/ / Example9 . 7 

//Lapalce Transform x(t) = exp ( — b . abs ( t ) ) . u ( t ) , 0<t< 

T 

//x(t) = exp( — bt ) . u ( t )+exp ( bt ) . u( — t ) 

syms t s ; 
b = 3; 

y = laplace ( ’%e~ ( — b* t )— %e " ( b* t )’, t , s ) ; 
disp (y) 
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//Result 

// [l/(s+b)]-[l/(s-b)] 


Scilab code Exa 9.8 clear 

//clear// 

// Example9 . 8 : Inverse Lapalce Transform 

//X(S) = 1/ ( ( s + 1) ( s+2) ) 

s =°/ 0 s ; 

G =syslin( ’c ’ , ( 1/ ( ( s + 1) * ( s + 2) ) ) ) ; 

disp(G,”G( s )=”) 
plzr (G) 
x=denom(G) ; 

d i s p ( x , ” Ch a r a c t e r i s t i c s Polynomial—’ ) 
y = roots (x) ; 

di sp (y , ” P oles of a system=” ) 

//Result 
// —1 and —2 

Scilab code Exa 9.8 

clear 

//clear// 

// Example9 . 8 : Inverse Lapalce Transform 
//X(S) = 1/ ( ( s + 1) ( s+2) ) 

s =°/oS ; 

G = syslin ( ’ c ’ , (1/ ( ( s + 1) * ( s + 2) ) ) ) ; 

disp(G,”G( s )=”) 
plzr (G) 
x=denom(G) ; 

d i s p ( x , ” Ch a r a c t e r i s t i c s Polynomial^” ) 
y = roots (x) ; 

di sp (y , ” P oles of a system—’ ) 
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/ / Result 
// —1 and —2 


Scilab code Exa 9.9 clear 


Scilab code Exa 9J9 //clear// 

// Example9 . 9 : Inverse Lapalce Transform 

//X(S) = 1/ ( ( s + 1) ( s+2) ) 

s =%s ; 
syms t ; 

[A]=pfss(l/((s+l)*(s+2))) //partial fraction of F ( s ) 
FI = i laplace ( A ( 1 ) , s , t ) 

F2 = ilaplace (A (2) , s , t) 

F=F1 +F2 ; 
disp(F,”f(t)=”) 

/ / Result 

// (%e~ — t ) — (%e~ — (2*t ) ) 

clear 

//clear// 

//Example9 . 9 : Inverse Lapalce Transform 
//X(S) = 1/ ( ( s + 1) ( s+2) ) 

s =°/oS ; 

syms t ; 

[A]=pfss(l/((s+l)*(s+2))) //partial fraction of F ( s ) 
FI = i laplace ( A ( 1 ), s , t ) 

F2 = ilaplace (A (2) , s , t) 

F=F1 +F2 ; 

di sp (F , ” f ( t )=” ) 

//Result 

// (%e" — t ) — (%e' — (2*t ) ) 
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Scilab code Exa 9.10 Inverse Lapalce Transform 

//clear// 

// Example9 . 1 0 : Inverse Lapalce Transform 
//X(S) = 1/ ( ( s + 1) ( s+2) ) Re ( s ) < -l,Re(s)< -2 

s =°/oS ; 

syms t ; 

[A]=pfss(l/((s+l)*(s+2))) //partial fraction of F ( s ) 
FI = i laplace ( A ( 1 ) , s , t ) 

F2 = ilaplace (A (2) , s , t) 

F = -F1-F2; 
disp(F,”f(t)=”) 

/ / Result 

// %e~ — (2*t )— %e~ — t 

Scilab code Exa 9.10 

Inverse Lapalce Transform 

//clear// 

//Example9 . 1 0 : Inverse Lapalce Transform 
//X(S) = l/(( s + 1) ( s+2) ) Re ( s ) < -l,Re(s)< -2 

s =%s ; 
syms t ; 

[A]=pfss(l/((s+l)*(s+2))) //partial fraction of F ( s ) 
FI = ilaplace (A ( 1 ), s , t ) 

F2 = ilaplace (A (2) , s , t) 

F = -F1-F2; 
disp(F,”f(t)=”) 

/ / Result 

// %e" -(2*t )— %e~ — t 
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Scilab code Exa 9.11 Inverse Lapalce Transform 

//clear// 

// Example9 . 1 1 : Inverse Lapalce Transform 
//X(S) = l/( ( s + 1) ( s+2) ) — 2< Re ( s ) < -1 

s =%s ; 
syms t ; 

[A]=pfss(l/((s+l)*(s+2))) //partial fraction of F ( s ) 
FI = i laplace ( A ( 1 ) , s , t ) 

F2 = ilaplace (A (2) , s , t) 

F = -F1+F2; 
disp(F,”f(t)=”) 

/ / Result 

// — (%e~ — t ) — (%e~ — (2*t ) ) 

Scilab code Exa 9.11 

Inverse Lapalce Transform 

//clear// 

//Example9 . 1 1 : Inverse Lapalce Transform 
//X(S) = 1/ ( ( s + 1) ( s +2) ) — 2< Re ( s ) < -1 

s =%s ; 
syms t ; 

[A]=pfss(l/((s+l)*(s+2))) //partial fraction of F ( s ) 
FI = ilaplace (A ( 1 ), s , t ) 

F2 = ilaplace (A (2) , s , t) 

F = -F1+F2; 
disp (F , ” f ( t ) — ’ ) 

//Result 

// — (%e~ — t)— (%e~— (2*t)) 
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Scilab code Exa 9.12 Inverse Lapalce Transform 


Scilab code Exa 9. 112 //clear// 

// Example9 . 1 2 : Inverse Lapalce Transform 
//X(S) = l/( s + (1/2) ) Re(s)> -1/2 
s =%s ; 

G = sysl in ( ’ c ’ , ( 1/ ( s + 0 . 5) ) ) 
disp(G,”G( s )=”) 
plzr (G) 


Inverse Lapalce Transform 

//clear// 

// Example9 . 1 2 : Inverse Lapalce Transform 
//X(S) = 1 / ( s + ( 1/2) ) Re(s)> -1/2 

s =°/oS ; 

G = syslin ( ’ c ’ , ( 1/ ( s + 0 . 5) ) ) 
di sp (G , ”G( s )=”) 
plzr (G) 


Scilab code Exa 9.13 Inverse Lapalce Transform 

//clear// 

/ / Example9 . 1 3 

//Inverse Lapalce Transform 
//X1(S) = l/(s+l) Re ( s ) > -1 
//X2(S) = 1/ ( ( s + 1) ( s+2) ) Re ( s ) >— 1 
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s =7,s ; 
syms t ; 

G1 = syslin (’ c ’,( 1/ ( s + 1) )) ; 
disp (Gl ,”G( s )=”) 
f igure 
plzr (Gl) 

G2 = syslin ( ’ c ’ ,(l/((s+l)*(s+2)))); 
disp (G2 , ”G( s )=”) 
f igure 
plzr (G2) 

G3 = syslin (’ c (1/ ( s+1) )-( 1/ (( s+1) *( s+2) ))) ; 
disp (G3 , ”G( s )=”) 
f igure 
plzr ( G3 ) 

Scilab code Exa 9.13 

Inverse Lapalce Transform 

//clear// 

/ / Example9 . 13 

//Inverse Lapalce Transform 
//X1(S) = l/(s+l) Re ( s ) > -1 
//X2(S) = 1/ ( ( s+1) ( s + 2) ) Re ( s ) >— 1 

s =%s ; 
syms t ; 

Gl = syslin (’ c ’,( 1/ ( s + 1 ))) ; 
disp (Gl , ”G( s )=”) 
f igure 
plzr (Gl) 

G2 = syslin ( ’ c ’ ,(l/((s+l)*(s+2)))); 
disp (G2 , ”G( s )=”) 
f igure 
plzr (G2) 

G3 = syslin (’ c (1/ ( s+1) )-( 1/ (( s+1) *( s+2) ))) ; 
disp (G3 , ”G( s )=”) 
f igure 
plzr ( G3 ) 
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Scilab code Exa 9.14 Lapalce Transform 

//clear// 

// Example9 .14: Lapalce Transform 
//x(t) = t.exp(— at), t >0 
//x(t) = ( t ~ 2 ) /2 . exp (— at ) , t >0 

s = %s ; 
syms t ; 
a =10; 

xl = laplace(’t *%e~ ( — 1 0 * t ) ’ , t , s) ; 
disp (xl ) 

x2 = laplace ( ’ ( ( t ~ 2) /2) *%e~ ( — 10* t ) ’ , t , s) ; 
disp (x2) 

/ / Result 

//!/(( s + 10) '2) 

// l/((s+10) -3) 

Scilab code Exa 9.14 

Lapalce Transform 

//clear// 

/ / Example9 .14: Lapalce Transform 
//x(t) = t.exp(— at), t >0 
//x(t) = ( t ~ 2 ) /2 . exp (— at ) , t >0 

s = %s ; 
syms t ; 
a =10; 

xl = laplaceC’t *%e~ ( — 10* t ) ’ , t , s) ; 
disp (xl ) 

x2 = laplace ( ’ ( ( t ~ 2) /2) *%e" ( — 10* t ) ’ , t , s) ; 
disp (x2) 
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//Result 
// l/((s + 10) 2) 
// 1 / ( ( s + 10) ~3) 


Scilab code Exa 9.15 Inverse Lapalce Transform 


Scilab code Exa 9.115 //clear// 

// Example9 . 1 5 : Inverse Lapalce Transform 
//X(S) = (2s~2+5s+5) /(( s + 1) "2) (s+2)) Re(s)>-1 

s =%s ; 
syms t ; 

[A]=pfss ((2*(s~2) + 5* s + 5) / ( ( ( s+1 ) ~2) * ( s+2) ) ) ; // 
partial fraction of F ( s ) 

FI = i laplace ( A ( 1 ) , s , t ) 

F2 = ilaplace (A (2) , s , t) 

//F3 = i 1 ap 1 a c e (A( 3 ) ,s,t) 

F = F1+F2; 
disp(F,”f(t)=”) 

/ / Result 

// (2* t * (%e~ — t ) ) — (%e / ' — t ) +(3*%e~ — (2*t ) ) 

Inverse Lapalce Transform 


//clear// 

// Example9 . 1 5 : Inverse Lapalce Transform 
//X(S) = (2s'2+5s+5) / (( s + 1) "2) (s+2) ) Re(s)>-1 

s =°/oS ; 

syms t ; 

[A]=pfss ((2*(s~2) +5*s+5) / (((s+1) "2) *(s+2) ) ) ; // 
partial fraction of F ( s ) 

FI = i laplace ( A ( 1 ), s , t ) 

F2 = ilaplace (A (2) , s , t) 

//F3 = i 1 ap 1 a c e (A( 3 ) ,s,t) 
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F = F1+F2 ; 
disp(F,”f(t)=”) 

/ / Result 

// (2* t * (%e~ — t ) ) — (%e~ — t ) +(3*%e A — (2*t ) ) 


Scilab code Exa 9.16 Initial Value Theorem of Lapalace Transform 


Scilab code Exa 9.116 //clear// 

// Example9 . 1 6 : Initial Value Theorem of Lapalace 
Transform 
syms s ; 

num =poly([12 5 2], , s , ,’coeff’) 
den =poly([20 14 4 l],’s’,’coeff’) 

X = num/ den 
disp (X , ”X( s )=” ) 

SX = s*X ; 

Initial_Value = limit(SX,s,°/„inf) ; 
disp ( Initial_Value ,”x(0)=”) 

/ / Result 

/ / (2*%inf A 3+5*%inf'2 + 12*%inf) /(%inf A 3+4*%inf A 2 + 14* 
%inf + 20) =2 


Initial Value Theorem of Lapalace Transform 

//clear// 

// Example9 . 1 6 : Initial Value Theorem of Lapalace 
Transform 
syms s ; 

num =poly([12 5 2], , s , ,’coeff’) 
den =poly([20 14 4 l],’s’,’coeff’) 

X = num/ den 
disp (X , ”X( s )=” ) 

SX = s*X ; 
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Initial_Value =limit(SX,s,°/„inf) ; 
disp ( Initial_Value ,”x(0)=”) 

/ / Result 

/ / (2*%inf~3 + 5*%inf~2 + 12*%inf) /(%inf~3+4*%inf~2 + 14* 
%inf + 20) =2 


Scilab code Exa 9.17 Analysis and Characterization of LTI System 

//clear// 

// Example9 . 1 7 : Analysis and Characterization of LTI 
System 

//Lapalce Transform h(t) = exp(— t).u(t) 

syms t s ; 

h = laplace ( ’%e~ ( — t ) ’ , t , s ) ; 
disp (h) 

//Result 

//l/(s+l) 

Scilab code Exa 9.17 

Analysis and Characterization of LTI System 

//clear// 

// Example9 . 1 7 : Analysis and Characterization of LTI 
System 

//Lapalce Transform h(t) = exp( — t).u(t) 

syms t s ; 

h = laplace ( ’%e~ ( — t )’, t , s ) ; 
disp (h) 

/ / Result 
//l/(s+l) 
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Scilab code Exa 9.18 Analysis and Characterization of LTI System 


//clear// 

// Example9 . 1 8 : Analysis and Characterization of LTI 
System 

//Lapalce Transform x(t) = exp( — abs(t)) 

//x(t) = exp(— t ) . u ( t )+exp ( t ) . u(— t ) 

syms t s ; 

y = lap lace ( ’%e~ ( — t )— %e ~ ( t ) ’ , t , s ) ; 
disp (y) 

/ / Result 

// (l/(s + l))-(l/(s-l)) 

Scilab code Exa 9.18 

Analysis and Characterization of LTI System 

//clear// 

// Example9 . 1 8 : Analysis and Characterization of LTI 
System 

//Lapalce Transform x(t) = exp( — abs(t)) 

//x(t) = exp(— t ) . u ( t )+exp ( t ) . u(— t ) 

syms t s ; 

y = lap lace ( ’%e~ ( — t )— %e ~ ( t ) 7 , t , s ) ; 
disp (y) 

//Result 

// ( 1 / ( s + 1) ) - (l/(s — 1) ) 


Scilab code Exa 9.19 Analysis and Characterization of LTI System 

//clear// 


198 


2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 


// Example9 . 1 9 : Analysis and Characterization of LTI 
System 

//Inverse Lapalce Transform 
//X(S) = ( e ~ s ) / ( s + 1) 

syms s t ; 

hi = exp (-1) ;// Inverse Laplace Transform of exp ( s ) 

H2 =l/(s+l) ; 

h2 = ilaplace (H2 , s , t) 

h = hi *h2 ; 

disp (h , ”h ( t )=” ) 

//Result 

// (18089*(%e~ — t) ) /49171 = 0 . 3 6 78 794 ( %e~-t ) 

Scilab code Exa 9.19 

Analysis and Characterization of LTI System 

//clear// 

// Example9 . 1 9 : Analysis and Characterization of LTI 
System 

//Inverse Lapalce Transform 
//X(S) = (e*s)/(s + l) 

syms s t ; 

hi = exp (-1) ;// Inverse Laplace Transform of exp ( s ) 

H2 =l/(s+l) ; 

h2 = ilaplace (H2 , s , t) 

h = hi *h2 ; 

disp (h , ”h( t )=” ) 

/ / Result 

// ( 18089* (%e* — t )) /49171 = 0 . 3 6 78 794 ( %e A -t ) 


Scilab code Exa 9.20 Scilab code Exa 9.20 Inverse Lapalce Transform 
Inverse Lapalce Transform 
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//clear// 

//Example9 . 2 0 : Inverse Lapalce Transform 
//X(S) = ( ( s — 1) / ((s+l)*(s— 2)) ) 

s =%s ; 
syms t ; 

[A] = pf ss (s/ ( (s + 1) * (s-2) ) ) ; 

[B] = pf ss (l/((s + l) *(s-2) ) ) ; 

FI = i laplace ( A ( 1 ), s , t ) 

F2 = ilaplace (A (2) , s , t) 

F3 = i laplace (B ( 1 ), s , t ) 

F4 = ilaplace (B (2) , s , t) 

F = F1+F2-F3-F4; 

di sp (F , ” f ( t )=” ) 

/ / Result 

// f ( t )= 33333329999999* exp (2* t) 

/ 999999 70000000 + 66 6666 4*%e~ — t /9999997 
// i . e . f(t) = 0.3333334* exp(2*t) + 0 . 6666666* %e' ( — t ) 
//Refer equation 9.120. (l/3)=0.3333 and (2/3) = 

0.66 666 


//clear// 

//Example9 . 2 0 : Inverse Lapalce Transform 

//X(S) = ( ( s — 1) / ((s+l)*(s— 2)) ) 

s =°/oS ; 

syms t ; 

[A] = pf ss (s/ ( ( s + 1) * (s-2) ) ) ; 

[B] = pf ss ( 1/ ( ( s + 1) * (s-2) ) ) ; 

FI = i laplace ( A ( 1 ), s , t ) 

F2 = ilaplace (A (2) , s , t) 

F3 = i laplace (B ( 1 ), s , t ) 

F4 = ilaplace (B (2) , s , t) 

F = F1+F2-F3-F4; 
disp(F,”f(t)=”) 

//Result 

// f ( t )= 33333329999999* exp (2* t) 

/ 999999 70000000 + 6666664* %e" — t /9999997 
//i.e. f(t) =0.3333334* exp(2*t) + 0 . 6666666 *%e" ( — t ) 
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//Refer equation 9.120. (l/3)=0.3333 and (2/3) = 

0.66 666 


Scilab code Exa 9.21 Analysis and Characterization of LTI System 

//clear// 

// Example9 . 2 1 : Analysis and Characterization of LTI 
System 

//Lapalce Transform h(t) = exp(2t)u(t), Re ( s ) >2 

syms t s ; 

X = lap lace ( ’ %e ~ (2 * t ) ’ , t , s ) ; 
disp (X) 

//Result 
// 1 / ( s -2) 

Scilab code Exa 9.21 

Analysis and Characterization of LTI System 

//clear// 

// Example9 . 2 1 : Analysis and Characterization of LTI 
System 

//Lapalce Transform h(t) = exp(2t)u(t), Re ( s ) >2 

syms t s ; 

X = lap lace ( ’ %e ~ (2 * t ) ’ , t , s ) ; 
disp (X) 

//Result 

/ / 1/ ( s 2) 


Scilab code Exa 9.25 Finding Transfer function H(S) of LTI system 
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//clear// 

// Example9 . 2 5 : LTI Systems Characterized by Linear 
Constant 

//Coefficient differential Equation 
//Finding Transfer function H( S ) of LTI system 
//x(t) = exp (— 3t ) . u ( t ) 

// y(t) = [ exp ( — t )— exp ( — 2 1 ) ] . u ( t ) 

syms t s ; 

X = laplace ( ’ %e~ ( — 3* t ) ’ , t , s ) ; 

Y = lap lace ( ’ %e~ ( — t )— %e~ ( — 2* t ) ’ , t , s ) ; 

H = Y/X; 

disp (H) 

/ / Result 

/ / (s+3)*(l/(s + l)— l/(s+2)) 

Scilab code Exa 9.25 

Finding Transfer function H(S) of LTI system 

//clear// 

// Example9 . 2 5 : LTI Systems Characterized by Linear 
Constant 

//Coefficient differential Equation 
//Finding Transfer function H(S) of LTI system 
//x(t) = exp (— 3t ) . u ( t ) 

//y(t) = [ exp ( — t )— exp ( — 2 1 ) ] . u ( t ) 

syms t s ; 

X = laplace ( ’%e~ (— 3* t ) ’ , t , s) ; 

Y = lap lace ( ’ %e~ ( — t )— %e~ ( — 2* t ) ’ , t , s ) ; 

H = Y/X; 

disp (H) 

//Result 

/ / (s+3)*(l/(s + l)— l/(s+2)) 
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Scilab code Exa 9.31 Partial Fraction 


//clear// 

// Example9 . 3 1 : Causal LTI Systems described by 
differential equations 
//and Rational System functions 
//Partial Fraction 
//H(S) = ( ( s — 1) / ((s+l)*(s— 2)) ) 
s =°/oS ; 

syms t ; 

[A] = pfss ( (2*s ~2+4*s-6) /(s~2+3*s+2) ) ; 
di sp ( A , ”H( S )=” ) 

//Result H( S ) = 

lilt ~ 8 

// 1 + s 

// 6 

// 2 + s 

// 

// 2 

Scilab code Exa 9.31 

Partial Fraction 

//clear// 

// Example9 . 3 1 : Causal LTI Systems described by 
differential equations 
//and Rational System functions 
//Partial Fraction 
//H(S) = ( ( s — 1) / ((s+l)*(s— 2)) ) 
s =°/oS ; 

syms t ; 

[A] = pf ss ( (2*s ~2 + 4*s-6) /(s ~2 + 3*s + 2) ) ; 
di sp ( A , ”H( S )=” ) 

//Result H( S ) = 

II tt ~ 8 
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// 

// 1 + s 

// 6 

// 

// 2 + s 

it 

// 2 


Scilab code Exa 9.33 Unilateral Laplace Transform 

//clear// 

// Example9 . 3 3 : U n i 1 at er al Laplace Transform : Time 
Shifting Property 
//x(t) = exp(— a(t+l) ) . u( t + 1) 

syms t s ; 
a = 2 ; 

X = lap lace ( ’ %e~ ( — a * ( t +1) ) ’ , t , s ) ; 
disp (X) 

//Result 
//%e~ — a/ ( s+a) 

Scilab code Exa 9.33 

Unilateral Laplace Transform 

//clear// 

// Example9 . 3 3 : U n i 1 at er al Laplace Transform : Time 
Shifting Property 
//x(t) = exp(— a(t + l)) .u(t + l) 

syms t s ; 
a = 2 ; 

X = laplace ( ’%e~ ( — a * ( t +1) ) ’ , t , s ) ; 
disp (X) 

//Result 
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/ / %e'~ — a / ( s+a) 


Scilab code Exa 9.34 Unilateral Laplace Transform 

//clear// 

// Example9 . 3 4 : U n i 1 at er al Laplace Transform 
//x(t) = s ( t ) T2u ( t ) Te t . u ( t ) 

syms t s ; 
a = 2 ; 

X = laplace ( 7 2+%e~ ( t ) 7 , t , s) ; 

Y = 1+X ; 
disp (X) 
di sp ( Y) 

//Result 

// (2/ s ) +(1/ (s —1) )+l 

Scilab code Exa 9.34 

Unilateral Laplace Transform 

//clear// 

//Example9 . 34 : Unilateral Laplace Transform 
/ / x ( t ) = s ( t ) +2u ( t ) +e ~ t • u ( t ) 

syms t s ; 
a = 2 ; 

X = laplace ( 7 2+%e~ ( t ) 7 , t , s) ; 

Y = 1+X; 
disp (X) 
di sp ( Y) 

/ / Result 

// (2/s)+(l/(s— 1))+1 
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Scilab code Exa 9.35 clear 


//clear// 

// Example9 . 3 5 : U n i 1 at er al Inverse Laplace Transform 
//X(S) = 1/ ( ( s + 1) ( s+2) ) 

s = °/ 0 s ; 
syms t ; 

X = 1/ ( ( s + 1) * ( s + 2) ) ; 
x = ilaplace (X , s , t) ; 
disp (X) 
disp (x) 

/ / Result 

// (%e~ — t ) — (%e~ — (2*t ) ) 

Scilab code Exa 9.35 

clear 

//clear// 

// Example9 . 3 5 : U n i 1 at er al Inverse Laplace Transform 
//X(S) = l/(( s + 1) ( s+2) ) 

s = °/ 0 s ; 
syms t ; 

X = 1/ ( ( s + 1) * ( s + 2) ) ; 
x = ilaplace (X , s , t) ; 
disp (X) 
disp (x) 

//Result 

// (%e~-t)-(%e'-(2*t)) 


Scilab code Exa 9.36 clear 
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//clear// 

// Example9 . 3 6 : U n i 1 at er al Laplace Transform 
//X(S) = ( ( s ~ 2 ) —3) / ( s +2) 

s = °/ 0 s ; 
syms t ; 

[X] = pf ss ( ( ( s ~2) -3) / (s+2) ) ; 
disp (X) 

Scilab code Exa 9.36 

clear 

//clear// 

// Example9 . 3 6 : Un i lat er al Laplace Transform 
//X(S) = ( ( s " 2 ) —3) / ( s +2) 

s = 7o s ; 
syms t ; 

[X] = pf ss ( ( ( s ~2) -3) / (s + 2) ) ; 
disp (X) 


Scilab code Exa 9.37 clear 

//clear// 

// Example9 . 3 7 : U n i 1 at er al Laplace Transform : Solving 

Differential Equation 
//Y(S) = alpha /( s ( s+1) ( s+2) ) 

s = 7,s ; 
syms t ; 

alpha = 1; //Alpha value assigned as some constant 
one 

[A] = pf ss ( alpha/ ( s * ( s+1 ) * ( s +2) ) ) ; 

FI = i laplace ( A ( 1 ) , s , t ) 

F2 = ilaplace (A (2) , s , t) 

F3 = ilaplace (A (3) , s , t) 
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F = F1+F2+F3 
disp (F) 

/ / result 

// (— %e~ — t )+((%e~— (2*t) ) / 2)+(l/2 ) 

Scilab code Exa 9.37 

clear 

//clear// 

// Example9 . 3 7 : Un i lat er al Laplace Transform : Solving 

Differential Equation 
//Y(S) = alpha /( s ( s+1) ( s+2) ) 

s = °/ 0 s ; 
syms t ; 

alpha = 1; //Alpha value assigned as some constant 
one 

[A] = pf ss ( alpha/ ( s * ( s+1 ) * ( s +2) ) ) ; 

FI = i laplace ( A ( 1 ) , s , t ) 

F2 = ilaplace (A (2) , s , t) 

F3 = ilaplace (A (3) , s , t) 

F = F1+F2+F3 
disp (F) 

II result 

// (— %e~ — t )+((%e~— (2*t) ) / 2)+(l/2 ) 


Scilab code Exa 9.38 clear 

//clear// 

// Example9 . 3 8 : Un i lat er al Laplace Transform : Solving 
Differential Equation 

//Y(S)=[beta(s+3)/((s + l)(s+2))] + [gamma / ((s+2)(s+2)) 
] + [ alpha / ( s ( s + 1) (s+2) ) ] 

s = °/o s ; 
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syms t ; 

alpha = 2; //input constant 
beta_B = 3; //intial condition 
gamma.v = -5; //initial condition 

Y 1 = 1/s; 

Y2 = l/(s + l) ; 

Y3 = 3/ ( s + 2) ; 

Y = Y1-Y2+Y3; 
di sp ( Y) 

y = ilaplace (Y , s , t) 
disp (y) 

/ / result 

// ( — %e~( — t))+3*(%e~( — (2*t)))+l 

Scilab code Exa 9.38 

clear 

//clear// 

// Example9 . 3 8 : U n i 1 at er al Laplace Transform : Solving 
Differential Equation 

//Y(S ) = [beta ( s+3) / ( ( s + 1) (s+2) ) ] + [gamma/ ( ( s+2) (s+2) ) 
] + [ alpha / ( s ( s + 1) (s+2) ) ] 

s = °/o s ; 
syms t ; 

alpha = 2; //input constant 
beta_B = 3; //intial condition 
gamma_v = -5; //initial condition 

Y 1 = 1/s; 

Y2 = l/(s + l) ; 

Y3 = 3/ (s+2) ; 

Y = Y1-Y2+Y3; 
di sp ( Y) 

y = ilaplace (Y , s , t) 
disp (y) 

II result 

// ( — %e" ( — t ) ) +3* (%e "( — (2*t)))+l 


209 


Chapter 10 
The Z Transform 


210 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

2 

3 

4 

5 

6 

7 


Chapter 11 

The Z Transform 


Scilab code Exa 10.1 Scilab code Exa 10.1 Ztransform of x[n] 

//clear// 

// ExamplelO . 1 : Ztransform of x[n] = (a) A n.u[n] 

syms n z ; 
a = 0.5; 
x = ( a) ~ n 

X = symsum (x* ( z ~ ( -n) ) , n , 0 , °/ 0 inf ) 
disp (X , ” ans=” ) 

//Result 

/ / 1 .0* (2 " ( — %inf — 1) *z A ( — %inf — 1) — 1) / (1/ (2* z ) —1) 
//Equivalent to — l/(0.5*(z~ — 1)— 1) 

Ztransform of x[nj 

//clear// 

// ExamplelO . 1 : Ztransform of x[n] = (a)"n.u[n] 

syms n z ; 
a = 0.5; 
x = ( a) ~ n 

X = symsum (x* ( z ~ ( -n) ), n , 0 , %inf ) 
disp (X , ” ans=” ) 
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//Result 

// 1.0*(2" ( — %inf — 1) *z A ( — %inf — 1) — 1) / (1/ (2*z) — 1) 
//Equivalent to — l/(0.5*(z~— 1) — 1) 


Scilab code Exa 10.2 Scilab code Exa 10.2 Z transform of x[n] =-a n .u[— n— 

1 ] 


//clear// 

//Example 1 0 . 2 : Z transform of x[n] = — a"n. u[ — n — 1] 
//a = 0.5 

clear ; 
close ; 
clc ; 

syms n z ; 
a = 0.5; 
x=- (0 . 5) ~ (-n) 

X=symsum(x*(z“(n)) , n , 1 , %inf ) 
disp (X , ” ans=” ) 

/ / Result 

// — 1.0*(2''(%inf + l)*z~(%inf + l)— 2*z) / (2* z — 1) 
//Equivalent to — 1*— 2*z / (2 * z — 1) = 1/(1 —0.5* z ~ — 1) 


Scilab code Exa 10.3 Scilab code Exa 10.3 Z transform of x[n] 

//clear// 

//Example 10.3: Z transform of x[n] = 7. ( 1/3) "n . u [n 
] — 6.(l/2)"n.u[n] 
syms n z ; 


212 


4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

1 

2 

3 

4 


xl = (0 . 33) ~ (n) 

Xl = symsum(7*xl*(z''(-n)) ,n,0, 0 /oinf) 
x2 = ( 0 . 5) ~ (n) 

X2 = symsum(6*x2*(z''(-n)) ,n,0, 0 /,inf) 

X = X1-X2 
disp (X , ” ans=” ) 

//Result 

// — 6.0*(2 ~ ( — %inf— 1) *z A ( — %inf — 1)— 1) / (1/ (2 * z ) —1) 

//Equivalent to — 6* — 1/(0. 5* z ~ — 1 —1) 

//The Region of Convergence is | z | > 1 / 2 

Z transform of x[nj 

//clear// 

//Example 1 0 . 3 : Z transform of x[n] = 7 . ( 1 / 3 ) ~ n . u [ n 
] — 6.(l/2) A n.u[n] 
syms n z ; 
xl = (0 . 33) ~ ( n) 

XI = syrasura (7* xl * ( z ~ ( - n) ) , n , 0 , °/„ i n f ) 
x2 = ( 0 . 5) ~ (n) 

X2 = syrasura (6*x2* (z~ ( - n) ) , n , 0 , °/„ i n f ) 

X = X1-X2 
disp (X , ” ans=” ) 

//Result 

// — 6.0*(2~( — %inf— l)*z“(- %inf — 1)— 1) / (1/ (2 * z ) —1) 

//Equivalent to — 6* — 1/(0. 5*z~ — 1 —1) 

//The Region of Convergence is | z | > 1 / 2 


Scilab code Exa 10.4 Z-transform of sine signal 

//clear// 

// ExamplelO . 4 : Z— transform of sine signal 

syms n z ; 

Wo =°/„pi/4; 
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a = (0 . 33) “n ; 

xl = °/ 0 e~(sqrt(-l)*Wo*n) ; 
Xl=symsum(a*xl*(z~(-n)) , n , 0 , 0 /,inf ) 
x2 = %e ~ ( - sqrt (-l)*Wo*n) 

X2 = symsum(a*x2*(z''(-n)) , n , 0 , 0 /,inf ) 

X =(l/(2*sqrt (-1)))*(X1-X2) 
disp (X , ” ans=” ) 

Scilab code Exa 10.4 

Z-transform of sine signal 

//clear// 

// ExamplelO . 4 : Z— transform of sine signal 

syms n z ; 

Wo =°/„pi/4; 

a = (0 . 33) ~n ; 

x 1 = % e “ ( s qr t (-l)*Wo*n) ; 

XI = syrasum(a*xl*(z''(-n)) ,n,0, 0 / 0 inf) 
x2 = %e ~ ( - sqrt (-l)*Wo*n) 

X2 = syrasum(a*x2*(z''(-n)) ,n,0, 0 / 0 inf) 

X =(l/(2*sqrt (-1)))*(X1-X2) 
disp (X , ” ans=” ) 


Scilab code Exa 10.5 Z-transform of Impulse Sequence 

//clear// 

// ExamplelO . 5 : Z— transform of Impulse Sequence 
syms n z ; 

X = symsui(l*(z‘(-n)) ,n,0,0) ; 
disp (X , ” ans=” ) 

/ / Result 

// 1 
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Scilab code Exa 10.5 


Z-transform of Impulse Sequence 

//clear// 

// ExamplelO . 5 : Z— transform of Impulse Sequence 
syms n z; 

X = symsum(l*(z‘(-n)) , n , 0 , 0) ; 
disp (X , ” ans=” ) 

//Result 

// 1 


Scilab code Exa 10.6 Scilab code Exa 10.6 Z transform of x[n] Z trans- 
form of x[n] 

//clear// 

//Example 10.6:Z transform of x[n] = a 'n , 0 < n < N 
-1 

syms n z; 
a = 0.5; 

N =6; 
x=(a) ~(n) 

X=symsum(x*(z~(-n)) , n , 0 , N ) 
disp (X , ” ans=” ) 

//Result 

// 0. 5/ z +0.25/ z" 2 + 0. 125/ z '3 + 0. 0625/ z '4 + 0.03125/ z 
"5 + 0.015625 / z " 6 + 1 .0 

//clear// 

//Example 10.6: Z transform of x[n] = a " n , 0 < n < N 
-1 

syms n z; 
a = 0.5; 
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N =6; 
x=(a) ~(n) 

X=symsum(x*(z~(-n)) , n , 0 , N ) 
disp (X , ” ans=” ) 

/ / Result 

// 0. 5/ z +0.25/ z A 2 + 0. 125/ z '3 + 0. 0625/ z "4 + 0.03125/ z 
"5 + 0.015625 / z " 6 + 1 .0 


Scilab code Exa 10.7 Scilab code Exa 10.7 Z transform of x[n] 

//clear// 

//Example 1 0 . 7 : Z transform of x[n] = b " n . u [ n] + b" — n . u 
[-n-l] 
syms n z ; 
b = 0.5; 
xl=(b) ~(n) 
x2= (b) “ ( -n) 

Xl = symsum(xl*(z~(-n)),n,0, 0 /oinf) 

X2=symsum(x2*(z~(n)) , n , 1 , °/„inf ) 

X = X1+X2; 
disp (X , ” ans=” ) 

//Result 

/ / +1.0* (2 " ( — %inf — 1) *z " ( — %inf — 1) — 1) / (1/ (2*z) —1) 
//Equivalent to — 1/(0. 5*z"— 1 — 1) 

//Region of Convergence | z | >0.5 

Z transform of x[nj 

//clear// 

//Example 1 0 . 7 : Z transform of x[n] = b " n . u [ n] + b"— n . u 
[-n-l] 

syms n z ; 
b = 0.5; 
xl = (b) ~ (n) 
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x2=(b) ~(-n) 

Xl=symsum(xl*(z~(-n) ) ,n,0 , %inf ) 
X2=symsum(x2*(z~ (n) ) ,n, 1 , °/„inf ) 

X = X1+X2; 
disp (X , ” ans=” ) 

/ / Result 

/ / + 1.0* (2 * ( — %inf — 1) *z A ( — %inf — 1) — 1) / (1/ (2*z) — 1) 
//Equivalent to — 1/(0. 5*z'— 1 — 1) 

//Region of Convergence | z | >0.5 


Scilab code Exa 10.9 clear 

//clear// 

//ExamplelO . 9 : Inverse Z Transform :ROC | z | > 1 / 3 

z = °/ 0 z ; 

syms n zl;//To find out Inverse z transform z must 
be linear z = zl 

X =z*(3*z-(5/6) )/((z-(l/4))*(z-(l/3) )) 

XI = denom (X) ; 
zp = roots (XI ) ; 

XI = zl*(3*zl-(5/6))/((zl-(l/4))*(zl-(l/3))) 

FI = Xl*(zl ~ (n-1) ) *(zl-zp (1) ) ; 

F2 = Xl*(zl" (n-1) ) *(zl-zp (2) ) ; 

hi = limit (FI , zl , zp ( 1) ) ; 

disp (hi , ’hi [ n] = ’ ) 

h2 = limit (F2 , zl , zp (2) ) ; 

disp (h2 , ’ h2 [n]= ’ ) 

h = hi +h2 ; 

disp (h , ’h [n]= ’ ) 

// // Result 

//h[ii]= (1/4) ' n + (2/3) ~n 


Scilab code Exa 10.9 clear 
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//clear// 

//ExamplelO . 9 : Inverse Z Transform :ROC | z | > 1 / 3 

z = °/ 0 z ; 

syms n zl;//To find out Inverse z transform z must 
be linear z = zl 

X =z*(3*z-(5/6) )/ ((z-(l/4) ) *(z-(l/3) )) 

XI = denom (X) ; 
zp = roots (XI ) ; 

XI = zl*(3*zl-(5/6) )/ C(zl— (1/4) )*(zl-(l/3) )) 

FI = Xl*(zl~(n-l))*(zl-zp(l)) ; 

F2 = Xl*(zl~(n-l))*(zl-zp(2)) ; 

hi = limit (FI , zl , zp ( 1) ) ; 

disp (hi , ’hi [ n] = ’ ) 

h2 = limit (F2 , zl , zp (2) ) ; 

disp (h2 , ’ h2 [n]= ’ ) 

h = hi +h2 ; 

disp (h , ’h [n]= ’ ) 

// // Result 

//h[n]= (1/4) ~ n + (2/3) * n 


Scilab code Exa 10.10 Inverse Z Transform 

//clear// 

// ExamplelO . 1 0 : Inverse Z Transform :ROC l/4<|z| <1/3 

z = °/ 0 z ; 

syms n zl;//To find out Inverse z transform z must 
be linear z = zl 

X =z*(3*z-(5/6) )/((z-(l/4))*(z-(l/3) )) 

XI = denom (X) ; 
zp = roots (XI ) ; 

XI = zl*(3*zl -(5/6) ) / ((zl-(l/4) ) *(zl-(l/3) ) ) 

FI = Xl*(zl ~ (n-1) ) * (zl-zp (1) ) ; 

F2 = Xl*(zl “ (n-1) ) * (zl-zp (2) ) ; 
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hi = limit (FI , zl , zp ( 1) ) ; 
disp (hi * ’u(n) ’ , ’hi [n]= ’) 
h2 = limit (F2 , zl , zp (2) ) ; 
disp((h2)* ’u(-n-l) ’ , ’ h2 [n]= ’ ) 
disp((hl)*’u(n) ’ -(h2) * ’u(n-l) ’ , ’h[n]= ’) 

// // Result 

// h[n]= u(n) /4 A n— 2*u(n — l)/3 A n 

//Equivalent to h[n] = (1/4) "n . u [n] — 2* ( 1/3) ~n . u[ — n — 1] 

Scilab code Exa 10.10 

Inverse Z Transform 

//clear// 

// ExamplelO . 1 0 : Inverse Z Transform :ROC l/4<|z| <1/3 

z = °/ 0 z ; 

syms n zl;//To find out Inverse z transform z must 
be linear z = zl 

X =z*(3*z-(5/6) )/((z-(l/4) ) * (z- ( 1/3) )) 

XI = denom (X) ; 
zp = roots (XI ) ; 

XI = zl*(3*zl-(5/6) )/ ((zl— (1/4) )*(zl-(l/3) )) 

FI = Xl*(zl ~ (n-1) ) *(zl-zp (1) ) ; 

F2 = Xl*(zl~ (n-1) ) *(zl-zp (2) ) ; 

hi = limit (FI , zl , zp ( 1) ) ; 

disp (hi * ’u(n) ’ , ’hi [n]= ’) 

h2 = limit (F2 , zl , zp (2) ) ; 

disp ( (h2) * ’u(-n-l) ’ , ’ h2 [n]= ’ ) 

disp ( (hi) * ’u(n) ’ - ( h2 ) * ’ u ( n-1) ’ , ’h[n] = ’) 

// // Result 

// h[n]= u(n) /4 A n— 2*u(n — l)/3 A n 

//Equivalent to h[n] =(1/4) A n . u [ n] — 2* ( 1/3) * n . u[ — n — 1] 


Scilab code Exa 10.11 Inverse Z Transform 
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//clear// 

// ExamplelO . 1 1 : Inverse Z Transform :ROC | z | < 1 / 4 

Z = °/ 0 z ; 

syms n zl;//To find out Inverse z transform z must 
be linear z = zl 

X =z*(3*z-(5/6) )/((z-(l/4))*(z-(l/3) )) 

XI = denom (X) ; 
zp = roots (XI ) ; 

XI = zl*(3*zl-(5/6) )/ C(zl— (1/4) )*(zl-(l/3) )) 

FI = Xl*(zl~(n-l))*(zl-zp(l)) ; 

F2 = Xl*(zl"(n-l))*(zl-zp(2)) ; 
hi = limit (FI , zl , zp ( 1) ) ; 
disp (hi* ’u(-n-l) ’ , ’hi [n]= ’ ) 
h2 = limit (F2 , zl , zp (2) ) ; 
disp ( (h2) * ’u(-n-l) ’ , ’ h2 [n]= ’ ) 

di sp ( - ( hi ) * ’u(-n-l) ’-(h2)* ’u(-n-l) ’ , ’h[n] = ’) 

// // Result 

// h[n]= — u( — n — l)/4~n — 2*u( — n — l)/3~n 

//Equivalent to h [n] = — (l/4)*n.u[— n — 1]— 2*(l/3)"n.u[— 
n — 1] 


Scilab code Exa 10.11 Inverse Z Transform 

//clear// 

//ExamplelO . 1 1 : Inverse Z Transform :ROC | z | < 1 / 4 

z = °/ 0 z ; 

syms n zl;//To find out Inverse z transform z must 
be linear z = zl 

X =z*(3*z-(5/6) )/ ((z-(l/4) ) *(z-(l/3) )) 

XI = denom (X) ; 
zp = roots (XI ) ; 

XI = zl*(3*zl -(5/6) ) / ((zl-(l/4) ) *(zl-(l/3) ) ) 

FI = Xl*(zl ~ (n-1) ) *(zl-zp (1) ) ; 

F2 = Xl*(zl~ (n-1) ) *(zl-zp (2) ) ; 
hi = limit (FI , zl , zp ( 1) ) ; 
disp (hi* ’u(-n-l) ’ , ’hi [n]= ’ ) 
h2 = limit (F2 , zl , zp (2) ) ; 
disp ( (h2) * ’u(-n-l) ’ , ’ h2 [n]= ’ ) 
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disp ( — (hi) * ’u(-n-l) ’-(h2)* ’u(-n-l) ’ , ’h[n] = ’) 

// // Result 

// h[n]= — u(— n — l)/4 A n- 2*u(— n — l)/3*n 

//Equivalent to h[n] = — (1/4) ~ n . u[ — n — 1] — 2* ( 1/3) " n . u[ — 
n — 1] 


Scilab code Exa 10.12 Inverse z tranform 

//clear// 

// ExamplelO . 1 2 : Inverse z tranform : For Finite 
duration discrete sequence 

syms z ; 

X = [4* z ~ 2 0 2 3*z ~ -1] ; 

n = -2:1; 

for i = l:length(X) 
x(i) = X(i)*(z“n(i)) ; 
end 

di sp (x , ’ x [ n ] = ’ ) 

Scilab code Exa 10.12 

Inverse z tranform 

//clear// 

// ExamplelO . 1 2 : Inverse z tranform : For Finite 
duration discrete sequence 

syms z ; 

X = [4* z ~ 2 0 2 3*z ~ -1] ; 

n = -2:1; 

for i = l:length(X) 

x ( i ) = X(i)*(z“n(i)) ; 
end 

di sp (x , ’ x [ n] = ’ ) 
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Scilab code Exa 10.13 Inverse z tranform of InFinite duration discrete 
sequence 

//clear// 

// ExamplelO . 1 3 : Inverse z tranform oflnFinite 
duration discrete sequence 
//Power Series Method (OR) //Long Division Method 
z = °/ 0 z ; 
a = 2 ; 

X = ldiv (z , z-a , 5) 

Scilab code Exa 10.13 

Inverse z tranform of InFinite duration discrete sequence 

//clear// 

// ExamplelO . 1 3 : Inverse z tranform oflnFinite 
duration discrete sequence 
//Power Series Method (OR) //Long Division Method 
z = °/ 0 z ; 
a = 2 ; 

X = ldiv (z, z-a, 5) 


Scilab code Exa 10.18 Ztransform-Differentiation Property 

//clear// 

// ExamplelO . 1 8 : Ztransform — D i f f e re n t i at i o n Property 
/ / x [ n ] = ( a ) A n . u [ n ] 
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syms n z ; 
a = 0.5; 
x = ( a) ~ n 

X = symsum (x* ( z ~ ( -n) ) , n , 0 , %inf ) 

XI = -l/((l/(2*z))-l) / / z transform of 0 . 5 ~ n . u [ 

11 ] 

Y = -z*diff(X,z) //Differentiation property of z— 
transform 

disp (X , ” ans=” ) 
disp (Y , ” ans=” ) 

//Result 

//X(z) = 1 . 0* (2 " ( — %inf— 1) *z ~ ( — %inf — 1)— l)/(l/(2*z)— 1) 
//Y(z) = — 1.0*( — %inf — 1) *2 ~ ( — %inf — 1) *z ~ ( — %inf — 1) 

/ ( 1 / ( 2* z ) — 1) 

/ / Y1 ( z ) = l/(2*(l/(2*z) -1) " 2* z ) 

//Equivalent to Y 1 ( z ) = 0 . 5 * z " — 1 / ( ( 1 — 0 . 5 * z ~ — 1 ) ~ 2 ) 

Scilab code Exa 10.18 

Ztransform-Differentiation Property 

//clear// 

// ExamplelO . 1 8 : Ztransfor m — D i f f e r e n t i a t i o n Property 
/ / x [ n ] = ( a ) ~ n . u [ n ] 

syms n z ; 
a = 0.5; 
x = ( a) ~ n 

X = symsum (x* ( z ~ ( -n) ), n , 0 , %inf ) 

XI = — 1/( (l/(2*z) ) — 1) / / z transform of 0 . 5 ~ n . u [ 

11 ] 

Y = -z*diff(X,z) //Differentiation property of z— 
transform 

disp (X , ” ans=” ) 
disp (Y , ” ans=” ) 

/ / Result 

//X(z) = 1 . 0* (2 A ^ %inf — 1) *z "( — %inf — 1) — 1) / ( 1/ (2* z ) — 1) 
//Y(z) = — 1.0*( — %inf — 1) * 2 " ( — %inf — 1) *z " ( — %inf — 1) 
/(l/(2*z) —1) 
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//Yl(z) = 1/ ( 2 * ( 1 / ( 2 * z ) — l)"2*z) 

//Equivalent to Y 1 ( z ) = 0 . 5 * z " — 1 / ( ( 1 — 0 . 5 * z ~ — 1 ) " 2 ) 


Scilab code Exa 10.19 Z Transform : Initial Value Theorem 

//clear// 

// ExamplelO . 1 9 : Z Transform : Initial Value Theorem 

z = °/ 0 z ; 

syms n zl;//To find out Inverse z transform z must 
be linear z = zl 

X =z*(z-(3/2) ) / ((z— (1/3) ) *(z-(l/2) ) ) 

XI = denom (X) ; 
zp = roots (XI ) ; 

XI = zl*(zl-(3/2) )/ C(zl— (1/3) ) *(zl-(l/2) )) 

FI = Xl*(zl ~ (n-1) ) *(zl-zp (1) ) ; 

F2 = Xl*(zl“ (n-1) ) *(zl-zp (2) ) ; 

xl = limit (FI , zl , zp ( 1) ) ; 

x2 = limit (F2 , zl , zp (2) ) ; 

x = x 1 +x2 ; 

di sp (x , ’x[n]= ’ ) 

x_initial = limit (x , n , 0) ; 

di sp ( x_ initial , ’ x [0] = ’ ) 

// // Result 

//x[n]= 7/3 " n — 3*2 A (1 — n) 

/ / x [ 0 ] = 1; Initial Value 

Scilab code Exa 10.19 

Z Transform : Initial Value Theorem 

//clear// 

// ExamplelO . 1 9 : Z Transform : Initial Value Theorem 

z = °/ 0 z ; 
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syms n zl;//To find out Inverse z transform z must 
be linear z = zl 

X =z*(z-(3/2) ) / ((z— (1/3) ) *(z-(l/2) ) ) 

XI = denom (X) ; 
zp = roots (XI ) ; 

XI = zl*(zl-(3/2) )/ C(zl— (1/3) )*(zl-(l/2) ) ) 

FI = Xl*(zl ~ (n-1) ) * (zl-zp (1) ) ; 

F2 = Xl*(zl~ (n-1) ) *(zl-zp (2) ) ; 

xl = limit (FI , zl , zp ( 1) ) ; 

x2 = limit (F2 , zl , zp (2) ) ; 

x = x 1 +x2 ; 

di sp (x , ’x[n]= ’ ) 

x_initial = limit (x , n , 0) ; 

di sp ( x_ initial , ’ x [0] = ’ ) 

// // Result 

//x[n] = 7/3 A n — 3*2 A (l— n) 

/ / x [ 0 ] = 1; Initial Value 


Scilab code Exa 10.23 Inverse Z Transform H(z) =z/z-a 

//clear// 

// ExamplelO . 2 3 : Inverse Z Transform H(z) =z/z— a 
/ / z = %z ; 

syms n z ; 
a = 2 ; 

H = z/ (z-a) ; 

F = H*z “ (n-1 ) * (z-a) ; 
h = limit (F , z , a) ; 
disp (h , ’h [ n] = ’ ) 

Scilab code Exa 10.23 

Inverse Z Transform H(z) =z/z-a 
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//clear// 

// ExamplelO . 2 3 : Inverse Z Transform H(z) =z/z— a 
/ / z = %z ; 

syms n z ; 
a = 2 ; 

H = z/(z-a) ; 

F = H*z “ (n-1 ) * (z-a) ; 
h = limit (F , z , a) ; 
disp (h , ’h [n]= ’ ) 


Scilab code Exa 10.25 Coefficient Difference equations 

//clear// 

// ExamplelO . 2 5 : LTi Systems characterized by Linear 
Constant 

//Coefficient Difference equations 
//Inverse Z Transform 
/ / z = %z ; 

syms n z ; 

HI = z/(z-(l/2) ) ; 

H2 = ( 1/3) / (z- ( 1/2) ) ; 

FI = Hl*z~(n-1) *(z-(l/2) ) ; 

F2 = H2*z~(n-1) *(z-(l/2) ) ; 
hi = limit (FI , z , 1/2) ; 
disp (hi , ’hi [ n] = ’ ) 
h2 = limit (F2 , z , 1/2) ; 
disp (h2 , ’ h2 [n]= ’ ) 
h = hi +h2 ; 
disp (h , ’h [n]= ’ ) 

//Result 

/ / h [ n ] = [(l/2)Dr] + [2'(l-n)]/3 

//Which is Equivalent to h[n] = [(l/2)*n] + [(l/2)*(n 

-l)]/3 
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Scilab code Exa 10.25 Coefficient Difference equations 

//clear// 

// ExamplelO . 2 5 : LTi Systems characterized by Linear 
Constant 

//Coefficient Difference equations 
//Inverse Z Transform 
/ / z = %z ; 

syms n z ; 

HI = z/(z-(l/2) ) ; 

H2 = ( 1/3) / (z- (1/2) ) ; 

FI = Hl*z~(n-1) *(z-(l/2) ) ; 

F2 = H2*z~(n-1) *(z-(l/2) ) ; 
hi = limit (FI ,z , 1/2) ; 
disp (hi , ’hi [ n] = ’ ) 
h2 = limit (F2 , z , 1/2) ; 
disp (h2 , ’ h2 [n]= ’ ) 
h = hi +h2 ; 
disp (h , ’h [n]= ’ ) 

//Result 

/ / h [ n ] = [(1/2) " n ] + [2 A (l— n)]/3 

//Which is Equivalent to h[n] = [(l/2)*n] + [(l/2)*(n 
-l)]/3 


Scilab code Exa 10.33 Differentiation Property of Unilateral Ztransform 

//clear// 

// ExamplelO . 3 3 : D i f f e r e n t i at i on Property of 
Unilateral Ztransform 
// x [n] = ( a ) ~ ( n + 1) . u [ n + 1] 

syms n z ; 
a = 0.5; 
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6 x = ( a) ~ ( n + 1 ) 

7 X = symsum(x*(z~(-n)),n,-l, 0 / 0 inf) 

8 disp (X , ” ans=” ) 

9 //Result 

10 //X(z)= 0 . 5* (2 A ( — %inf — 1) *z A ( — %inf — 1) — 2*z ) / ( 1 / (2 * z ) 

-1) 

11 //Equivalent to z /(I — 0.5* z ~ — 1) 

Scilab code Exa 10.33 

Differentiation Property of Unilateral Ztransform 

1 //clear// 

2 // ExamplelO . 3 3 : D i f f e r e n t i at i on Property of 

Unilateral Ztransform 

3 // x [ n ] = ( a) ~ (n + 1) . u [ n + 1] 

4 syms n z ; 

5 a = 0.5; 

6 x = ( a) " (n+1 ) 

7 X = symsum(x*(z“(-n)),n,-l, 0 /oinf) 

8 disp (X , ” ans=” ) 

9 //Result 

10 //X(z)= 0 . 5* (2 A ( — %inf — 1) *z A ( — %inf — 1) — 2*z ) / ( 1 / (2 * z ) 

-1) 

11 //Equivalent to z /(I — 0.5*z ~ — 1) 


Scilab code Exa 10.34 Unilateral Ztransform- partial fraction 

1 //clear// 

2 // ExamplelO . 34 : Unilateral Ztransform— partial 

fraction 

3 // X(z) =(3 (5/6) *(z A — 1) ) / ((1 — (1/4) *(z A — 1)) * ( 1 — (1/3) 

*( z ~-!)) ) 

4 z = y o z ; 
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s = y„s ; 
syms n t ; 
a = 0.5; 

[A]=pfss((3-(5/6)*(z~-l))/((l-(l/4)*(z~-l))*(l-(l/3) 
*(z~-l)))) 
xl = horner (A (1) , z) 
x2 = horner (A (2) , z) 
x3 = A (3) 
x = xl+x2+x3 
disp (xl , ” ans=” ) 
disp (x2 , ” ans=” ) 
disp (x3 , ” ans=” ) 
di sp (x , ” ans=” ) 

/ / Result 

// 0.6666667 


// 

// - 0.3333333 + z 


// 


0.25 


// 

// - 0.25 + z 


//3 


/ / sum 

// 


of these , gives the original value 


2 

— 0.8333333 z + 3z 


// 


// 


// 


2 


// 0.0833333 - 0.5833333z + z 


Scilab code Exa 10.34 


Unilateral Ztransform- partial fraction 


//clear// 
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// ExamplelO . 34 : Unilateral Ztransform— partial 
fraction 

// X(z) =(3 — (5/6) * ( z A — 1) ) / ((1 — (1/4) *(z~— 1)) * ( 1 — (1/3) 
*(z~-l))) 

Z = °/ 0 z ; 
s = °/ 0 s ; 
syms n t ; 
a = 0.5; 

[A]=pfss((3-(5/6)*(z~-l))/((l-(l/4)*(z~-l))*(l-(l/3) 
*(z~-l)))) 
xl = horner (A (1) , z) 
x2 = horner (A (2) , z) 
x3 = A (3) 
x = x 1 + x2 + x3 
disp (xl , ” ans=” ) 
disp (x2 , ” ans=" ) 
disp (x3 , ” ans=” ) 
di sp (x , ” ans=" ) 

//Result 

// 0.6666667 

// 

// - 0.3333333 + z 

// 0.25 

// 

// - 0.25 + z 

// 3 

/ / sum of these, gives the original value 
// 2 
// - 0.8333333 z + 3z 

// 

// 2 

// 0.0833333 - 0.5833333z + z 
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Scilab code Exa 10.36 Output response of an LTI System 

//clear// 

//Example 10.36:To find output response of an LTI 
System 

syms n z ; 

H = z/(z+3) 

X = z/(z-l) 

Y = X*H 

FI = Y*(z~(n-1) ) * (z-1) ; 
yl = limit (FI ,z,l) ; 

F2 = Y* (z~ (n-1) ) *( z +3) ; 
y2 = limit (F2 , z , -3) ; 
disp(yl*”u(n)”+y2*”u(n)” , ’ y [ n ] = ’ ) 

/ / Result 

/ / y [ n ] = u(n)/4-(-3) A (n + 1) *u (n) /4 
//Equivalent to = (1/4). u[n]— (3/4) ( — 3 ) " n . u [ n ] 

Scilab code Exa 10.36 

Output response of an LTI System 

//clear// 

//Example 10.36:To find output response of an LTI 
System 

syms n z ; 

H = z/(z+3) 

X = z/(z-l) 

Y = X*H 

FI = Y* (z~ (n-1) )*(z-l) ; 
yl = limit (FI , z , 1) ; 

F2 = Y* (z~ (n-1) ) *(z+3) ; 
y2 = limit (F2 , z , -3) ; 
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d i s p ( y 1 * ” u ( n ) ” + y 2 * ” u ( n ) ” , ’ y [n]= ’ ) 

/ / Result 

/ / y [ 11 ] = u ( n ) / 4 — ( — 3 ) *(n + l)*u(n)/4 
//Equivalent to = (1/4). u[n]— (3/4) ( — 3 ) A n . u [ n ] 


Scilab code Exa 10.37 Output response of an LTI System 

//clear// 

//Example 10.37:To find output response of an LTI 
System 

syms n z ; 

alpha = 8; //input constant 

beta_b = 1; //initial condition y [ — 1 ] = 1 

Y1 = - ( (3*beta_b*z) / (z+3) ) 

Y2 = (alpha*z~2/ ( (z+3) *(z-l) ) ) 

FI = Yl*(z~(n-1) ) * (z+3) ; 
yl = limit (FI , z , -3) ; 

F2 = Y2* (z~ (n-1) )*(z+3) ; 
y2 = limit (F2 , z , -3) ; 

F3 = Y2*(z~ (n-1) ) *(z-l) ; 

y3 = limit (F3 , z , 1) ; 

disp((yl + y2 + y3)* ’ u (n) ’ , ’ y [n]= ’ ) 

/ / Result 

/ / y [ 11 ] = (2-(-3) A (n+l))*u(n) 

Scilab code Exa 10.37 

Output response of an LTI System 

//clear// 

//Example 10.37:To find output response of an LTI 
System 

syms n z ; 

alpha = 8; //input constant 
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beta_b = 1; //initial condition y [ — 1 ] = 1 

Y1 = - ( (3*beta_b*z) / (z+3) ) 

Y2 = (alpha*z~2/ ( (z+3) *(z-l) ) ) 

FI = Yl*( z ~ (n-1) ) *( z +3) ; 
yl = limit (FI , z , -3) ; 

F2 = Y2*(z~(n-l))*(z + 3) ; 
y2 = limit (F2 , z , -3) ; 

F3 = Y2* (z~ (n-1) )*(z-l) ; 

y3 = limit (F3 , z , 1) ; 

disp ( (yl+y2 + y3) * ’u(n) ’ , ’ y [n]= ’ ) 

//Result 

/ / y [ n ] = (2 — ( — 3) A ( n + 1) ) *u (n) 
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Chapter 13 

Linear Feedback Systems 


Scilab code Exa 11.1 Root locus Analysis of Linear Feedback Systems 

//clear// 

// Examplell . 1 : Root locus Analysis of Linear Feedback 
Systems 

//Continuous Time Systems 

//Refer figure 11.12(a) in Openhiem &Willksy page 
840 

s = y„s ; 

H = sys 1 in ( ’ c ’ , [ 1/ ( s + 1 ) ] ) ; 

G = sys 1 in ( ’ c ’ , [ 1/ ( s + 2) ] ) ; 

F = G*H ; 
elf ; 

evans (F , 3) 

Scilab code Exa 11.1 

Root locus Analysis of Linear Feedback Systems 

//clear// 

// Examplell . 1 : Root locus Analysis of Linear Feedback 
Systems 


235 


3 //Continuous Time Systems 

4 //Refer figure 11.12(a) in Openhiem &Willksy page 

840 

5 s = '/,s; 

6 H = syslin ( ’ c ’ , [1/ ( s + 1) ] ) ; 

7 G = syslin (’ c [1/ ( s + 2) ]) ; 

8 F = G*H ; 

9 elf ; 

10 evans (F , 3) 


Scilab code Exa 11.2 Continuous Time Systems 

1 //clear// 

2 // Examplell . 2 : Root locus Analysis of Linear Feedback 

Systems 

3 //Continuous Time Systems 

4 //Refer figure 11.14(a) in Openhiem &Willksy page 

844 

5 s = '/,s; 

6 G = syslin (’ c ’,[( s - 1 )/(( s + 1 )*( s +2) )]) ; 

7 elf ; 

8 evans (G , 2) 

Scilab code Exa 11.2 

Continuous Time Systems 

1 //clear// 

2 // Examplell . 2 : Root locus Analysis of Linear Feedback 

Systems 

3 //Continuous Time Systems 

4 //Refer figure 11.14(a) in Openhiem &Willksy page 

844 

5 s = °/ 0 s ; 
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6 G = sysl in ( ’ c ’ , [ ( s -1 ) / ( ( s+1 ) * ( s+2) ) ] ) ; 

7 elf ; 

8 evans (G , 2) 


Scilab code Exa 11.3 Discrete time system 

1 //clear// 

2 // Examplell . 3 : Root locus Analysis of Linear Feedback 

Systems 

3 ////Discrete time system 

4 //Refer figure 11.16(a) in Openhiem &Willksy page 

846 

5 z = °/ 0 z ; 

6 G = syslin ( ’d ’ , [z/ ( (z-0 . 5) * (z-0 . 25) ) ] ) ; 

7 elf ; 

8 evans (G , 2) 

Scilab code Exa 11.3 

Discrete time system 

1 //clear// 

2 // Examplell . 3 : Root locus Analysis of Linear Feedback 

Systems 

3 ////Discrete time system 

4 //Refer figure 11.16(a) in Openhiem &Willksy page 

846 

5 z = °/,z ; 

6 G = syslin ( ’d ’ , [z/ ( (z-0 . 5) * (z-0 . 25) ) ] ) ; 

7 elf ; 

8 evans (G , 2) 
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Scilab code Exa 11.05 Nyquist criterion for Continuous Time Systems 

//clear// 

//Example 11.5: Nyquist criterion for Continuous Time 
Systems 

//Nyquist Plot 
s = y„s ; 

//Open Loop Transfer Function 

G = syslin (’ c [1/ ( s + 1) ]) ; 

H = syslin ( ’c [1/(0 . 5*s + l) ]) ; 

F = G*H ; 
elf ; 

nyquist (F) 

show_margins (F , ’nyquist ’) 

Scilab code Exa 11.05 

Nyquist criterion for Continuous Time Systems 

//clear// 

//Example 11.5: Nyquist criterion for Continuous Time 
Systems 

//Nyquist Plot 
s = °/ 0 s ; 

//Open Loop Transfer Function 
G = syslin (’ c [1/ ( s + 1) ]) ; 

H = syslin ( ’c [1/(0 . 5*s + l) ]) ; 

F = G*H ; 
elf ; 

nyquist (F) 

show_margins (F , ’nyquist ’) 
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Scilab code Exa 11.5 Bode Plot 


// clear // 

//Example 11.5:Nyquist criterion for Continuous Time 
Systems 
//Bode Plot 
s = °/ 0 s ; 

//Open Loop Transfer Function 

G = syslin (’ c [1/ ( s + 1) ]) ; 

H = syslin (’ c [1/ (0 . 5*s + l) ]) ; 

F = G*H ; 

elf ; 

bode (F , 0 . 01 , 100) 
show_raargins (F) 

Scilab code Exa 11.5 

Bode Plot 

//clear// 

//Example 11.5:Nyquist criterion for Continuous Time 
Systems 
//Bode Plot 
s = °/ 0 s ; 

//Open Loop Transfer Function 
G = syslin (’ c [1/ ( s + 1) ]) ; 

H = syslin ( ’c [1/(0 . 5*s + l) ]) ; 

F = G*H ; 

elf ; 

bode (F , 0 . 01 , 100) 
show_margins (F) 
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Scilab code Exa 11.6 Nyquist Plot 


//clear// 

//Example 11.6: Nyquist criterion for Continuous Time 
Systems 

//Nyquist Plot 
s = y„s ; 

//Open Loop Transfer Function 

F = syslinC ’c ’,[( s+1) / ( (s-1) * (0 . 5* s+1) ) ] ) 
elf ; 

nyquist (F) 

show_margins (F , ’nyquist ’) 


Scilab code Exa 11.6 Nyquist Plot 

//clear// 

//Example 11.6: Nyquist criterion for Continuous Time 
Systems 

//Nyquist Plot 
s = °/ 0 s ; 

//Open Loop Transfer Function 

F = syslinC ’c \ [( s + 1) / ( (s-1) * (0 . 5* s + 1) ) ] ) 
elf ; 

nyquist (F) 

show_margins (F , ’nyquist ’) 


Scilab code Exa 11.7 Nyquist Plot 

//clear// 

//Example 11.7 
//Nyquist Plot 
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s = y„s ; 

T =1; 

//Open Loop Transfer Function 

G = syslin ( ’ c ’ , [-°/ 0 e ~ ( -s *T) ] ) ; 
elf ; 

nyquist (G) 

show_margins (G , ’nyquist ’) 

Scilab code Exa 11.7 

Nyquist Plot 

//clear// 

//Example 11.7 
//Nyquist Plot 
s = °/ 0 s ; 

T =1; 

//Open Loop Transfer Function 
G = syslin (’ c [-°/ 0 e ~( -s*T) ]) ; 
elf ; 

nyquist (G) 

show_margins (G , ’nyquist ’) 


Scilab code Exa 11.8 Nyquist Plot 

//clear// 

//Example 11.8: Nyquist criterion for Discrete Time 
Systems 

//Nyquist Plot 
//Discrete Time System 
z = °/ 0 z ; 

//Open Loop Transfer Function 

F = syslin( ’d’ , [1/ (z* (z+0 . 5) ) ] ) 
elf ; 
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nyquist (F) 

show_margins (F , ’nyquist ’) 


Scilab code Exa 11.8 

Nyquist Plot 

//clear// 

//Example 11.8: Nyquist criterion for Discrete Time 
Systems 

//Nyquist Plot 
//Discrete Time System 
z = °/ 0 z ; 

//Open Loop Transfer Function 

F = syslin( ’d’ , [1/ (z* (z + 0 . 5) ) ] ) 
elf ; 

nyquist (F) 

show_margins (F , ’nyquist ’) 


Scilab code Exa 11.09 Root locus analysis of Linear feedback systems 

//clear// 

// Figure 1 1 . 9 : Root locus analysis of Linear feedback 
systems 

s = °/ 0 s ; 
beta_bl = 1 ; 
beta_b2 = -1; 

G1 = sysl in ( ’ c ’ , [2* bet a_b 1 /s ] ) ; 

G2 = syslin ( ’ c ’ , [2*beta_b2/s] ) ; 

H = sy s 1 in ( ’ c ’ , [s / ( s -2) ] ) ; 

FI = G1 *H ; 

F2 = G2 *H ; 
elf ; 

evans (FI ,2) 
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f igure 
evans (F2 ,2) 


Scilab code Exa 11.09 

Root locus analysis of Linear feedback systems 

//clear// 

// Figure 1 1 . 9 : Root locus analysis of Linear feedback 
systems 

s = °/ 0 s ; 
beta_bl = 1 ; 
beta_b2 = -1; 

G1 = sy si in ( ’ c ’ , [2* bet a_b 1 /s ] ) ; 

G2 = syslin ( ’ c ’ , [2*beta_b2/s] ) ; 

H = syslin (’ c [s/ (s-2) ]) ; 

FI = G 1 * H ; 

F2 = G2 *H ; 
elf ; 

evans (FI ,2) 
f igure 
evans (F2 ,2) 


Scilab code Exa 11.9 Gain and Phase Margins 

//clear// 

//Example 11.9: Gain and Phase Margins and their 
//associated cross over frequencies 

s =poly (0 , ’ s 7 ) ; // Define ss as polynomial variable 
//Create s transfer function in forward path 

F = sy s 1 in ( ’ c ’ , [ (4* ( 1 + 0 . 5* s ) ) / ( s * ( 1 + 2* s ) * ( 1 + 0 . 05* s 
+ (0 . 1 25* s ) ~ 2 ) ) ] ) 

B = syslinC ’c ’ , (l + 0*s)/ (l + 0*s)) 

0L = F*B ; 
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frain = 0.01; // Min freq in Hz 

fmax = 10; // Max freq in Hz 

scf ( 1) ; 

// elf ; 

// Plot frequency response of open loop transfer 
function 

bode (0L ,0.01,10); 

// display gain and phase margin and cross over 
frequencies 

show_margins (0L) ; 

[gm,frl] = g_margin(0L) 

[phm,fr2] = p_margin(0L) 
dispCgm, ’gain margin in dB ’ ) 

disp(frl,’gain cross over frequency in Hz’) 

disp(phm, ’phase margin in dB ’ ) 

disp(fr2, ’phase cross over frequency in Hz’) 


Scilab code Exa 11.9 Gain and Phase Margins 

//clear// 

//Example 11.9: Gain and Phase Margins and their 
//associated cross over frequencies 

s =poly (0 , ’ s ’ ) ; // Define ss as polynomial variable 
//Create s transfer function in forward path 

F = sys 1 in ( ’ c ’ , [ (4* ( 1 + 0 . 5* s ) ) / ( s * ( 1 + 2* s ) * ( 1 + 0 . 05* s 
+ (0 . 1 25* s ) “2) ) ] ) 

B = syslin ( ’ c ’ , (l + 0*s) / ( 1 + 0* s) ) 

0L = F*B ; 

frain = 0.01; // Min freq in Hz 

fra ax = 10; // Max freq in Hz 

scf (1) ; 

// elf ; 

// Plot frequency response of open loop transfer 
function 

bode (0L ,0.01,10); 

// display gain and phase margin and cross over 
frequencies 

show_margins (0L) ; 
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17 [gm,frl] = g_margin(OL) 

18 [phm,fr2] = p_margin(OL) 

19 disp(gm, ’gain margin in dB ’ ) 

20 disp(frl,’gain cross over frequency in Hz’) 

21 disp(phm, ’phase margin in dB ’ ) 

22 disp(fr2, ’phase cross over frequency in Hz’) 
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